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I.  Summary  of  Work 

We  developed  a  new  algorithm  to  adaptively  optimize  binary  and  quaternary 
signatures  for  code-division  multiple-access  (CDMA)  communications  over  mul¬ 
tipath  channels  and/or  asynchronous  channel  access.  Using  the  observed  signal 
autocorrelation  matrix,  the  algorithm  attempts  to  maximize,  over  the  binary  or 
quaternary  antipodal  sequence  field,  the  signal-to-interference-plus-noise  ratio 
(SINR)  at  the  output  of  the  maximum  SINR  linear  filter.  While  this  maximiza¬ 
tion  problem  is  NP-hard,  the  algorithm  is  seen  to  produce  in  short  polynomial 
time  highly  desirable  solutions  that  approach  in  performance  the  theory-only 
complex/real-field  optimal  signature  vectors.  Signature  adaptation  may  be  car¬ 
ried  out  in  either  a  single  or  a  multi-user  mode.  Simulation  studies  offer  direct 
performance  comparisons  with  other  known  binary  signature  set  designs  and 
the  theoretical  complex/real-valued  optimal  vectors. 

We  derived  new  bounds  on  the  periodic  (cyclic)  total  squared  correlation 
(PTSC)  of  binary  antipodal  signature  sets  for  any  number  of  signatures  K  and 
any  signature  length  L.  Optimal  designs  that  achieve  the  new  bounds  are  then 
developed  for  several  (K,  L)  cases.  As  an  example  that  arguably  may  be  against 
common  expectation,  it  is  seen  that  neither  the  Gold  nor  the  Kasami  sets  are 
PTSC  optimal. 

We  developed  a  new  method  for  the  optimization  of  binary  spreading  codes 
under  a  rank-2  approximation  of  the  inverse  interference  autocovariance  matrix 
where  the  rank-2-optimal  binary  code  is  obtained  in  lower  than  quadratic  com¬ 
plexity.  Significant  SINR  performance  improvement  was  demonstrated  over  the 
common  binary  hard-limited  eigenvector  design  which  was  seen  to  be  equivalent 
to  the  rank-l-optimal  solution. 

Our  recent  advances  in  the  area  of  binary  sequences  for  code-division  multi¬ 
plexing  provide  us  with  minimum  total-squared-correlation  (TSC)  optimal  sig¬ 
nature  sets  for  (almost)  all  signature  lengths  L  and  set  sizes  K.  The  sets  are 
scalable  as  long  as  K  <  L  (underloaded  systems)  and  non-scalable  when  K  >  L 
(overloaded  systems);  in  general  the  non-scalable  case  requires  signature  re- 
design/re-assignment  as  users  enter  or  exit.  We  derive  new  lower  bounds  on 
the  conditional  TSC  of  overloaded  binary  signature  sets  built  on  fixed  full-load 
TSC-optimal  sets.  Overloading  is  allowed  to  be  as  high  as  100%.  Scalable 
designs  that  achieve  the  new  bounds  are  then  developed.  To  evaluate  the  per¬ 
formance  of  the  proposed  designs,  we  compared  the  TSC  of  our  constructions 
to  the  unconditionally  minimum  achievable  TSC  values. 

We  derived  new  bounds  on  the  aperiodic  total  squared  correlation  (ATSC)  of 
binary  antipodal  signature  sets  for  any  number  of  signatures  K  and  any  signa¬ 
ture  length  L.  We  then  presented  optimal  designs  that  achieve  the  new  bounds 
for  several  ( K ,  L)  cases.  As  an  example,  we  proved  that  the  familiar  Gold  and 
(small  or  large)  Kasami  designs  are  not  ATSC-optimal  in  general.  The  optimal 
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signature  set  designs  that  we  provide  are,  in  this  sense,  better  suited  for  asyn¬ 
chronous  and/or  multipath  code-division  multiplexing  applications. 

We  developed  a  binary  signature  design  procedure  to  scale  upwards  over¬ 
loaded  minimum  total-squared-correlation  (TSC)  binary  signature  sets.  The 
quality  of  the  design  is  measured  against  the  recently  published  binary  TSC 
bounds. 

We  revisited  the  classical  problem  of  detecting  a  complex  signal  of  unknown 
amplitude  in  colored  Gaussian  noise  in  the  context  of  adaptive  detection  with 
limited  training  data  via  the  auxiliary-vector  (AV)  filter  estimation  algorithm. 
Based  on  statistical  conditional  optimization  criteria,  the  iterative  AV  algo¬ 
rithm  starts  from  the  target  vector  and,  adding  non-orthogonal  auxiliary  vector 
components,  generates  an  infinite  sequence  of  tests  that  converges  to  the  ideal 
matched  filter  (MF )  processor  for  any  positive  definite  input  autocorrelation  ma¬ 
trix.  Computationally,  the  algorithm  is  a  simple  recursive  procedure  that  avoids 
explicit  matrix  inversion,  decomposition,  or  diagonalization  operations.  When 
the  input  autocorrelation  matrix  is  replaced  by  a  conventional  sample-average 
estimate,  the  algorithm  effectively  generates  a  sequence  of  MF  estimators;  their 
bias  converges  rapidly  to  zero  and  the  covariance  trace  rises  slowly  and  asymp¬ 
totically  to  the  covariance  trace  of  the  familiar  adaptive  matched  filter  (AMF). 
For  finite  data  records,  the  generated  sequence  of  estimators  offers  favorable 
bias/covariance  balance  and  members  of  the  sequence  are  seen  to  outperform 
in  probability  of  detection  (for  any  given  false  alarm  rate)  all  known  and  tested 
adaptive  detectors  (for  example  AMF  and  the  multistage  Wiener  filter  algo¬ 
rithm)  .  These  issues  are  addressed  in  the  context  of  joint  space-time  adaptive 
processing  for  array  radar. 

We  developed  a  new  subspace  direction-of-arrival  (DOA)  estimation  proce¬ 
dure  that  utilizes  a  non-eigenvector  basis.  Computation  of  the  basis  is  carried 
out  by  a  modified  version  of  the  orthogonal  auxiliary- vector  (AV)  filtering  algo¬ 
rithm.  The  procedure  starts  with  the  linear  transformation  of  the  array  response 
scanning  vector  by  the  input  autocorrelation  matrix.  Then,  successive  orthogo¬ 
nal  maximum  cross-correlation  auxiliary  vectors  are  calculated  to  form  a  basis 
for  the  scanner-extended  signal  subspace.  As  a  performance  evaluation  example, 
our  studies  for  uncorrelated  sources  demonstrated  a  gain  in  the  order  of  15dB 
over  MUSIC,  7dB  over  ESPRIT,  and  3dB  over  the  grid-search  maximum  likeli¬ 
hood  DOA  estimator  at  probability  of  resolution  0.9  with  a  ten-element  array 
and  reasonably  small  observation  data  records.  Results  for  correlated  sources 
are  reported  as  well. 

An  alternative  form  of  the  8x8  two-symbol  decodable  quasi-orthogonal  space- 
time  block  code  (QO-STBC)  that  can  be  transmitted  across  either  4  or  8  an¬ 
tennas  with  full  rate  and  full  diversity  order  is  presented.  For  the  8  transmit 
antenna  system,  we  derive  a  new  expression  for  the  rotation  angles  that  max¬ 
imize  the  diversity  (eigenvalue)  product.  In  addition,  we  show  that  the  previ¬ 
ously  proposed  sum-eigenvalue  maximization  criterion  for  the  design  of  rotation 
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angles  is  not  relevant /applicable  and  we  suggest,  as  an  alternative,  minimum 
eigenvalue  maximization.  Finally,  working  directly  with  the  pairwise-error- 
probability  (PEP)  upper-bound  expression,  we  obtain  new  true  PEP-upper- 
bound  optimal  rotation  angles.  For  4  transmit  antenna  systems  and  correlated 
channel  fading  conditions,  we  modify  our  PEP-upper-bound  to  account  for  chan¬ 
nel  correlation.  Using  the  new  PEP-upper-bound  we  obtain  rotation  angles  that 
maximize  the  diversity  product  and  find,  contrary  to  previous  results,  that  the 
optimized  angles  are  independent  of  the  correlation  coefficient.  Simulation  stud¬ 
ies  initiated  herein  demonstrate  the  advantage  of  using  the  proposed  codeword 
across  4  transmit  antennas  when  compared  with  other  4x4  QO-STBC  transmis¬ 
sion  schemes.  For  8  transmit  antennas,  the  studies  compare  the  three  selected 
rotation  angle  optimization  criteria  (diversity  product,  minimum  eigenvalue, 
PEP-upper-bound) . 

Finally,  we  considered  the  problem  of  signature  waveform  design  for  code 
division  medium-access-control  (MAC)  of  wireless  sensor  networks  (WSN).  In 
contrast  to  conventional  randomly  chosen  orthogonal  codes,  we  developed  an 
adaptive  signature  design  strategy  under  the  maximum  pre-detection  SINR  (sig¬ 
nal  to  interference  plus  noise  ratio)  criterion.  The  proposed  algorithm  utilizes 
slowest  descent  cords  of  the  optimization  surface  to  move  toward  the  optimum 
solution  and  exhibits,  upon  eigenvector  decomposition,  linear  computational 
complexity  with  respect  to  signature  length.  Numerical  and  simulation  studies 
demonstrate  the  performance  of  the  proposed  method  and  offer  comparisons 
with  conventional  signature  code  sets. 


11 


II.  Adaptive  Optimization  of  Binary/Quaternary 
CDMA  Signatures  in  Asynchronous  Multi- 
path  Environments 

Part  of  the  work  reported  in  this  section,  namely  adaptive  optimization  of  binary 
signatures,  has  been  presented  at  IEEE  MILCOM  2005,  Atlantic  City,  NJ. 

In  direct-sequence  code-division  multiple-access  (DS-CDMA)  communica¬ 
tion  systems,  multiple  users/signals  with  individual  identifying  signature  wave¬ 
forms  occupy  the  same  channel  in  frequency  and  time.  Alleviation  of  the  result¬ 
ing  multiple-access-interference  (MAI)  problem  relies  on  proper  design  of  the 
user  signatures  and  receivers  under  a  separate  or  joint  consideration. 

In  the  past,  several  optimization  metrics  and  system  configurations  were  con¬ 
sidered.  For  example,  [1] ,  [2]  minimize  iteratively  the  total-squared-correlation 
(TSC)  of  the  signature  set;  convergence  results  are  reported  in  [3].  A  gen¬ 
eralized  version  of  the  problem  with  a  multiple  base  station  configuration  is 
considered  in  [6]  and  [5].  Signatures  that  maximize  user  capacity  for  a  given 
signal-to-interference-plus-noise  (SINR)  level  are  sought  in  [6].  Under  a  mul¬ 
tipath  channel  assumption,  [7]  carries  out  signature  adaptation  to  maximize 
the  SINR  at  the  output  of  the  RAKE-filter  receivers.  The  same  problem  is 
considered  in  [8]  with  signature  optimization  over  a  reduced-rank  vector  space 
and  minimum-mean-square-error  (MMSE)  receivers.  Signature  adaptation  for 
multiuser  MMSE  minimization  is  described  in  [9]. 

In  digital  communication  systems,  it  is  necessary  to  have  signature  sets  that 
are  defined  over  a  finite  alphabet.  Yet,  all  previous  design  efforts  described 
above  deal  with  complex  (or  real)  field  signatures.  Recently,  TSC-optimal  binary 
signature  sets  were  reported  for  nearly  all  signature  lengths  and  set  sizes  [10]- [4]; 
their  sum  capacity  is  identified  in  [5] .  The  sum  capacity  of  several  other  binary 
designs  is  calculated  in  [6]. 

A  binary  signature  that  has  minimum  sum  of  squared  correlations  with  the 
other  multiple-access  signature  codes  is  also  maximum  SINR  optimal  under 
MMSE  filtering  and  synchronous,  equal-power  CDMA  transmissions  over  (ad¬ 
ditive  white  Gaussian  noise)  ideal  Nyquist  channels.  This  is  certainly  not  the 
case  for  asynchronous  multiple  access  protocols  and/or  multipath  channels.  To 
the  best  of  the  knowledge  of  the  authors,  the  only  finite-alphabet  adaptive  sig¬ 
nature  optimization  effort  for  asynchronous/multipath  channels  reported  in  the 
literature  is  [15].  In  this  present  work,  we  attempt  to  build  on  and  improve 
upon  [15]  working  in  the  binary  {±1}  or  quaternary  {±l,±i}  field.  Specifi¬ 
cally,  we  develop  a  new  adaptive  signature  design  procedure  that  maximizes  the 
output  SINR  of  the  MMSE  receiver  filter.  At  nominal  increase  in  computa¬ 
tional  cost,  the  proposed  procedure  offers  significant  performance  improvement 
over  [15].  In  fact,  as  seen  in  the  studies  included  herein,  our  designs  over  the 
restricted  binary/quaternary  field  exhibit  little  performance  loss  in  comparison 
with  the  theoretical  complex- field  optimal  solutions. 
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A.  Signal  model 

We  consider  a  DS-CDMA  system  with  K  users  and  processing  gain  (sig¬ 
nature  length)  L.  The  uplink  (downlink)  transmission  from  (to)  user  k,  k  = 
0, 1, . . . ,  K  —  1,  is  denoted  by 

OO 

uk  =  ^2  bk(m)2~EkSk  ( t  -  mT )  e (1) 

m—1 

where  bk(m)  €  {— 1,  +1}  is  the  mth  data  bit  (binary  phase-shift-keying  data 
modulation),  Ek  is  the  total  transmission  energy,  and  (f>k  is  the  carrier  phase 
with  carrier  frequency  /c;  sk(t)  is  the  normalized  unit-energy  user  signature 
waveform  with  duration  T  given  by 

L—l 

Sk{t)  =  ^  ~  ITc)  (2) 

1—0 

where  s k(l),  l  =  0, 1, . . . ,  L  —  1,  is  the  value  of  the  Ith  chip  of  the  spreading-code 
vector  of  the  kth  user  (sk  is  in  {±1}L  if  binary  or  {±l,±i}L  if  quaternary), 
ip(t)  is  the  chip  waveform,  and  Tc  =  ^  is  the  chip  period. 

The  user  signals  propagate  over  multipath  additive  white  Gaussian  noise 
channels.  The  path  coefficients  are  modeled  as  complex  Gaussian  random  vari¬ 
ables  (Rayleigh  amplitude  and  uniform  phase)  that  are  independent  across  paths 
and  user  signals  (if  uplink  transmissions  are  considered)  and  remain  constant 
during  the  signature  adaptation  period  of  several  symbol  intervals  (quasi-static 
fading).  The  compound  received  signal  due  to  all  users  after  channel  “process¬ 
ing”  and  carrier  demodulation  is  given  by 


K- 1  JV-l 

r(t)  =  ^2  M  m)  ^2  ®k,nSk  ( t  —  mT  -  nTc  —  rk)  +  n(t)  (3) 

k— 0  m  n— 0 

where  N  denotes  the  number  of  resolvable  paths,  a/c,n?  k  =  0,1, . . . ,  K  —  1,  n  = 
0,1,...,  IV  —  1,  is  the  nth  path,  fcth  user,  fading  coefficient,  Tk  £  [0,  T)  is  the 
relative  delay  of  user  k  with  respect  to  user  0  with  tq  =  0,  and  n(t)  is  a  complex 
additive  white  Gaussian  noise  (AWGN)  process. 

Assuming  synchronization  with  the  signal  of  the  user  of  interest  k,  after 
chip-matched  filtering  of  r(t)  and  sampling  at  the  chip  rate  over  a  multipath 
extended  bit  period  we  obtain  the  data  vector  rk(m)  gCl+jV_1  of  the  form 

rfe(m)  =  \/Ekbk(ni)Hksk  +  zk  +  ik  +  nk,  m=  1,2,...,  (4) 

where  zk  £CL+N~1  represents  comprehensively  multiple-access-interference  (MAI) 
to  user  k,  ik£CL+N~1  denotes  channel  induced  inter-symbol-interference  (ISI), 
and  rife £CL+N~1  is  a  zero  mean  complex  Gaussian  noise  vector  with  covariance 
matrix  E  {nfe.n^}  =  <r2I  (£?{•}  denotes  statistical  expectation  and  H  is  the  con¬ 
jugate  transpose  operator).  In  (4),  ER-Sfe  £  CL+JV~1  is  the  channel  processed 
signature  of  user  k  where 
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ak,  0  0 

ak,  1  ak,  0 


Hfc 


(L  +  N-l)xL 


ak,N-  1  ak,N  —  2 

0  “fc.iV-l 


(5) 


ak,  0 
“fc.l 

“fc.iV-l  _ 

For  data  bit  detection  purposes,  we  consider  linear  minimum-mean-square- 
error  (MMSE)  filter  receivers  with  sign-real-part  detectors, 

h  =  sgn  (Re  {wMMSE>krk})  ,  k  =  0, 1, . . . ,  K  -  1,  (6) 


0  0 
0  0 

0  0 


where  w MMSE,k  =  cR/71Hfcsfc  G  CL+N~1,  Rfc  =  E{ rferf  }  and  c  >  0.  The 
output  SINR  of  the  (maximum  SINR)  filter  w MMSE,k  is  a  direct  function  of 
the  binary/quaternary  signature  s k: 


SINR.MMSE,fc(Sfc) 


E< 

( 

l 

wMMSE,k  (v^fc^fcHfcSfc) 

1 

E\ 

wMMSE,fc  (zfc  +  V-  +  Hfc) 

1 

Ek  (sfHfR^H*sfc) 


(7) 

where  Rfc  is  the  autocorrelation  matrix  of  the  disturbance-only  part  of  the  input 
vector  defined  by 


Rfc  =  E  {(zfc  +  ifc  +  rifc)(zfc  +  ifc  +  nk)H }  .  (8) 

Viewing  the  pre-detection  SINR  expression  in  (7)  as  a  function  of  the  signa¬ 
ture  Sfc  in  {±1}L  or  {±1,  ±.?}L  motivates  a  maximization  search  over  the  binary 
or  quaternary  field,  respectively.  In  the  following  section,  we  present  such  an 
adaptive  signature  optimization  scheme. 

B.  The  algorithm 

A  rather  realistic  assumption  for  CDMA  communication  systems  is  that 
the  ISI  influence  is  negligible  in  comparison  with  the  MAI  and  AWGN  dis¬ 
turbance.  In  this  context,  we  can  safely  approximate  Rfc  in  (8)  by  Rfc  ~ 
E  {(zfc  +  rifc)(zfc  +  nfc)^},  which  makes  Rfc  independent  of  sk  and  greatly  sim¬ 
plifies  the  problem  of  maximizing  SINR.MMS.E,fc(sfc)  with  respect  to  sk.  The 
maximization  criterion  takes  the  form 

Sfc  =  argmax  j  sk  Hk  R^HfcSfc  \  (9) 

sfc 6{±1}I/  or  ^  ^ 
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and  Hf  R^ 1Hfc  is  treated  as  an  Lx L  matrix  term  that  does  not  involve  .  Since 
Rfc  is  Hermitian  4positive  definite  by  definition,  HfRjfHfc  is  also  Hermitian 
positive  definite  and  can  be  Cholesky  decomposed  as  first  suggested  in  [15]: 

Af  Ak  =  Hf  R^Hfc  (10) 

where  Ak  is  an  L  x  L  upper  triangular  matrix. 

Upon  Cholesky  decomposition  of  Hf  RjfHfc,  the  quantity  under  maximiza¬ 
tion  in  (9)  becomes 

sf  Af  Afcsfc  =  ||  AfcSfcH2.  (11) 

Therefore,  the  objective  is  the  selection  of  the  signature  (binary  or  quaternary) 
that  maximizes  the  norm  of  the  vector 


Uk  =  AfcSfc  =  a  k,i  afc,2  ■  •  •  a  k,L 


where  a k,j  represents  the  jth  column  of  A;-,  j  =  1, . . . ,  L.  The  norm-square  in 
question  is  given  by 

llufcl|2  =  llafc,Z,||2  +  l!afc,z,-l  l|2  +  ...  +  ||afe,i||“  + 

2ite{af L_1a.k,Lsk(L  —  l)sk(L)}  +  2ite{af L_2&k ,ls%(L  ~  2)sfc(L)}  +  •  •  ■  + 
2ite{afiafciis^(1)sfc(L)}  +  2i?e{af L_2a.k,L-is*k(L  -  2)s k(L  -  1)}+ 

2ite{af  £_3afc,£_lsfc(-t/  —  3)s k(L  —  1)}  +  . . .  + 

2iMaf iafc,L-lsfc(1)sfc(i  -!)}  +  ■••  +  2i?e{af  1at.j2s^(l)sfc(2)} 

(13) 

where  *  denotes  conjugation. 

The  proposed  algorithm  works  in  a  sequential  manner  assigning  values  to 
chips  sk[L),sk(L  —  1), . . . ,  Sfc(l)  one  at  a  time,  in  an  attempt  to  maximize  the 
norm  in  (13).  Since  (13)  is  quadratic,  we  can  start  by  assigning  s k(L)  =  1  with¬ 
out  loss  of  optimality.  Having  set  Sfc(L),  the  value  of  s k(L  —  1)  that  maximizes 
(13)  is  given  by 

Sfc(i-l)  =  argmax  \  Re{a^L_1ak^Lsl(L  -  l)sk(L)}+ 

sfc(^  —  l)^{d=l}Or  {d=l,=hi}  *■ 

Re{ak,L-2ak,L-lsk(L  —  2)Sfc(L  —  1)}  +  .  •  •  + 

Ke(afiafe,i— iSfc(l)s fc(L  -  1)}  }  .  (14) 

If  we  treat  Sfc(l), . . . ,  s/fL  —  2)  as  fixed  unknown  parameters,  the  value  of 
Sfc(L  —  1)  £  {±1}  that  maximizes  (14)  is 

s k(L  -  1)  =  sgn  {Re{akL_1aktL}sk(L)}  ;  (15) 

the  value  of  sk(L  —  1)  £  {±1,  ±i}  that  maximizes  (14)  is 

if  |Re{afi_1afeiisfc(L)}|  > 

/m(afi,_iafc,l,sfc(i)}|  (16) 

otherwise. 


S k(L  -  1)  = 


C  ssn|/?e{afZ/_1afciZ/sfc(L)}|  , 
sgn  |/m{afZ/_1afciZ/sfc(L)}|  i, 
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Sfc(l) 

Sfc(2) 

s  k(L) 


(12) 


We  proceed  similarly  to  calculate  the  values  of  the  remaining  chips  s k{L  — 
j),  j  =  2, 3, . . . ,  L  —  1,  treating  the  lower  indexed  chips  Sfc(l),  Sfc(2), . . . , s k(L — 
j  —  1)  as  unknown  fixed  parameters.  If  s *.(L  —  j)  G  {±1},  j  =  2, 3, . . . ,  L  —  1, 
then 

sk(L  —  j)  =  sgn  ^Re{a^L_jak}L}sk(L)  +  Re{a^L_jaktL_1}sk(L  -  1)  +  . . .  + 

Re{ak,L-ja>=,L-j+l}sk(L  -  j  +  1)  j  ; 

(17) 

lfs k(L-j)  G  {±1,±*},  j  =  2,3,..., L-  1,  then 

sgn^Re  {a£  L_.ak}Ls k(L)  +  a^L_jak,L-lsk(L  -  1)  +  .  . .  + 

a^L-jafc,t-j+1sfc(L- j  +  l)}|  , 

if  Re  ^a^L_jaktLsk(L)  +  . . .  +  a^L_jaktL_j+1s k(L  -  j  +  1)}|  > 

| Im  ^a^L_jak}Ls k(L)  +  . .  .  +  a^L_jak}L_j+1s k(L  -  j  +  1)} 
sgn^Im  ^a^L_jak}Ls k(L)  +  a^L_jaktL_1s k(L  -  1)  +  . . .  + 
ak}L-jak,L-j+isk(L  -  j  +  1)}}  i,  otherwise. 

(18) 

The  above  sequential  scheme  provides  the  first  signature,  sj,1  \  in  what  we  call 
the  signature  candidate  set.  The  complete  candidate  set  is  created  as  follows. 
For  binary  signatures  and  n  =  2,3 , ...  ,L  —  1,  initialize 

s<n)(L-»+l)  =  -s^1}(L-n+l)  and  s  {k\j)  =  s(k\j),  j  =  L  -  n  +  2, . . . ,  L, 

(19) 

and  optimize  recursively  s  —  n),s  —  n  —  1), . . . ,  sj™^(l)  by  (17).  For 
quaternary  sequences  and  n  =  2,3,...,A  —  l,p  =  l,2,3,  initialize 


g(3(n  2)+P+ i)^L  _  n  +  i)  =  -  n  +  1)  and  s^3("  2)+p+1\j)  = 

j  =  L  —  n  +  2, ...  ,L, 

(20) 

and  optimize  recursively  s^n  2'l+p+1'> _  n),s^n  2',+p+1'1  (T  —  n  —  1), . . . , 
Sf(n"2)+P+i)(1)  by  (18). 

The  described  algorithm  gives  a  set  of  L  —  1  candidate  binary  signatures 
or  3 (L  —  2)  +  1  quaternary  signatures.  Among  them,  we  choose  the  one  that 
maximizes  (13). 


C.  Computational  complexity 

The  computational  cost  to  obtain  the  first  signature  of  the  candidate  set, 
s^\  is  0(L3).  At  first  sight,  it  may  seem  that  to  fill  in  the  candidate  set 
with  L  —  1  binary  signatures  or  3  (A  —  2)  +  1  quaternary  signatures  amounts 
to  0(L4).  Yet,  taking  advantage  of  the  inherent  redundancy  in  (19)  or  (20) 
we  can  keep  the  complexity  at  0(L3).  Consider,  for  example,  the  binary  sig¬ 
nature  case  and  calculate  the  above-the-diagonal  part  of  the  L  x  L  matrix 
T  =  Re{Ak  Afc},  T/  j ,  i,j  =  l,...,L,i<j.  The  entry  T)j,  i  <  j,  corresponds 
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to  the  criterion  function  term  i?e{a^!;afcj}sfc(*)sfc(j)  in  (13).  Maximization  of 
(13)  is  equivalent  to  finding  an  assignment  for  Sfc(Z),  l  =  1, . . . ,  L,  that  maximizes 
the  sum  of  the  matrix  elements  ^i<j  ^i,j  (this  problem  falls  under  the  class  of 
quadratic  integer  programming  problems  which  are  known  to  be  NP-harcl  [16]). 
Once  Tij ,  i  <  j,  are  calculated,  we  can  form  the  first  signature  of  the  candidate 
set  via  the  iteration 

41}(L  _  0  =  sgn{TL-i:L-i+1s^\L  —  l  +  1)  +  . . .  +  TL_lt  Ls^\L)}, 
l  =  l,2,...,L-l,  s«(L)  =  l. 

(n) 

Then,  the  other  candidate  set  signatures  sk  ,  n=2,3,...,L  —  1,  are  calculated 
as  follows: 

4n)(j)  =41}G')’  j  =  L-n  +  2,...,L,  s{k\L-n  +  l)  =  -s{k\L  -  n+  1), 
s  in)(i  -  0  =  sgn{TL_UL_l+1^\L  -1  +  1)  +  ...  +  TL_h 
l  =  n, . . . ,  L  —  1. 

(22) 

The  total  computational  cost  of  the  algorithm  as  implemented  above  is  the 
cost  of  the  Cholesky  decomposition,  plus  calculation  of  Thj ,  i  <  j,  plus  the  cost 
of  iterations  (21)  and  (22).  That  is, 

Computational  complexity  =  ^-  +  L(L  +  ~  +  L<'L~  l)  +  L(-L  ~  l)^2L  ~  H  =  0{L3). 

(23) 

For  quaternary  signatures  the  cost  of  (22)  increases  by  a  factor  of  3,  but  the 
overall  complexity  remains  0(L3).  In  either  case,  binary  or  quaternary,  the 
trade-off  in  keeping  the  computational  complexity  at  0(L 3)  (instead  of  0(L4)) 
is  the  space  requirement  of  0(L2)  for  the  matrix  T.  We  conclude  that  the 
proposed  signature  optimization  procedure  has  little  overhead  in  comparison 
with  [15].  In  fact,  at  a  higher  level  of  system  abstraction  both  schemes  are 
bounded  only  by  the  cost  of  inverting  R*,,  O  ((L  +  N  —  l)3). 

D.  Simulation  studies 

In  this  section,  we  compare  the  performance  of  the  proposed  signature  opti¬ 
mization  algorithm  (binary  and  quaternary)  against  the  following  benchmarks: 
(i)  The  maximum-eigenvalue  eigenvector  of  H^R^H/c,  denoted  by  max-EV, 
which  is  the  theoretical  maximum  SINR  signature  solution  over  the  CL  held  (cf. 
(7));  (ii)  the  binary  and  (iii)  quaternary  minimum  Euclidean  distance  quantized 
versions  of  max-EV;  and  (iv)  the  algorithm  in  [15]  where  binary  conditional 
maximization  of  a  modified  version  of  (11)  was  considered1. 

In  Fig.  1,  we  consider  the  uplink  multipath  DS-CDMA  signal  model  of  Sec¬ 
tion  II  with  spreading  gain  L  =  16  and  we  assume  the  presence  of  K  =  8  users. 
Each  user  signal  experiences  N  =  3  paths  with  coefficents  (zero-mean  complex 
Gaussian  random  variables)  of  equal  power.  Then,  following  the  notation  of 

1In  [15],  ad  hoc  posterior  bit  flipping  tests  after  conditional  maximization  were  executed 
that  improved  the  final  reported  SINR  performance.  This  is  not  considered/accounted  for 
herein. 
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Section  II,  the  total  average  received  SNR  for  user  k,  k  =  0,1, ...  ,7,  over  all 


paths  is 


OATT)  A^fcELo^{Kn|2}  EkE{  ||afc||2} 

oiNKfc  =  - - -  =  - - - . 


We  set  SNR0-2  =  8dB,  SNR3_5  =  9dB,  SNR6j7  =  lOdB.  We  initialize  the  signa¬ 
ture  set  so,  si, . . . ,  S7  arbitrarily  and  execute  each  signature  set  design  algorithm 
sequentially  user-after-user  in  what  we  call  a  multiuser  adaptation  cycle.  Several 
multiuser  adaptation  cycles  are  carried  out.  A  fixed  Walsh-Hadamard  signature 
assignment  is  also  included  in  the  study  to  challenge,  potentially,  the  notion  of 
signature  adaptivity.  In  Fig.  1,  we  plot  the  pre-detection  SINR  of  a  representa¬ 
tive  user  per  adaptation  cycle  averaged  over  1,000  random  channel  realizations 
and  initial  signature  assignments.  All  algorithms  are  seen  to  converge  in  about 
three  cycles.  The  proposed  procedure  offers  superior  performance  for  both  bi¬ 
nary  and  quaternary  alphabets.  The  (non-adaptive)  Walsh-Hadamard  signature 
assignment,  as  expected  arguably,  exhibits  rather  poor  performance.  It  is  inter¬ 
esting  to  observe  that  under  the  proposed  optimization  scheme,  upgrading  the 
signature  domain  from  binary  to  quaternary  provides  close  to  0.8dB  gain,  while 
upgrading  the  domain  from  quaternary  to  the  full  complex  field  provides  only 
0.52dB  additional  gain. 

In  Fig.  2,  we  repeat  the  same  study  for  a  system  with  processing  gain 
L  =  31  and  K  =  16  users  with  SNR  values  SNR0-3  =  8dB,  SNR.4_7  =  9dB, 
SNRg-n  =  lOdB,  SNR12-15  =  lldB.  A  (non-adaptive)  Gold  signature  assign¬ 
ment  is  included  in  the  comparisons  which,  arguably  contrary  to  common  belief, 
fails  disappointingly  to  withstand  the  asynchronous  multipath  environment.  In 
terms  of  the  adaptive  signature  design  schemes,  the  conclusions  are  similar  to 
the  ones  drawn  from  Fig.  1. 

In  Fig.  3,  we  repeat  the  study  of  Fig.  2  under  finite  sample  support  esti¬ 
mation  of  the  signal  autocorrelation  matrix  and  plot  the  average  pre-detection 
SINR  against  data  record  size  after  the  third  multiuser  adaptation  cycle.  The 
dominance  of  the  proposed  signature  design  scheme  (that  does  not  rely  on  eigen¬ 
vector  estimates)  is  now  even  more  prominent.  In  fact,  both  Cholesky-based 
algorithms  (the  proposed  binary/quaternary  and  the  binary  conditional  maxi¬ 
mization  procedure  of  [15])  experience  lower  performance  degradation  than  the 
eigen-decomposition-based  schemes. 

Finally,  to  illustrate  the  performance  of  the  signature  design  algorithms  un¬ 
der  single-user  adaptation,  in  Fig.  4  we  plot  the  bit-error-rate  (BER)  of  a  repre¬ 
sentative  user  versus  SNR  when  the  signal  autocorrelation  matrix  is  estimated 
from  3  (L  +  N  —  1)  samples  and  all  interferers  have  fixed,  non-adapted  Gold 
signature  assignments.  The  comparative  findings  are  no  different  than  those 
obtained  in  Figs.  1-3  under  the  SINR  metric  and  multiuser  adaptation. 


E.  Conclusions 

We  considered  the  problem  of  adaptive  signature  optimization  for  CDMA 
communications  in  asynchronous  mode  of  operation  and/or  multipath  signal 
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propagation  environments.  Contrary  to  past  and  current  research  literature,  we 
pursued  finite-alphabet  signature  designs,  binary  and  quaternary,  in  particular. 

Finite-alphabet  optimization  of  signatures  is  NP-hard.  The  proposed  subop- 
timal  algorithm  described  and  studied  in  this  report  has  computational  complex¬ 
ity  0(L3)  and  storage  complexity  0(L2)  where  L  is  the  signature  length  under 
consideration.  The  algorithm  can  be  viewed  as  a  truly  practical  means  for  either 
single-user  signature  optimization  or  multi-user  signature  set  optimization  since 
no  exchange  of  information  between  users  is  required. 

Performance-wise,  it  is  satisfying  -and  somewhat  surprising  arguably-  to  see 
that  the  quaternary  adaptive  designs  are  less  than  ldB  away  in  attained  SINR 
from  the  theoretical  complex-field-optimal  signatures.  Being  a  Cholesky-based 
procedure,  the  algorithm  has  favorable  small-sample-support  behavior;  3 L  or  4L 
data  vectors  suffice  for  effective  adaptation,  which  is  well  within  the  coherence 
time  of  common  commercial  wireless  channels. 
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Pre-detection  SINR  (dB) 


Multiuser  adaptation  cycle 


Fig.  1:  Average  pre-detection  SINR  for  a  representative  user  versus  multiuser 
adaptation  cycle  (L  =  16,  I\  =  8). 
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Pre-detection  SINR  (dB) 


Multiuser  adaptation  cycle 


Fig.  2:  Average  pre-detection  SINR  for  a  representative  user  versus  multiuser 
adaptation  cycle  (L  =  31,  I\  =  16). 
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Pre-detection  SINR  (dB) 


Data  record  size  [x  (L+N-1)] 


Fig.  3:  Average  pre-detection  SINR  for  a  representative  user  versus  sample 
support  ( L  =  31,  K  =  16,  three  multiuser  adaptation  cycles). 
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BER 


Fig.  4:  BER  versus  SNR  under  single-user  adaptation  ( L  =  31,  K  =  16,  fixed 
Gold  MAI  signatures,  3 (L  +  N  —  1)  sample  support). 
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III.  New  Bounds  on  the  Periodic  Total  Squared 
Correlation  of  Binary  Signature  Sets  and 
Optimal  Designs 

This  work  has  been  presented  at  IEEE  MILCOM  2005,  Atlantic  City,  NJ. 

In  code-division  multiplexing  applications,  for  example  direct-sequence  code¬ 
division  multiple-access  (DS-CDMA)  cellular  communications  systems,  each  of 
the  K  participating  signals/users  is  equipped  with  a  unique  identifying  signature 
vector  sfc  <E  CL,  ||sfc||  =  1,  k  =  1,2, . . . ,  K.  All  signatures  put  together  in  the 
form  of  a  matrix  define  what  we  call  the  signature  set 

S  =  [Sl  S2...SA-]  eCLxK.  (25) 

In  synchronous  code-division  multiplexing  transmissions  over  well  behaved 
Nyquist  channels,  we  are  interested  in  using  a  signature  set  with  the  smallest 
possible  total  squared  correlation  (TSC)  value  [2]- [7] 

TSC(S)  =  EE  k"s/  (20 

i=  1  i= 1 

where  H  denotes  the  Hermitian  operator.  For  complex/real-valued  signature 
sets  5  e  CLxK  or  RLxK,  if  K  >  L ,  TSC (5)  >  ^  [1]  (of  course,  TSC (5)  >  K 
if  K  <  L).  Overloaded  ( K  >  L)  sets  with  TSC  equal  to  have  been  known 
as  Welch-bound-equality  (WBE)  sets.  Algorithms  and  studies  for  the  design  of 
complex  or  real  WBE  signature  sets  can  be  found  in  [3]- [7].  In  digital  trans¬ 
mission  systems,  however,  it  is  necessary  to  have  finite-alphabet  signature  sets. 
Recently,  new  bounds  were  derived  on  the  TSC  of  binary  antipodal  signature 
sets  together  with  optimal  designs  for  arbitrary  signature  lengths  and  set  sizes 
[2]- [4].  The  sum  capacity,  total  asymptotic  efficiency,  and  maximum  squared 
correlation  of  minimum-TSC  optimal  binary  sets  were  evaluated  in  [5].  The 
sum  capacity  of  several  other  signature  set  designs  under  potentially  a  binary 
or  quaternary  alphabet  was  examined  in  [6].  All  developments  that  follow  in 
this  present  report  refer  to  binary  signature  sets. 

When  asynchronous  code-division  multiplexing  is  attempted  and/or  the  chan¬ 
nel  exhibits  multipath  behavior,  apart  from  the  total  squared  correlation  be¬ 
tween  signatures  we  are  also  concerned  about  the  individual  periodic  (cyclic) 
auto-correlation  and  the  periodic  (cyclic)  cross-correlation  values  [8] .  For  nota- 
tional  simplicity,  let  sk\i  denote  the  cyclic  right-shifted  version  of  s*,  €  {±1}L  ,  k  = 
1,2,...,  A",  by  l  bit  positions,  l  =  0,1,2,...  (hence,  sfc|0  =  sk\L  =  ...  =  sk). 
First,  we  define  the  cyclic  extension  matrix  Sc  €  {±l}LxAL  of  the  signature  set 
5e  {±i}LxK 
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*5c  —  [si|0  S2|0  •  ■  •  S/f|o  S;l | !  s2|l  •  ■  •  sx n 


S1|L-1  s2|i-l  •  ■  •  sK|i— l]  •  (27) 


Then,  we  define  the  periodic  total  squared  correlation  (PTSC)  of  the  signature 
set  S  as  the  total  squared  correlation  (TSC)  of  the  matrix  Sc 


PTSC(S)  =  TSC(Sc). 


(28) 


Since, 

si\hsi\l2  =  si\osj\h-h  =sj\osi\L-(l2-h),  h  —  ^2>  h,h  =  0, 1,  2, . .  . ,  L—  1,  i,j  =  1,2,. 
where  T  denotes  the  transpose  operator,  we  can  calculate 


. ,  K, 
(29) 


PTSC(S)  = 


A'L3  +  2L£f=i 

EA  lsfcsfc|i|2  + 

2LT,r=iT,?=i,i<j'EiJo  \sTsi\i 

2  ,  L  =  1  (mod  2) 

A-L3  +  L£Ei 

2EA1  lsfcsfc|i|2+  sfesfc|A 

+ 

2iEAiEf=i,i<JEf=01I^K 

2  ,  L  =  0  (mod  2) 

The  first  two  terms  of  the  PTSC  expressions  (for  odd  or  even  L)  contain  all 
periodic  auto-correlation  contributions.  The  third,  triple  summation  term  con¬ 
tains  all  periodic  cross-correlation  contributions.  Minimizing  PTSC  addresses 
the  problem  of  minimizing  periodic  auto-  and  cross-correlations  jointly.  In  the 
present  work,  we  derive  new  lower  bounds  on  the  PTSC  of  binary  signature  sets 
for  all  possible  values  of  K  (number  of  signatures)  and  L  (signature  length). 
Then,  we  derive  optimal  designs  for  several  ( K ,  L)  pair  cases.  The  designs  are 
based  on  Hadamarcl  matrix  transformations  (an  approach  followed  in  [2]  and 
[14],  for  example)  and  serve  as  proof-by-construction  for  the  tightness  of  the 
pertinent  PTSC  bounds. 


A.  New  Bounds  on  the  PTSC  of  Binary  Antipodal  Signa¬ 
ture  Sets 

Consider  the  cyclic  extension  matrix  Sc  £  {±l}ixAL  in  (137)  and  denote  its 
zth  row  by 

=  [si|o(*)  •  ■  •  snr|o(*)  si|i (*)  •  •  -s/qi(*)  .  spi-ib)  •  •  -s/fiL-d*)] T  =  b  2,  •  ■  ■ ,  L, 

(31) 

where  k  =  1, 2, . . . ,  K,  l  =  0,  2, . . . ,  L  —  1,  refers  to  the  zth  element  of 

vector  Sfe| Then,  by  (28)  and  the  “row-column  equivalence”  for  the  TSC  metric 
of  matrices  [7], 


PTSC(S)  =  ^3  +  E^E^|dfd,|2.  (32) 

Our  goal  is  to  obtain  a  bound  on  the  term  E-[=i  EjA  j^i  |  df  dj  | 2  for  all  (K,  L) 
cases.  For  this  purpose,  we  will  make  use  of  the  following  theorem. 
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Theorem  1  ( On  the  Properties  of  Cyclic  Extension  Matrices )  Let  Sc  = 

[di  d2  . . .  d l]T,  di  G  {±l}Ai  ,  i  =  1, 2, . . . ,  L,  be  the  cyclic  extension  matrix 
of  S  =  [s-i  s2  . . .  s K] ,  sfc  G  {±1}L,  k  =  1,2,...,  K. 

(i)  Then, 

K 

di  d j  ^  ^  s fc  s^i  j_i ,  i)  j  1,2,...,L,  (33) 

k= 1 

and 


E  E  d^ 


2LE,: 


-+1 


dfdjl 


2EJL2|dfdjr  + 


L  =  1  (mod  2) 
L  =  0  (mod  2). 


(34) 


(ii)  If  KL  ^  0  (mod  4),  |dfdj|  ^  0  for  all  j  =  2,3 Specifically,  if 
KL  =  1  (mod  2),  |df  dj  |  =  1  (mod  2)  for  all  j  =  2, 3, . . . ,  L;  if  KL  =  2  (mod  4), 
|df  dj  |  =2  (mod  4)  for  all  j  =  2,  3, . . . ,  L.  □ 


From  (32)  and  Theorem  1,  Part  (i),  the  PTSC  of  a  binary  signature  set  S  equals 


PTSC(S)  = 


K2L3  +  2LJ2- 


!-+! 


d(  d, 


K2L 3  +  L 


2  Ej^=2  |di  dj 


L  =  1  (mod  2) 

(35) 

L  =  0  (mod  2) 


and  — interestingly—  depends  only  on  the  sum  of  the  periodic  auto-correlations 
of  the  signatures  for  each  shift  (dfdj  =  Et=i  sk] sfc|y— l ,  j  =  2,3  By 

(35)  and  Theorem  1,  Part  (ii), 

(  K2 L3  +  2L  ,  L=  1  (mod  2),  K  =  1  (mod  2) 

PTSC(S)  >  <  A'2L3  +  8L(LzL)  ,  L  =  1  (mod  2),  A' =  2  (mod  4) 

{  K2L3  +  8L(L  -  1)  +4A,  Lee  2  (mod  4),  A'  =  1  (mod  2). 

(36) 

Simplifying/compressing  (36)  and  adding  the  entry  KL  =  0  (mod  4),  we  con¬ 
clude 

(  I<2L 3,  A'L  =  0  (mod  4) 

PTSC(S)  >  <  K2L3  +  4L(L  —  1),  KL  =  2  (mod  4)  (37) 

[  K2L3  +  L(L~-  1),  A'L  =  1  (mod  2). 

Expression  (37)  defines  our  new  bounds  on  the  periodic  total  squared  correlation 

of  binary  antipodal  signature  sets.  In  the  following  section,  we  develop  signature 
set  designs  that  meet  these  bounds  with  equality. 


B.  Design  of  Minimum  PTSC  Binary  Antipodal  Signature 
Sets 

In  view  of  the  formulation  introduced  in  Section  I,  minimum  PTSC  design  of  a 
signature  set  S  is  equivalent  to  minimum  TSC  design  of  an  extension  matrix  Sc. 
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In  the  following,  H^r  refers  to  a  normalized  Hadamard  matrix2  of  size  N  with 
columns  hj,  i  =  1,2 Vec[A]  denotes  the  standard  column- by-column 
vectorization  operation  on  matrix  A  and  diagjb}  is  the  square  diagonal  matrix 
with  the  elements  of  vector  b  on  its  diagonal.  Finally,  we  find  it  useful  to  in¬ 
troduce  the  column  replication  operator  A p,  p  =  1,2,....  If  A  =  [ai  a2  . . .  aj], 

then  AAp  =  [ai  ai . . .  ai  a2  a2  . . .  a2  . . .  a,/  a j  . . .  ajj . 

p  p  p 

Below,  we  provide  PTSC-optimal  signature  set  designs  under  several  K,  L 
cases  for  both  underloaded  ( K  <  L)  and  overloaded  (K  >  L)  systems. 


Case  1:  Underloaded  Systems  (K  <  L) 

Case  1:  K  =  0  (mod  2),  L  =  2",  n=  1,2,... 

Calculate  N  =  4  \Jj\  {K  is  rounded  down  to  the  nearest  multiple  of  four). 
Decompose  N  into  N  =  Ni  +  _/V2  +  . . .  +  Nj  where  Ni  >  A2  >  . . .  >  Nj  are 
powers  of  2.  Obtain  a  Hadamard  matrix  H  v,  and  create  the  initial  template 
matrix 

&LxNl  =  [nNlnNl...HNl)T.  (38) 

v - v - ' 

L 

77T 

Define  the  diagonal  “correction”  matrices 


Ci 


diag  <  Vec 


Vec  [H 


JViJ 


T/\JV  n 
1 


and  the  final  diagonal  correction  matrix 


1,2,...,  | logNlL-  1J  , 
(39) 


X  =  diag 


Vec 


r 

I 

Vec 

hi  h2  ...  h  l 

i 

_ 

. 

V 

(40) 


Calculate 


S'  =  ©' 


^  Lj°9JVi  L— lj 

n  c* 

•-1 


x. 


(41) 


Repeat  (38)-(41)  for  {Nj}'.=2  to  create  {Sj}j=2-  ^:i i  3  =  1,  2,  ■  •  ■  J,  is  an 

L  x  Nj  PTSC-optimal  set.  Concatenate  to  form 


Slxn  =  [<5i  <S2  . . .  Sj } .  (42) 

2  We  recall  that  a  Hadamard  matrix  Hjv  is  an  N  x  N  matrix  with  elements  +1  or  —1  and 
orthogonal  columns,  H^H^v  =  Nljy  where  I  at  is  the  size- AT  identity  matrix.  A  necessary 
condition  for  a  Hadamard  matrix  to  exist  is  that  its  size  is  a  multiple  of  four,  except  for  the 
trivial  cases  of  size  one  or  two. 
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If  K  =  0  (mod  4),  then  S  in  (42)  is  an  L  x  K  PTSC-optimal  signature  set.  If 
K  =  2  (mod  4),  design  PTSC-optimal  sets  S LxrK_2)  and  $Lx2  as  above  and 


concatenate3; 


5'  5" 


is  an  L  x  K  PTSC-optimal  signature  set. 


□ 


Case  2a:  K  =  0  (mod  8)  and  I\  >  8 


viz 

8 


L  ^  2  (mod  4)  and  L  ^  2n,  n  = 


1,2,... 


2b:  K  =  0  (mod  8)  and  8 

~  vzz~ 
8 

<K  <  4  [-§]  ,  L  =  2  (mod  4) 

Calculate  N  =  ^,  q  =  ^ 

,  and  r  =  L  —  N  .  Obtain  a  Hadamard 

matrix  Hjy  and  create  the  initial  template  matrix 


&lxn  =  [HjyHjv__.  Hjy  hx  h2  . . .  hr}T.  (43) 

<? 


If  r  ^  0,  define  the  diagonal  correction  matrices 


C i  =  diag 

|  Vec 

1 

cr 

tr 

to 

ir 

(-l)i+1  (-l)i+1 . . 

■(-i  )i+1 

T~ 

> 

i 

1 

r 

- 

> 

i  =  1,2. 


(44) 


If  r  =  0,  define  the  diagonal  correction  matrices  as  follows: 

C4  =  diag  { Vec  [  [(-l)i+1hx]  h2  . . .  hq]  }  ,  i  =  1, 2.  (45) 

Calculate  Sj  =  @TCj,  i  =  1,2,  and  concatenate  to  form  the  PTSC-optimal 
set 


Slxk  =  [5i  52] . 


(46) 


□ 


Vl 

4 


Case  3a:  K  =  0  (mod  4)  and  4 

3b:  K  =  4  (mod  8)  and  K  >  8 


<K< 


vTl 


,  L  =  0  (mod  K) 


8 


,  L  =  0  (mod  A') 


Calculate  q  =  jr.  Obtain  a  Hadamard  matrix  H k  and  create  the  initial 
template  matrix 

@LxK  =  [H K  Ha-  ■  ■  ■  (47) 

<? 

Define  the  diagonal  correction  matrix 


C  =  diag  {Vec  [hi  h2  . . .  h9]} 


(48) 


and  calculate  the  PTSC-optimal  set 

ST  =  ©TC.  (49) 

3 Notice  that  repetitive  concatenation  of  the  SLx2  PTSC-optimal  set  would  yield,  math¬ 
ematically,  PTSC-optimal  sets  for  any  L  x  K-even  dimensions;  it  would  do  so,  however,  via 
signature  repetition  which  is  unacceptable  application-wise. 
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□ 


Case  4:  K  =  0  (mod  4)  and  K  >  4  ,  L  =  2  (mod  4) 


Calculate  r  =  L  —  K .  Obtain  a  Hadamard  matrix  and  create  the  initial 
template  matrix 

®t,xk  =  [Hjf  hi  h2  . . .  hr]T. 


"LxK 

Define  the  diagonal  correction  matrix 

T 


C  =  diag  ^Vec 


Calculate  the  PTSC-optimal  set 


1-11  -1...1  -1 


S1  =  &1  C. 


(50) 


(51) 


(52) 

□ 


Some  of  our  designs  that  follow  utilize  known  cyclic  Hadamard  sequences 


g m  G  {±1}M  of  length  M  =  3  (mod  4)  that  have  ideal  two- level  periodic  auto¬ 
correlation, 


gAf|m 


M,  m  =  0 

— 1,  m  =  1,  2,  1, 


(53) 


and  have  been  used  in  the  past  to  construct  cyclic  Hadamard  matrices  [16]- [18]. 
In  particular,  if 

A4  =  {M  :  (i)  M  prime  congruent  to  3  (mod  4)  or 

(ii)  M  congruent  to  3  (mod  4)  and  product  of  “twin  primes”  p  and  p  +  2  or 
(in)  M  =  2n  —  1,  n  =  2,  3, . . .  }  , 

(54) 

then  there  exists  at  least  one  systematic  method  [16]- [18]  to  construct  such  a 
sequence  g m  when  M  £  M.  We  continue  with  the  presentation  of  the  design 
cases. 


Case  5:  K  =  2  (mod  8)  and  K  >  8 


viz 

8 


+  2,  L&M 


Obtain  the  signature  set  SLx(K_2)  from  (43)-(46)  and  form  the  PTSC- 
optimal  set 


Slxk  — 


S  g L  gL|l 


(55) 

□ 


Case  6:  K  =  1  (mod  8),  K  >  8 


+  1  and  K^L-2,  LeM 


Obtain  the  signature  set  SLx(K_i\  from  (43)-(46)  and  form  the  PTSC- 
optimal  set 


Slxk  — 


S  g  L 


(56) 
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□ 


Case  7:  K  =  L,  LeM 

Set  TV  =  K  +  1,  obtain  a  Hadamard  matrix  H^r  and  remove  its  first  row 
and  column.  Call  the  resulting  size- K  matrix  H^.  “Correct”  with  the  diagonal 
matrix 

C  =  diag  {gi}  (57) 

to  obtain  the  PTSC-optimal  set 

ST  =  HnC.  (58) 

□ 


Case  8:  K  €  {L  -  4,  L  -  2}  ,  L  €  M 
Design  directly  the  PTSC-optimal  set 


- 

T 

11  ...1  - 1 

gL  gL|l  •  •  •  &L\K-2 

L—l 

(59) 

□ 


Next,  we  proceed  with  the  presentation  of  our  optimal  designs  for  overloaded 
systems. 

Case  2:  Overloaded  Systems  (I<  >  L) 

Case  1:  K  =  0  (mod  4) 

Obtain  a  Hadamard  matrix  H^-  and  trivially  design  the  PTSC-optimal  set 

Slxk  =  [hi  h2  . . .  h l]t  .  (60) 

Hence,  interestingly,  overloaded  direct  Hadamard  designs  are  optimal  under 
both  the  PTSC  and  TSC  [2]  metrics.  □ 


Case  2:  K  =  2  (mod  4),  L  £  M 

Set  TV  =  K  —  2,  obtain  a  Hadamard  matrix  Hjv  and  keep  only  the  first  L 
columns.  Call  the  resulting  matrix  Hy  Then, 


Slxk 


H\  g L  gi|l 


(61) 


is  PTSC-optimal. 


□ 


Case  3:  K  =  1  (mod  4),  L  £  M 

Set  TV  =  K  —  1,  obtain  a  Hadamard  matrix  H  ,v  and  keep  only  the  first  L 
columns.  Call  the  resulting  matrix  Hy  Then, 


Slxk 


Hjv  g  L 


(62) 
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is  PTSC-optimal. 


□ 


Case  4:  K  =  3  (mod  4),  L  G  M 

Set  N  =  K  +  1,  obtain  a  Hadamard  matrix  H^v  and  remove  the  first  row 
and  N  —  L  columns.  Call  the  resulting  matrix  H^.  “Correct”  with  the  diagonal 
matrix 

C  =  diag  {gi}  (63) 

to  obtain  the  PTSC-optimal  set 

ST  =  HnC.  (64) 

□ 

In  the  following  section,  we  discuss  these  design  findings  and  present  some 
examples. 

C.  Discussion  and  Examples 

The  optimal  design  cases  presented  in  the  previous  section  constitute  proof-by¬ 
construction  of  the  tightness  of  the  corresponding  PTSC  bounds  developed  in 
Section  II.  To  acquire  a  quantitative  feeling  of  the  coverage  of  the  presented 
designs,  if  we  restrict  the  domain  of  K,  L  to  {1, 2, ... ,  256}  (at  present,  it  does 
not  appear  of  practical  interest  to  consider  code-division  applications  outside 
this  parameter  range),  we  can  calculate  that  Underloaded  Cases  1  through  8 
and  Overloaded  Cases  1  through  4  together  represent  26.61%  of  all  possible 
combination  pairs  (K,  L )  €  {1,2,...,  256}2.  Certainly,  tightness  of  the  bounds 
and  optimal  PTSC  designs  under  the  remaining  cases  is  an  important  open 
research  problem. 

Direct  comparison  of  our  PTSC  optimal  designs  with  the  TSC  bounds  and 
optimal  sets  in  [2]  shows  that  Underloaded  Case  1  when  K  is  a  power  of  2, 
Underloaded  Case  3  and  Overloaded  Cases  1,  3,  and  4  are  doubly,  both  PTSC 
and  TSC,  optimal.  Furthermore,  we  can  now  establish  that  the  familiar  Gold 
sets  [15],  which  have  been  widely  used  for  their  periodic  correlation  properties 
[8],  are  not  PTSC  optimal  (ironically,  Gold  sets  were  shown  to  be  TSC  opti¬ 
mal  in  [2]).  To  that  respect,  we  recall  [8],  [15]  that  Gold  sets  are  defined  for 
I<  <L  +  2,  L  =  2"  —  1,  n  >  3  and  n  ^  0  (mod  4),  and  for  every  i  ^  j,  i,j  = 
1  the  periodic  signature  cross-correlations  |sfsj|j|  ,  l  =  Q,...,L  —  1, 

have  value  —  1 ,  —  2  L “2“  J  —  1 ,  or  2  L  J  —  1  with  frequency  of  occurrence  ( “cor¬ 
relation  spectrum”  [8])  2"  —  2n~a  —  1,  2n~a~ 1  —  2  s  and  2n~a~ 1  +  2  3  , 

respectively,  where  a  =  1  if  n  =  1  (mod  2)  and  a  =  2  if  n  =  2  (mod  4).  Con¬ 
sider  such  a  Gold  set  Glxk  =  [si  s2  . .  .s^-]  where  si,  s2  are  the  two  preferred 
m-sequences  [21]  and  call  its  cyclic  extension  matrix  Qc.  By  (28)  and  (30), 

K  K  K  L—l 

PTSC{G)  =  TSC(Gc)  =  KL3+2L  H  X!  Isfcs*rl  +2LH 

k—1  1=1  i= 1  j= 1,  %<j  1=0 

(65) 
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and  utilizing  the  cross-correlation  spectrum  information 


-  2  2" 


Ef^l^nl2  =  [(2"-2™-«-l)(-l)2+(2™- 

^—2 J  -  l)2  +  (2n~a~1  +  2j^2'J  -  iy 


=  (2n  —  l)2  +  2”  —  2. 

Regarding  the  auto-correlation  contributions,  if  k  =  1,  2, 


(66) 


2^  sjfs*|,  =  2”  —  2; 


(67) 


if  k  >  2, 


2  |^fc|;|2  >  (-1)2  +  (-2ffl  -  I)'  +  (2ffl  -  l)2  >  2™+2  +  3.  (68) 


i=i 


For  K  <  2,  combining  (65),  (66),  and  (67)  we  calculate 

PTSC(Q)  =  KL 3  +  KL  (2n  -  2)  +  LK(K  -  1)  [(2"  -  l)2  +  2n  -  2] 
=  K2L3 +KL(L-1)+LK(K 


(69) 


For  K  >  2,  combining  (65),  (66),  and  (68)  we  obtain 


PTSC{G )  >  KL3  +  2L(2"  -  2)  +  (AT  -  2)L  (2"+2  +  3) 

+LK(I<  -  1)  [(2n  -  l)2  +  2n  -  2] 

=  K2L3  +  2 L(L  -  1)  +  (A'  -  2)L(4A  +  7)  +  LK{K  -  1  )(L  -  1) 
>  K2L3  +  K2L(L  -  1). 


(70) 

Expression  (69)  for  I\  =  2  and  expression  (70)  for  I\  >  2  are  strictly  greater 
than  the  lower  bounds  of  Section  II  which  are  achievable  for  some  Gold-compatible 
values  (AT,  L)  as  seen  in  Section  III.  Fig.  5  shows  as  an  example  a  (16,  31)  Gold 
set  which  has  PTSC  =  8213760  together  with  our  optimal  (16,  31)  design  (under 
Underloaded  Case  2a)  with  minimum  PTSC  =  K2L3  =  7626496. 

Other  signature  sets  well  known  for  their  periodic  correlation  properties 
are  the  small  and  large-set  Kasami  designs  [8],  [16].  We  recall  that  small-set 
Kasami  designs  are  defined  for  lengths  L  =  2™  —  1,  n  =  0  (mod  2),  and  have  size 
I\  <  2?.  The  attained  periodic  cross-correlation  values  are —1,  —2^—1,  2^— 1, 
but  their  frequency  of  occurrence  is  not  known  in  closed  form  as  a  function  of 
n.  Fig.  6  presents  a  (2, 15)  small-set  Kasami  design  with  calculated  PTSC  = 
15300  together  with  our  optimal  (2, 15)  design  (under  Underloaded  Case  5)  with 
minimum  PTSC  =  K2L 3  +  4 L(L  —  1)  =  14340.  We  directly  conclude  that,  in 
general,  small-set  Kasami  designs  are  not  PTSC  optimal. 

Large-set  Kasami  designs  [8],  [16]  are  defined  for  L  =  2"  —  1  and  even  n.  If 
n  =  2  (mod  4),  K  <  2^(2”  +  1);  if  n  =  0  (mod  4),  K  <  2%(2n  +  1)  -  1.  Fig. 
7(a)  shows  an  (8,15)  large-set  Kasami  design  that  has  PTSC  =  227040.  Fig. 
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7(b)  shows  our  PTSC-optimal  set  S°$x  8  designed  under  Underloaded  Case  2a 
with  minimum  PTSC  value  216000.  Hence,  in  general,  large-set  Kasami  are  not 
optimal  either. 

We  conclude  this  section  with  an  example  of  an  overloaded  PTSC-optimal 
design  Sgf*  42  given  in  Fig.  8.  The  set  is  designed  under  our  Overloaded  Case 
2  procedure  and  has  minimum  PTSC  value  52555044. 

D.  Conclusions 

We  derived  bounds  on  the  periodic  (cyclic)  total  squared  correlation  (PTSC)  of 
binary  signature  sets  for  any  signature  length  L  and  set  size  K  and  provided 
optimal  constructions  for  a  variety  of  K,  L  values  that  establish  the  tightness 
of  the  corresponding  bounds.  The  constructions  include  underloaded  ( K  <  L) 
and  overloaded  (K  >  L)  design  cases  and  cover,  as  an  example,  26.61%  of  all 
possible  combinations  of  K,  L  in  {1,2,...,  256}. 

Side  results  of  the  presented  research  include  derivation  of  a  lower  bound 
on  the  PTSC  of  Gold  sets  which  is  seen  to  be  strictly  larger  than  the  gen¬ 
eral  derived  bounds.  PTSC-optimal  constructions  described  herein  for  Gold- 
compatible  ( K ,  L)  pairs,  as  well  as  small  and  large-set  Kasami,  establish  that 
neither  the  Gold  nor  the  Kasami  sets  are  PTSC-optimal  in  general.  In  view 
of  these  findings,  the  developed  PTSC-optimal  sets  take  precedence  whenever 
the  periodic  correlation  properties  of  signatures  is  of  concern  in  code-division 
multiplexing  applications. 
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Fig.  5.  (a)  03ixi6  Gold  set  with  PTSC  =  8213760.  (b)  Optimal  signature 
set  designed  under  Underloaded  Case  2a  with  PTSC  =  (16) 2 (31)3  = 

7626496. 
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Fig.  6.  (a)  /Cffx2  small-set  Kasami  with  PTSC  =  15300.  (b)  Optimal  signa¬ 
ture  set  <S°5x 2  designed  under  Underloaded  Case  5  with  PTSC  =  (2)2(15)3  + 
4(15)(14)  =  14340. 
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Fig.  7.  (a)  ACj |x8  large-set  Kasami  with  PTSC  =  227040.  (b)  Optimal  <S°5*  8 
designed  under  Underloaded  Case  2a  with  PTSC  =  (8)2(15)3  =  216000. 
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Fig.  8.  Optimal  signature  set  iSgf  x 42  designed  under  Overloaded  Case  2  with 
minimum  PTSC  =  (42)2(31)3  +  4(31)(30)  =  52555044. 
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IV.  Rank-2-optimal  Adaptive  Design  of  Binary 
Spreading  Codes 

This  work  has  been  published  in  the  IEEE  Transactions  on  Information  Theory, 
Sept.  2007. 

The  performance  of  direct-sequence  code-division-multiple-access  (DS-CDMA) 
systems  is  determined  by  the  set  of  user  spreading  codes  in  conjunction  with 
channel  and  receiver  design  specifics.  In  recent  literature,  several  methods  have 
been  presented  for  the  design  of  real/complex- valued  [1]- [14]  or  binary  [8],  [15]- 
[24]  sets  of  DS-CDMA  spreading  codes  under  various  optimization  criteria. 
Among  them,  in  [l]-[6],  [14]- [22]  the  spreading  code  set  is  treated  as  a  sin¬ 
gle  matrix  parameter  to  be  optimized  while  the  works  in  [7]- [12],  [23],  [24] 

present  iterative  algorithms  that  update  the  codes  of  the  set  individually  in  a 
round-robin  fashion. 

In  particular,  sets  of  spreading  codes  that  minimize  the  total  squared  cor¬ 
relation  (TSC)  were  designed  in  [1],  [2],  while  distributed  algorithms  that  it¬ 
eratively  decrease  TSC  by  updating  one-by-one  the  individual  codes  of  the  set 
were  developed  in  [7]- [10].  Band-limited  sets  that  maximize  the  user  capacity 
of  synchronous  DS-CDMA  systems  were  constructed  in  [3].  Optimal  sets  for 
asynchronous  DS-CDMA  systems  were  designed  in  [4]- [6],  while  the  iterative 
method  of  [9],  [10]  was  modified  to  suit  asynchronous  systems  in  [11]  and  multi- 
path  channels  in  [12].  Collaborative  multibase  designs  were  studied  in  [13].  In 
all  the  above  developments,  each  user  spreading  code  was  allowed  to  take  any 
value  in  the  real  vector  space  subject  to  a  (unit-)norm  constraint. 

In  [17],  we  concentrated  on  binary  sets  of  spreading  codes  and  derived  new 
lower  bounds  on  the  TSC.  Optimal  binary  sets  that  meet  the  TSC  bound  with 
equality  were  also  designed  in  [17],  as  well  as  in  [18], [19],  while  their  sum 
capacity,  total  asymptotic  efficiency,  and  maximum  squared  correlation  were 
evaluated  in  [20] .  In  addition,  a  searching  algorithm  for  minimum-TSC  spread¬ 
ing  code  sets  with  low  cross-correlation  spectrum  was  presented  in  [21].  Inter¬ 
estingly,  if  we  consider  binary  sets  of  spreading  codes  for  transmission  over  a 
synchronous  DS-CDMA  channel  with  unequal  channel  gains  for  different  users, 
then  the  minimum-TSC  sets  of  [17]- [19]  also  minimize  the  corresponding  total 
weighted  squared  correlation  (TWSC)  [14]  as  long  as  the  system  is  underloaded 
and  the  spreading  code  length  (processing  gain)  is  odd  or  a  multiple  of  four. 

For  overloaded  systems,  however,  with  unequal  user  gains,  the  minimum- 
TWSC  optimality  of  minimum-TSC  binary  sets  is  lost;  the  design  of  optimal 
binary  sets  for  such  systems  is  still  an  open  problem.  In  an  attempt  to  design 
appropriate  binary  spreading  codes  for  overloaded  synchronous  transmissions 
with  unequal  received  power  values  as  well  as  asynchronous  transmissions  over 
-potentially-  multipath  channels,  suboptimal  distributed  algorithms  were  pre¬ 
sented  in  [8],  [23],  [24].  In  these  works,  the  user  codes  are  updated  one-by-one, 
similar  to  the  approach  followed  in  [7]- [12]  for  real- valued  spreading  codes,  sim¬ 
plifying  the  set  design  problem  to  the  adaptive  design  of  one  spreading  code  in 
the  presence  of  colored  interference. 
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In  this  report,  we  revisit  the  work  of  [8],  [23],  [24]  and  consider  again  the 
adaptive  design  of  binary  spreading  codes  that  -in  the  presence  of  colored 
interference-  maximize  the  signal-to-interference-plus-noise  ratio  (SINR)  at  the 
output  of  the  maximum-SINR  (MSINR)  filter.  The  optimal  code  is  a  function 
of  the  disturbance  (interference  plus  noise)  autocovariance  matrix  and  its  eval¬ 
uation  over  the  binary  field  is  NP-harcl  [24].  Instead,  in  this  present  work  we 
propose  to  eigendecompose  and  approximate  the  inverse  disturbance  autoco¬ 
variance  matrix  by  its  two  maximum-energy  components  alone.  Then,  we  show 
how  to  optimize  the  binary  spreading  code  under  the  rank-2  approximation  of 
the  inverse  disturbance  autocovariance  matrix  with  lower  than  quadratic  com¬ 
plexity.  We  demonstrate  the  significant  SINR  performance  superiority  of  the 
proposed  rank-2-optimal  adaptive  design  in  comparison  to  direct  hard-limiting 
of  the  maximum  eigenvector  of  the  inverse  disturbance  autocovariance  matrix 
(or  the  minimum  eigenvector  of  the  interference  autocovariance  matrix)  [8] ,  [24] , 
which  we  see  to  be  equivalent  to  rank-l-optimal  adaptive  binary  spreading  code 
design.  Moreover,  simulation  studies  indicate  that  the  proposed  rank-2-optimal 
spreading  code  can  be  the  exact  full-rank-optimal  solution  with  significantly 
higher  probability  than  the  rank-l-optimal  code  (for  example,  when  the  re¬ 
turned  rank-l-optimal  code  is  globally  optimal  with  probability  0.4,  the  rank- 
2-optimal  code  is  globally  optimal  with  probability  0.85).  In  fact,  simple  it¬ 
erative  Hamming-distance-1  steepest  descent  search  [22],  [24]  initialized  at  the 
proposed  rank-2-optimal  spreading  code  raises  the  probability  of  convergence 
to  the  full-rank-optimal  solution  to  as  high  as  0.95.  Certainly,  the  proposed 
adaptive  binary  code  design  can  serve  as  an  appropriate  initialization  point  for 
other,  potentially  more  sophisticated,  iterative  binary  search  procedures. 

A.  Signal  Model  and  Problem  Statement 

Consider  the  vector  signal  model 

r  =  b\f~Ps  +  y  (71) 

where  6e{±l}isa  uniformly  distributed  bit  random  variable,  s  £  RL  a  deter¬ 
ministic  vector  waveform  (spreading  code)  with  ||s||  =  1  such  that  all  collected 
energy  is  absorbed/represented  by  the  scalar  P  >  0,  and  y  €  RL  is  a  zero- 
mean  random  disturbance  vector  with  positive  definite  autocovariance  matrix 
Ft  =  if{yyT}  (-E{-}  denotes  statistical  expectation  and  T  is  the  transpose 
operator). 

The  general  signal  model  of  (71)  covers,  for  example,  synchronous  DS- 
CDMA  transmissions  where  a  user  of  interest  with  spreading  code  s  transmits 
over  an  additive  noise  channel  in  the  presence  of  K  interfering  users.  In  that 
case,  the  overall  additive  disturbance  term  takes  the  form  y  =  Y^k= l  ^k\fPk^k  + 
n  where  bk  £  {=h  1}  ,  Pk  >  0,  and  s k  £  RL  are  the  uniformly  distributed  user 
bit,  received  energy  per  bit,  and  normalized  deterministic  spreading  code  of  the 
fcth  interferer,  k  =  1,2,  and  n  represents  additive  zero-mean  channel 

noise.  Similarly,  (71)  can  model  asynchronous  DS-CDMA  transmissions  when 
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the  receiver  is  synchronized  with  the  signal  of  interest  or  can  be  extended  to 
cover  multipath  transmissions  if  s  is  replaced  by  the  effective  (channel  processed) 
user  signature. 

For  an  arbitrary  spreading  code  s,  the  linear  filter/receiver  w  that  exhibits 
maximum  SINR  at  its  output  has  the  form 

w(s)  =  cR^1s,  c  >  0,  (72) 

and  the  (maximum)  output  SINR  value  is 


SINR  (s) 


(wT  (s)  b\/Ps 

>1 

£ 

;{(wT  (s)z)2] 

=  PsTR_1s. 


(73) 


We  recall  that,  if  the  additive  disturbance  vector  y  is  Gaussian,  then  sgn  [wT  (s)  r] 
is  the  minimum  bit-error-rate  (BER)  detector  with  error  probability  given  by 

BER(s)  =  Q  (^\/PsTH^1sJ  where  Q  (a)  =  —=e~ti dt. 

Our  objective  is  to  design  the  spreading  code  s  so  that  the  corresponding 
SINR  (s)  value  is  maximized  (and  BER  (s)  is  minimized  under  the  Gaussian 
additive  disturbance  assumption) .  Since  the  disturbance  autocovariance  matrix 
R  is  positive  definite, 

L 

R  =  ^  Ajqjqf,  Ai  >  A2  >  . . .  >  AL  >  0,  ||q»||  =  1,  *  =  1, 2, . . . ,  L,  (74) 
*= 1 

represents  its  eigendecomposition  where  A;  and  q;  are  the  ith  eigenvalue  and 
normalized  eigenvector,  respectively,  of  R.  Then,  the  real-valued  spreading  code 
s  that  maximizes  SINR(s)  is  given  by  [10]- [12],  [24] 

sr-opt  =  arg  max  {PsTR_1s|  =  qL.  (75) 

s(EMl  ,||s||=1 

Therefore,  the  optimal  real-valued  code  sr_opt  can  be  obtained  with  complex¬ 
ity  that  is  dominated  by  the  complexity  of  the  eigendecomposition  of  the  L  x  L 
matrix  R. 

In  this  present  work,  we  are  interested  in  maximizing  SINR  (s)  subject  to 
the  constraint  that  s  is  binary.  The  problem  of  obtaining  the  optimal  binary 
spreading  code, 

sopt  =  arg  max  {PsTR_1sj,  (76) 

is  NP-hard  [23] ,  [24]  and  can  be  solved  through  exhaustive  search  over  all  possi¬ 
ble  L-bit  combinations.  In  the  next  section,  we  seek  a  computationally  efficient 

alternative  for  the  design  of  a  spreading  code  s  G  j  ±  -4=  j  . 
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B.  Rank-2-optimal  Design  of  Binary  Spreading  Codes 

Given  the  eigendecomposition  of  R  in  (74), 

R  1  =  V'  T~  >  y^—  >  •■•  >  y-  >  0,  ||qi||  =  1,  i  =  1, 2, . . . ,  L, 

M  Al  Al-1  a  1 

1=1 

(77) 

and  the  binary  spreading  code  optimization  problem  in  (76)  can  be  rewritten  as 
sopt  =  arg  max  \  V  -j-  (sTq.i)“  1  (78) 

l=A‘  J 

where  0  <  (sTq;)2  <  1,  i  =  1,2 Therefore,  the  optimal  binary  vector 
sopt  maximizes  the  sum  of  its  weighted  projections  on  the  eigenvectors  q;  with 
weights  1/A  j,  i  =  l,...,L. 

If  we  simplify  the  optimization  problem  by  keeping  only  the  strongest  term 
Yj-  (sTq^)  in  (78)  (which  is  equivalent  to  using  the  approximation  R-1  ~ 
Yj-q^qK  in  (76)),  we  obtain  the  rank- 1- optimal  spreading  code 

si  =  arg  max  ( (sTqL)“|  =  ±sgn  (qL)  (79) 


where  sgn  (•)  denotes  the  sign  operator.  Hence,  we  showed  that  the  rank-1- 
optimal  binary  code  is  simply  the  minimum-eigenvalue  eigenvector  of  R  passed 
through  a  sign  hard-limiter.  From  another  point  of  view,  if  we  follow  the  com¬ 
mon  approach  of:  (i)  first  relaxing  the  binary  constraint  by  allowing  s  £ 

||s||  =  1,  (ii)  then  solving  the  relaxed  problem  to  obtain  sr_opt  =  qL,  and 
finally  mapping  sr_opt  =  qL  to  the  binary  field  by  taking  the  sign  of  its 
components,  then  we  again  obtain  the  rank-l-optimal  binary  spreading  code 
sgn(sR_opT)  =  sgn(qi)  =  Sp  Therefore,  the  relaxation  approach  in  [8],  [24] 
is  equivalent  to  rank-1  approximation  of  R-1. 

To  obtain  a  binary  spreading  code  with  higher  SINR  than  SINR(si) ,  in 
this  present  work  we  propose  to  keep  the  two  strongest  terms  (sTqL)“  + 
^  1  (sTqT-i)2  in  (78)  or  -equivalently-  use  rank-2  approximation  of  R_1  ~ 
Yj-qLqJ  +  YyyQL-iqL-i-  Then,  our  optimization  problem  becomes 


Below,  we  show  that  the  rank-2-optimal  binary  vector  S2  in  (80)  can  always  be 
obtained  with  less  than  quadratic  complexity  (less  than  the  complexity  required 
for  the  eigendecomposition  of  R) . 

We  begin  by  defining  the  complex  vector 


A 

z  = 


(81) 
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Then,  the  problem  in  (80)  is  converted  to 


S2  =  arg  max 


sel  ± 


JlV 

vrj 


sTz|“  >  =  arg  max  ||sTz| 


(82) 


In  an  effort  to  solve  (82)  in  less  than  quadratic  complexity,  we  consider  the 
auxiliary-variable  techniq  ue  that  has  been  used  in  a  communications  theory 
context  before  in  [25]  and  then  again  in  [26].  Let  s  =  [si  S2  s.l]T  and 


z  =  [zt  z2 


.  Zl]  where 

Zi  =  \zi\e3<ki, 


7T  3?r  ,  ^ 

~  <  4>i  <  ?  i=  1, 2, . . . ,  L. 


(83) 


We  introduce  an  auxiliary  variable  <j>  £  [— 7r,  7 r)  and  rewrite  the  quantity  to  be 
maximized  as 


s  z  =  max 

</>G  [— 7T,7r' 


(Re  { 


sTze  — 


}}  =  max  <  sz  \zi\ cos  {4>  -  fa)  >  (84) 

</>£[— 7T, 7 r)  1  A J  1 


where  Re  {■}  extracts  the  real  part  of  a  complex  number.  Then, 


max  ||sTz|}=  max  max  <  Si  IzA  cos  ( <f)  —  fa) 

-  L  f  0G[— 7T, 7r)  I  ^ J 


sG 


sG 


=  max 

0G[— 7T,7r) 


^2\zi\  max  {si  cos  ((j>-  fa)} 
i=i  ^{±7?} 


(85) 


=  max 

0G[—  7T,7r) 


W  \cos((j)-fa)\ 


(86) 


since  for  any  </>  £  [— k,  n)  the  optimal  s,  value  in  (85)  is  s*  (0)  =  -isgn  (cos  (</>  —  fa)) . 
i  =  1, 2, . . . ,  L.  The  final  quantity  J2i=  i  \zi\  ^  lcos  (</>  ~  <&)|  in  (86)  is  maxi¬ 
mized  for  a  particular  value  <)>opt  £  [— 7r,  tt)  and  s  (</>opt)  =  [si  (^opt)  s2  (<Poy>t) 

...  Sl  (0opt)]  is  the  rank-2-optimal  binary  vector  we  are  searching  for  in  (12), 
i.e.  S2  =  s  (4>opt)  •  We  will  now  show  that  we  can  always  construct  a  set  of  L 

spreading  codes  U  =  {ui,  u2, . . . ,  u^}  ,  U;  £  j±^=j  ,  and  guarantee  that 

s  (^OPt)  £  lA.  Therefore,  the  maximization  in  (12)  can  be  carried  out  over  a 
set  of  L  candidates  only  without  loss  of  optimality. 

Partition  =  {1, 2 , ,L}  into 


h  =  { i  : 


7T  7T\  1 

"2’ 27/ 


and  I-  <  i  :  fa  £ 


5-  ¥  I  !•  =  -  n  (87) 


and  define  the  set  of  angles 

7  A 


fa,  i  €  h 

-  7T,  i  £  I2 


*=1,2, ...,L, 


(88) 
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such  that  <  fa  <  i  =  1, 2, . . . ,  L.  Define,  for  notational  simplicity,  the  vec¬ 
tor  operation  s  (</>)  =  [si  (</>)  s2  (</>)  Sl  (</>)] T  with  §i  (</>)  =  sgn  ^cos 
i  =  1,  2, . . . ,  L,  (f>  £  [—7 r,  7r) .  Then, 


i  £  I\ 


Consider  a  mapping  e  from  Z^  to  Z^  that  sorts  </>i,  cf> 2, . . . ,  4>l-  —  §  <  ^e(i)  < 
</>e(2)  <  ■  •  ■  <  <j>e(L)  <  f-  In  (85),  (86),  maximization  with  respect  to  </>  can  be 
carried  out  over  any  subinterval  of  length  7r.  We  choose  4>e(i)  —  f  and  +  f 
as  the  limits  of  such  a  subinterval  and  rewrite  the  original  optimization  problem 
of  (85)  as 


max  {|sTz|}=  max  1  /J  I2*  I  max  {si  cos  (0  —  </>j)}  >  . 

sC  |  =h  ~^=  |  L  0e[0e(l)-f  ,0e(l)  +  f )  ^  i=l  j 

(9°) 

By  examining  the  subintervals  0e(i)  -  f,^e(2)  -  f)  ,  4>e(2)  -  §,</>e(3)  -  f)  , 

•  •  ■ ,  $e(L-i)  -  f ,  4>e(L)  ~  f )  ,  4>e(L)  -  f ,  </>e(i)  +  f ) ,  we  conclude  that 

[®e(l)  ('5-0  ®e(2)  (0)  ■  •  ■  ^e(L)  (</*)] 


Pe(3)  2 


[+1  -1  -1  . 

•  •  - 1  - 1] , 

0e(l)  - 

[+1  +1  -1  . 

•  •  - 1  - 1] , 

0e(2)  - 

[+1  +1  +1  . 

•  •  +  1  —  1] ) 

fie(L-l) 

[+1  +1  +1  . 

..  +  1  +1], 

We  collect  the  L  binary  vectors  that  appear  in  the  L  cases  above, 

U;  =  [+1  . . .  +  1  —1  ...  —  1]T,  i  =  1, 2, . . . ,  L,  and  create  the  matrix 


U  =  [ui  u2  ...  u L]  (92) 

whose  itli  row  we  denote  by  d^,  i  =  1,2 ,...,L.  Then,  we  reorganize  U  to 

Af-  •  -  1T 

U  =  [ui  u2  ...  ul]  =  di  d2  ...  di  by  defining  the  binary  vectors 

di  de— 1  (*)  ?  ^  1)2,...,  Z/,  (98) 

where  e_1  :  Z /,  — ■>  Z^  is  the  inverse  sorting  mapping  (note  that  s  (</>)  G 
{up  u2, . . . ,  uL}  for  any  <j>  <E  Ue(1)  -  f ,  <^e(i)  +  f  ))■  Finally,  we  define 


A  |  *  e  h 

I  -4-di,  i  G  h 


,  i  1,2,...,  L, 
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and  construct 


U  =  [ui  u2  ...  uL]  =  [di  d2  ...  dL]T  .  (95) 


The  set  U  =  {ui,  u2, . . . ,  u^}  contains  s  for  any  (f>  £  |  <pen)  —  ^ ,  <pem  ,  2 

0OPT  €  |pe(i)  _  f  ,<f>e(i)  +  f  J>  which  implies  s  (</>opt)  =  s2  <E  U.  Hence, 
becomes 

s2  =  arg  max  { |  sTz  | }  . 


.  But 
(82) 

(96) 


We  conclude  that  the  maximization  task  in  (82)  has  been  converted  to  an 
equivalent  linear-complexity  maximization  problem  in  (96).  The  complexity  of 
the  construction  of  U  is  dominated  by  the  complexity  of  the  sorting  function  e 
which  is  of  order  O  ( L  log2  L) .  Therefore,  we  have  described  a  method  to  obtain 
the  rank-2-optimal  binary  vector  s2  in  (80)  with  complexity  that  is  dominated 
by  the  complexity  of  the  eigendecomposition  of  the  L  x  L  matrix  R  alone. 

As  a  brief  summary,  the  sequence  of  the  proposed  calculations  is  as  follows. 
The  autocovariance  matrix  R  is  eigendecomposed  as  in  (74)  followed  by  the 
computation  of  the  z  vector  by  Eq.  (81).  The  sets  I\  and  /2  are  constructed 
through  (87);  the  angles  ■  ■  ■  ,<I>l  are  calculated  by  (88)  and  sorted  to 

<  4>e{2)  <  •  ■  •  <  <t>e(L)-  Next,  U  is  constructed  as  in  (92)  and  reorganized 
to  U  via  (93);  then,  U  is  found  by  (95).  The  columns  of  U,  ui,  u2, . . . ,  Ul, 
are  used  to  calculate  the  quantities  |  uf  z  |  ,  |  u[T z  |  , . . . ,  |  u£  z  |  that  are  compared 
to  each  other  to  identify  the  maximum  value  among  them,  say  |  u J z  |  for  some 
j  £  {1,2,  The  latter  determines  the  rank-2-optimal  binary  spreading 

code  s2  =  Uj. 


C.  Simulation  Studies 

We  consider  a  DS-CDMA  system  where  the  user  of  interest  transmits  over  an 
additive  zero-mean  white  Gaussian  noise  channel  with  variance  a2  in  the  pres¬ 
ence  of  K  synchronous  interferers.  The  system  processing  gain  (spreading  code 
length)  is  L  =  16.  The  received  signal-to-noise  ratio  (SNR)  of  the  user  of  inter¬ 
est,  SNR  =  -Sr,  is  set  to  lOdB,  while  the  received  SNRs  of  the  K  interferers, 

SNR.;.  =  -pf-,  k  =  1, 2, . . . ,  K,  are  uniformly  spaced  between  8dB  and  lldB.  The 
interfering  spreading  codes  are  randomly  generated. 

In  our  studies,  we  compare  the  SINR  performance  of:  (i)  the  optimal  bi¬ 
nary  spreading  code  sopt  of  (76)  obtained  through  exhaustive  search  over  all 
possible  L-bit  combinations;  (ii)  the  rank- 1- optimal  binary  spreading  code  si  of 
(79)  obtained  by  applying  the  sign  operator  on  the  minimum-eigenvalue  eigen¬ 
vector  of  the  interference-plus- noise  autocovariance  matrix  [8],  [24];  and  (in) 
the  proposed  rank- 2- optimal  binary  spreading  code  s2  of  (80)  obtained  by  the 
procedure  developed  in  Section  III.  For  comparison  purposes,  we  evaluate  the 
SINR  loss,  SINR(sopt)  —  SINR(s),  of  s  =  si  or  s2  with  respect  to  the  SINR 
of  the  optimal  binary  spreading  code  sopt-  The  results  that  we  present  are 
averages  over  1,  000  randomly  generated  interference  signature-set  realizations. 
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In  Fig.  9,  we  plot  the  SINR  loss  of  the  rank-l-optimal  and  rank-2-optimal 
binary  spreading  codes  as  a  function  of  the  number  of  interferers  K.  Since  we  are 
particularly  interested  in  overloaded  systems,  we  vary  K  from  16  to  40  interfer¬ 
ers.  We  observe  that  the  proposed  rank-2-optimal  spreading  code  exhibits  less 
than  0.04dB  performance  loss  which  is  significantly  lower  than  the  performance 
loss  of  the  rank-l-optimal  code  (interestingly,  both  loss  values  decrease  as  K 
increases) . 

In  Fig.  10,  we  plot  the  probability  of  global,  full-rank,  optimality 
Pr  {s  =  sopt}  for  the  rank-l-optimal  and  rank-2-optimal  binary  spreading  codes 
as  a  function  of  the  number  of  interferers  K.  We  observe  that 
.23  <  Prjsi  =  sopt}  <  -43,  while  .75  <  Pr{s2  =  sopt}  <  -86  as  K  is  varied 
from  16  to  40  interferers.  Therefore,  with  the  proposed  optimization  of  the  bi¬ 
nary  spreading  code  under  the  rank-2  approximation  of  the  inverse  disturbance 
autocovariance  matrix,  we  have  practically  doubled  the  probability  that  the 
designed  spreading  code  is  full-rank  optimal  with  only  O  ( L  log2  L)  additional 
computational  cost. 

To  the  extend  that  neither  the  rank-l-optimal  nor  the  rank-2-optimal  binary 
spreading  code  design  is  globally  optimal  with  probability  one,  the  returned 
codes  can  potentially  be  used  as  the  initialization  point  of  an  iterative  steepest 
descent  search  [22],  [24]  that  may  converge  to  a  (suboptimal  in  general)  spread¬ 
ing  code  with  higher  SINR.  In  this  spirit,  we  repeated  the  studies  of  Figs.  9 
and  10  and  fed  the  returned  codes  Si  and  s2  to  a  Hamming-distance-1  steepest 
descent  search  to  produce  the  improved  spreading  codes  s^sd  and  s2jsd,  cor¬ 
respondingly.  In  Fig.  11  we  plot  the  SINR  loss  of  S13D  and  s2jsd;  hr  Fig.  12 
we  plot  the  corresponding  probabilities  of  full-rank  optimality.  The  SINR  su¬ 
periority  of  the  rank-2  approach  is  maintained.  Interestingly,  we  observe  that 
Pr{s2,SD  =  sopt}  >  -9  for  any  interference  load.  For  example,  when  K  =  20 
the  probability  of  global  optimality  is  .95  for  s2isd  (and  .79  for  s^sd)-  To  ex¬ 
amine  the  speed  of  convergence  of  the  steepest  descent  search,  in  Fig.  13  we 
plot  the  evolution  of  the  SINR  value  as  a  function  of  the  number  of  iterations. 
The  SINR  of  the  globally  optimal  code  sopt  is  included  as  reference.  We  ob¬ 
serve  that  when  the  search  is  initialized  at  the  proposed  rank-2-optimal  code, 
convergence  is  achieved  within  three  iterations  only. 

D.  Conclusions 

In  the  broad  context  of  spread-spectrum  communications  (for  example  CDMA 
systems)  or  chip-based  signal  waveform  design,  we  considered  the  problem  of 
adaptively  identifying  the  binary  code  that  maximizes  the  signal-to-interference- 
plus-noise  ratio  (SINR)  at  the  output  of  the  maximum-SINR  filter.  The  optimal 
code  is  a  function  of  the  disturbance  autocovariance  matrix  and  the  optimiza¬ 
tion  problem  is  NP-hard.  Instead,  we  developed  a  new  algorithm  of  less  than 
quadratic  complexity  for  the  computation  of  the  optimal  code  under  rank-2 
approximation  of  the  inverse  disturbance  autocovariance  matrix. 

As  an  illustration  of  practical  significance,  we  demonstrated  the  great  SINR 
performance  improvement  of  the  proposed  rank-2-optimal  adaptive  design  for 
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overloaded  CDMA  systems  with  unequal  user  power  over  the  rank-l-optimal 
adaptive  design  which  -as  we  showed-  is  equivalent  to  direct  hard-limiting  of 
the  minimum-eigenvalue  eigenvector  of  the  disturbance  autocovariance  matrix. 
In  fact,  the  proposed  rank-2-optimal  design  practically  doubles  the  probability 
that  the  returned  code  is  full-rank-optimal  and  when  the  rank-2-optimal  design 
is  used  to  initialize  trivial  Hamming-distance-1  steepest  descent,  the  full-rank 
optimal  code  is  reached  with  probability  greater  than  0.9  (within  three  itera¬ 
tions).  Certainly,  the  proposed  adaptive  binary  code  design  may  serve  well  as 
the  initialization  point  for  other,  potentially  more  sophisticated,  iterative  binary 
search  procedures. 
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SINR  Loss  (dB) 


Number  of  Interferers 


Fig.  9:  SINR  loss  of  rank-l-optimal  and  rank-2-optimal  binary  spreading  code 
designs  versus  number  of  interferers  (signature  length  L  =  16). 
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Probability  of  full-rank  optimality 


Fig.  10:  Probability  of  full-rank  optimality  of  rank-l-optimal  and  rank-2-optimal 
binary  spreading  code  designs  versus  number  of  interferers. 
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SINR  Loss  (dB) 


Fig.  11:  SINR  loss  of  steepest  descent  upon  convergence  with  rank-l-optimal 
and  rank-2-optimal  initialization. 
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Probability  of  full-rank  optimality 


Number  of  Interferers 


Fig.  12:  Probability  of  full-rank  optimality  of  steepest  descent  search  upon 
convergence  with  rank-l-optimal  and  rank-2-optimal  initialization. 
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SINR  (dB) 


Fig.  13:  SINR  of  steepest  descent  with  rank-l-optimal  and  rank-2-optimal  ini¬ 
tialization  versus  number  of  iterations. 
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V.  Scalable  TSC-optimal  Overloading  of  Binary 
Signature  Sets 

The  work  has  been  presented  at  IEEE  GLOBECOM  2006. 

In  code-division  multiplexing  applications,  for  example  direct-sequence  code¬ 
division  multiple-access  (DS-CDMA)  cellular  communications  systems,  each  of 
the  K  participating  signals/users  is  equipped  with  a  unique  identifying  signature 
vector  Sfe  G  CL,  ||sfc||  =  1,  k  =  1,2,..., K.  All  signatures  put  together  in  the 
form  of  a  matrix  define  what  we  call  the  signature  set 

S  =  [sr  S2...S*-]  eCLxK.  (97) 

In  synchronous  code-division  multiplexing  transmissions  over  well  behaved 
Nyquist  channels,  we  are  interested  in  using  a  signature  set  with  the  smallest 
possible  total  squared  correlation  (TSC)  value  [1]- [6] 

TSC(S)  =  EXXSE  (98) 

»=  1  3= 1 

where  H  denotes  the  Hermitian  operator.  For  complex/real-valued  signature 
sets  5  G  CLxK  or  RLxK ,  if  I<  >  L,  TSC(S)  >  ^  [1]  (of  course,  TSC(S)  >  K 
if  K  <  L).  Overloaded  (K  >  L)  sets  with  TSC  equal  to  have  been  known 
as  Welch-bound-equality  (WBE)  sets.  Algorithms  and  studies  for  the  design 
of  complex  or  real  WBE  signature  sets  can  be  found  in  [3]- [10].  While  not  an 
issue  for  underloaded  systems,  it  is  well  understood  that  as  new  users  enter/exit 
an  overloaded  system  the  existing  signature  set  has  to  be  re-designed  or  re-cast 
to  maintain  TSC  optimality  (WBE  sets  do  not  stay  optimal  if  signatures  are 
added/removed).  As  a  solution  to  the  scalability  problem  in  the  complex/real 
vector  domain,  two  subclasses  of  WBE  sets  are  proposed  in  [11]  that  provide 
satisfactory  TSC/signal-to-interference-plus-noise  ratio  (SINR)  performance  as 
the  number  of  users  changes  dynamically. 

In  digital  transmission  systems  it  is  necessary  to  have  finite-alphabet  signa¬ 
ture  sets.  Recently,  new  bounds  were  derived  on  the  TSC  of  binary  antipodal 
signature  sets  together  with  optimal  designs  for  arbitrary  signature  lengths  and 
set  sizes  [2] -[4].  Studying  the  designs  provided  therein  we  observe  that,  similar 
to  the  complex/real-field  case,  as  long  as  the  system  is  underloaded  users  can  be 
accommodated  in  and  out  without  having  to  change  the  existing  signature  set 
to  maintain  TSC  optimality.  In  the  overloaded  regime,  when  additional  users 
enter/exit  the  system  the  existing  signature  set  has  to  be  re-designed/re-cast.  In 
this  paper,  we  present  a  complete  solution  to  the  scalability  problem  for  binary 
sets  with  overloading  up  to  100%. 

Consider  a  (K  =  L  +  n,L),  n  G  {—1,0, 1,2},  TSC-optimal  fixed  initial 
matrix  5/  =  [si  s2  . . .  sL+„] , 

sfc  G  j±-^=  j  ,  k  =  1, 2, . . . ,  L  +  n,  n  G  {—1,  0, 1,  2},  given  by  [2]- [4] .  To  this 
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set,  we  will  append  N,  1  <  N  <  L  +  2,  carefully  designed  unique  signatures 
{sE+„_|_fc}^1 ,  n  £  {—1,  0, 1,  2},  that  define  the  extended  signature  matrix 


A  f  1  1  N 

Se  =  [Si+n+l  •  •  .  S£+n+jv]  £  ^=J 


(99) 


The  complete  ( K  =  L  +  n  +  N,L)  signature  matrix  S  can  then  be  expressed  as 


(  -i  -Lx  (L+n+iV) 

S  =  [S,  S,,\  £  |±^=j  .  (100) 


From  (98)  and  (100), 

TSC(S)  =TSC(5/)+TSC(5b)+F(5/,5b)  (101) 

where  F (Si,Se)  is  the  sum  of  the  absolute  squared  cross-correlations  between 
every  signature  in  set  Si  and  every  signature  in  set  Se, 

F(5/,5£)=2^|sfs,|2  (102) 

i  £  { 1 ,2, . . . ,L-\~ny 

j  €  {Z,  +  n  +  l,...,L  +  n  +  iV} 

(T  denotes  the  transpose  operator).  To  provide  scalable  signature  set  designs 
and  control  the  TSC  value,  in  this  paper  we  consider  the  problem  of  minimizing 
the  TSC  of  S  =  [5/  Se]  subject  to  Si  being  fixed  and  TSC-optimal.  In  past 
related  literature,  Si  and  Se  were  designed  to  be  orthogonal  binary  Hadamard 
sets  (which  requires  the  signature  length  I  to  be  a  multiple  of  four)  scrambled 
independently  by  complex  random  PN-sequences  [15].  The  mutual  interference, 
quantified  as  a  loss  in  SINR,  between  the  two  sets  of  signatures  was  suppressed 
through  iterative  detection.  In  [11]  it  was  shown  that  for  a  fixed  and  optimal 
S/,  TSC(S)  >  L  +  3N  for  all  S  £  CLx(L+N\  1  <  N  <  L.  Here,  we  derive  new 
tighter  bounds  on  the  conditional  TSC  of  binary  sets  for  all  overloaded-by-lV 
cases  where  1  <  N  <  L  +  2.  The  work  can  be  easily  extended  to  cover  higher 
overloading  (from  a  practical  perspective,  however,  overloading  by  more  than 
100%  may  not  be  of  much  interest  at  present) .  In  addition,  we  provide  optimal 
constructions  that  meet  the  new  conditional  TSC  bounds  with  equality  (hence, 
establish  their  tightness).  Included  numerical  comparison  studies  show  that  the 
proposed  sets  minimally  exceed  the  global  TSC  bound  of  [2]  and,  in  fact,  in 
certain  cases  meet  the  bound  with  equality  and  become  globally  TSC  optimal. 


A.  New  Bounds  on  the  Conditional  TSC  of  Overloaded 
Binary  Sets 

Our  goal  is  to  derive  lower  bounds  on  (101)  subject  to  Si  being  fixed  and 

LxN 


optimal  over  all  Se  € 


{H 


—  Si  >  where  1  <  N  <  L  +  2.  We  derive 


the  (conditional  TSC)  bounds  on  a  case- by-case  basis. 
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Case  1:  L  =  0  (mod  4),  1  <  N  <  L. 

Let  Si  be  an  optimal  (L,L)  initial  matrix  with  TSC(<S/)  =  L  as  observed  from 
Table  I4.  From  [2],  the  set  S  =  [5/  s]  is  optimal  for  any  s  €  {±-^=}L  and 


TSC(S')  -TSC(5/)  -  ||  s  || 2 


(L+l)2 


L  —  l 


-L-  1  =  2; 


(103) 


hence,  for  any  arbitrary  extension  matrix  Se  €  |±-^= j  ,  F(Si,Se) 
From  Table  I,  TSC(<Se)  >  N  and  we  calculate 


=  2N. 


TSC(S)  >L  +  3N. 


(104) 


Case  2:  L  =  2  (mod  4),  1  <  N  <  L  +  2. 

Let  Si  be  an  optimal  ( L  +  2,  L)  initial  matrix  with  TSC(<S/)  =  as  calcu¬ 

lated  from  Table  I.  From  [2],  the  set  S  =  [5/  s]  is  optimal  for  any  s  £  {±-l=}L 
and 


TSC(S')  -  TSC(Sj)  -  ||s||2  =  ^(L  +  3)  +  ^  -  {L+l2)  -  1  =  2  +  i; 

(105) 


hence,  for  any  arbitrary  extension  matrix  Se  £  >  F (Si.  Se) 

2 N  (^x^)-  From  Table  I, 


TSC(Sb)  > 
and  we  calculate 


N  +  2{N~il)2 ,  N  =  1  (mod  2) 
TV  +  2N^~1\  N  =  0  (mod  2) 


(106) 


TSCLS1  >  I  +  N  =  1  (mod  2) 

\  +N+  2N{^1]  +  2jv^+2) ,  N  =  0  (mod  2). 


(107) 


Case  3:  L  =  3  (mod  4),  1  <  IV  <  L  +  1. 

Let  Si  be  an  optimal  ( L  +  1  ,L)  initial  matrix.  By  Table  I,  TSC(<S/)  = 
By  [2],  S  =  [5/  s]  is  optimal  for  any  s  £  {±^}L  and 


TSC(S')  -  TSC(5/)  -  ||s||2  =  ^ L+l2)  +  ^j-^J  -  {L+Ll)  -  1  =  2  +  |; 

(108) 

4Table  I  of  [2]  is  reproduced  for  convenience  in  Table  I  herein. 
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r  x  LxN 

hence,  for  any  arbitrary  extension  matrix  Se  G  |±^=|  ,  F (Si,Se)  = 

2 N  From  Table  I, 


TSC(Sij)  >  N  - 


N(N  —  1) 

Z2 


and  we  calculate 


TSC(S)  > 


(L  +  l)2 


■  N  ■ 


N(N-l)  i  2N(L  +  1) 

Z2  +  Z  ' 


(109) 

(110) 


Case  4:  L  =  1  (mod  4),  1  <  N  <  L  —  2. 

Let  Si  be  an  optimal  ( L  —  1 ,  L)  initial  matrix5.  By  Table  I,  TSC(<S/)  = 
L  —  1  +  _  Sj  can  be  obtained6  from  an  optimal  (L  —  1,L  —  1)  set 

S  ,  i.e.  Si  contains  an  optimal  (hence,  orthogonal  [2])  submatrix  S  of  size 
(L  —  1,  L  —  1).  This  property  is  utilized  in  the  following  Lemma  and  Theorem, 
which,  in  turn,  are  used  to  derive  the  bound. 


Lemma  1:  Let  S  = 


S1  s2  ■  •  •  SL— 1 


{±^} 


(Z--l)x(L-l) 


,  L  =  1  (mod  4), 


f  1  L—l 

be  the  optimal  (orthogonal)  submatrix  of  Si.  Then,  for  any  x  €  j  , 


£;=1  x  si  e  (V--V)- 


'  L—l  L  —  l  \ 

L 

Proof:  The  vectors  sx , . . . ,  sL_1  form  an  orthogonal  basis  for  RL_1.  Then,  for 

L  —  l 


anyxe{±_^}  , 


X  =  ais,  +  a2s2  +  . . .  +  aL-iSL_1  (111) 

where  a*  G  K  and  at  least  one  a,:  ^  0,  i  =  1, 2, . . . ,  L  —  1.  Using  (111), 

L—l 


E 

i=l 


m  // 

X  s„-  = 


L—l 


(ai  +  a2  +  . .  •  +  cll—  l)  • 


(112) 


Let  si  (1),  i  =  1,. . .  ,L  —  1  denote  the  first  element  of  st  and  x{\)  the  first 
element  of  x.  Without  loss  of  generality,  over  the  binary  antipodal  signature 
field  si  (1)  =  -^=,  i  =  1, 2, . . . ,  L  —  1.  Since,  x(l)  G  {±-^j=},  Y^h=i  ai  =  ±1  and 
the  result  follows.  § 


Theorem  1:  Let  Si  be  an  optimal  ( L  —  1,1/),  L  =  1  (mod  4),  initial  matrix. 


Create  the  (L,L)  set  S  =  | Si  [±^= 
Then, 


_TlT 


for  an  arbitrary  x  G 


L- 1 


argimn 


TSC(5  )  -  TSC(S7) 


=  1  +  2 


(£-2)(£-l) 

L2 


5No  systematic  construction  of  (L,L),  L  =  1  (mod  4),  TSC-optimal  sets  is  known  [2]- [4] . 

6 To  the  best  of  the  authors’  knowledge,  no  other  method  is  known  to  systematically  obtain 
an  (L  —  1,L)  TSC-optimal  set. 
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Proof:  Design  an  (L,  L— 1)  optimal  set  S  =  |<S  xj  where  x  G  {=b-^=}L  1  and 
S  =  s-j  s2...si_1  ,  si  €  {-^j=}L_1,  is  an  optimal  (orthogonal)  (L  —  1,  L  —  1) 
set.  From  [2], 

TSC(S'")  -  TSC(5")  =  {Ll21)2  +  2  £  (xTs")2  =  3(L~21)2.  (113) 

i=  1 


Hence, 


E  (XE')‘ 


L2 


r  i  lxL  m 

Now  append  any  row  vector  j  ±  ~^=  >  to  the  top  of  S 
using  (114), 


(114) 

Then, 


B.  Scalable  Conditionally  TSC-optimal  Binary  Sets 

Following  our  formulation  and  notation  in  (100),  (101),  conditionally  TSC- 
optimal  design  of  an  overloaded  signature  set  S  is  equivalent  to  designing  an 
extension  matrix  Se  such  that  TSC(<Se)  and  F(<Sj,<Se)  are  jointly  minimized. 
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In  the  following,  Hp  with  columns  Iq,  i  =  1, . . . ,  P,  represents  a  Hadamard  ma¬ 
trix7  of  size  P.  Below,  we  provide  scalable  conditionally  TSC-optimal  signature 
set  designs  for  all  lengths  L  and  overloading  of  (about)  100%. 


Case  1:  L  =  0  (mod  4),  1  <  N  <  L. 

Let  Si  =  ^HL  be  the  (L,L)  optimal  [2]  initial  matrix.  Obtain  (an  equivalent 
Hadamard  matrix)  Mp  by  multiplying  any  row  of  by  —1  and  select  any 

N  columns  from  Mp  to  form  the  extension  matrix  Se  =  [mi  m2 . . .  mjv]-  It 
can  be  shown  that 

TSC(<S)  =  L  +  3./V  (119) 

which  is  equal  to  the  bound  in  (104).  Hence,  S  is  conditionally  optimal. 

It  is  interesting  to  note  that  the  designs  of  Vanhaverbeke,  Moeneclaey,  and 
Sari  in  [15]  fall  under  this  case,  meet  the  bound,  and  are  conditionally  optimal. 


Case  2:  L  =  2  (mod  4),  1  <  N  <  L  +  2. 

Obtain  Hp+2  and  omit  the  first  two  rows  to  create  the  ( L  +  2,  L)  optimal  [3] 
initial  matrix  Si  =  -^=Hp+2.  Obtain  Mp+2  by  multiplying  any  of  the  last 
L  —  2  rows  of  ^Hp+2  by  —1.  Omit  the  first  two  rows  of  Mp+2,  select  any 
Lf  J  columns  that  begin  with  [-^  -^=]T  and  any  columns  that  begin  with 
[~j=  -^=]T  to  create  the  extension  set  Se  =  ML+2.  It  can  be  shown  that 


(  +N+  2(N_if_  +  2WCL+2),  TV  =  1  (mod  2) 

1'SOuS'i  -  ^ 

(  (L^f  +N+  2JV(w-2)  +  2JVCL+2) ,  ^  =  0  (mod  2) 

(120) 

which  is  equal  to  our  bounds  in  (107).  Hence,  S  is  conditionally  optimal. 

Case  3:  L  =  3  (mod  4),  1  <  N  <  L  +  1. 

Obtain  Hp+i  and  omit  the  first  row  to  create  the  (L  +  l,  L)  optimal  [2]  initial 
matrix  Si  =  -^HL+1.  Obtain  Mp+i  by  multiplying  any  of  the  last  L  —  1  rows 
of  ^Hp+i  by  —1.  Omit  the  first  row  of  Mp+i  and  select  any  N  columns  to 
create  the  extension  matrix  Se  =  Mi+1.  It  can  be  shown  that 


TSC(S) 


(L+  l)2 

L 


+  N  + 


N(N~1) 

L2 


2N(L  +  1) 
L 


(121) 


which  is  equal  to  our  bound  in  (110).  Hence,  S  is  conditionally  optimal. 

Case  4:  L  =  1  (mod  4),  1  <  N  <  L  —  2. 

Obtain  ^H^_i  and  multiply  the  /cth  column,  k  G  {1,  2, . . . ,  L  —  1},  by  — 1 

7 We  recall  that  a  Hadamard  matrix  Hjv  is  an  N  x  TV  matrix  with  elements  +1  or  —1  and 
orthogonal  columns,  H^Hjv  =  NIn  where  I^v  is  the  size- AT  identity  matrix.  A  necessary 
condition  for  a  Hadamard  matrix  to  exist  is  that  its  size  is  a  multiple  of  four,  except  for  the 
trivial  cases  of  size  one  or  two.  Without  loss  of  generality,  the  first  row  and  the  first  column 
of  Hat  contain  only  +1. 
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to  create  HL_1  =  h,  h2  . . .  hL_1 


.  Append  the  row  vector 


1  1 


Vl  Vl'"  Vl 


L- 1 


as  the  first  row  to  HL1  to  create  the  initial  matrix  Si  =  HL_1,  which  is  an 
( L  —  1  ,L)  optimal  set  [2].  Multiply  the  first  row  of  HL_1  by  —1  and  exclude 
the  kth  column  to  create  ML_2  =  rrq  m2  . . .  mL_2  .  Finally,  append  the  row 

vector  [— =  —=  . . .  — =1  as  the  first  row  to  M,  „  and  keep  any  N  columns  to 

1Vl  Vl  Vl  l~ 2 


L—2 


create  the  extension  matrix  Se  —  ML_2  = 


1N 


.  We  calculate 


TSC(S)  =  L—1+ —  ^[L  2)+N+N(N  1)+^(A-l)(A-2) 

Lz  Lz  Lz 

(122) 

which  is  equal  to  our  bound  in  (118).  Hence,  S  is  conditionally  optimal.  The 
derivation  of  (122)  is  provided  in  Appendix  A. 

To  conclude  the  presentation  and  establish  that  all  set  designs  S  provided 
herein  contain  unique  signatures,  we  need  only  to  observe  that  for  any  matrix 

Mp  =  [mi  m 2  . . . ,  mp]  ,  m,  g  j±-l=  j  ,  i  =  1, ...  ,P,  obtained  by  multiply¬ 
ing  any  row  of  HP  =  [hi  h2  . . . ,  hP] ,  h.j  €  {^Tp  }  ,  ^  =  1 ,  -  -  - ,  ^  by  1 , 

if  d  ( ht ,  nij )  denote  the  Hamming  distance  between  ht ,  ,  i,j  =  1,2,...,P. 

Then, 


|hj'm,-|  = 

1  2d(h?;,mJ) 

\  l  J  \ 

'  P 

By  construction  of  Mp,  if  *  =  j,  d  (h?;,  mj)  =  1  and  if  i  ^  j,  d  (hi;  m;/ )  =  -j  ±  1. 
From  (123) 

|hfmJ|  =  {  p’  i^jj  ■  (124) 

Hence,  every  column  of  Mp  is  different  from  every  column  of  Hp. 

Regarding  the  scalability  characteristics  of  the  presented  designs  for  medium- 
access-control  (MAC),  as  users  enter  the  system  we  may  assign  signatures  at  will 
first  from  Si  and  then  from  Se  without  having  to  change  the  existing/deployed 
signature  set.  In  the  event  of  a  user  leaving  the  system  from  the  extended 
set  Se  no  action  is  needed  to  maintain  conditional  TSC  optimality.  If  a  user 
leaves  from  the  initial  set  Si  and  Se  is  non-empty  then  its  signature  should 
be  reassigned  to  an  Se  user  (or  simply  to  the  next  incoming  user)  to  regain 
conditional  optimality. 


C.  Numerical  Evaluation  and  Comparisons 

In  this  section  we  evaluate  numerically  the  TSC  quality  of  the  proposed  over¬ 
loaded  signature  sets.  Fig.  14  shows  as  an  example  a  [19, 16]  conditionally 
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TSC-optimal  set  (obtained  under  Case  1).  The  natural  benchmark  is  the  un¬ 
conditionally  minimum  attainable  TSC  value  in  [2].  Fig.  15  plots  the  TSC  value 
of  our  Case  1,  L  =  16  design  against  the  bound  for  all  1  <  N  <  16.  Figs.  16,  17, 
and  18  carry  out  similar  comparative  studies  for  L  =  14  (Case  2),  L  =  15  (Case 
3),  and  L  =  17  (Case  4)  designs,  respectively. 

A  broader  comment  on  the  TSC  quality  of  the  proposed  conditionally  opti¬ 
mal  sets  for  varying  lengths  L  is  offered  by  the  following  proposition. 
Proposition  1:  The  TSC  of  the  designs  described  in  Section  III  is  always 
within  j  +  -i  from  the  bounds  in  [2] . 

Proof :  Let  S  be  a  Section  III  design  and  let  minTSC  denote  the  corresponding 
TSC  lower  bound  in  [2].  We  can  show  that  (see  Appendix  B) 


(  k 

4  ’ 


TSC (S)  -  minTSC  <  { 


l 

4  L 


1 

4  L 


2 

L2 


1 

4  L2 


1 

4 T2 


1 

2  » 


1 

4’ 


1 

4’ 


L  =  0(mod  4), 
with  equality  at  N 
L  =  2  (mod  4), 
with  equality  at  N 
L  =  3  (mod  4), 
with  equality  at  N 
L  =  l(mod  4). 


L 

2 


L+2 

2 


L+ 1 
2 


(125) 

Somewhat  finer  bounds  depending  on  the  form  of  the  total  number  of  signatures 
K  are  derived  in  Appendix  B  and  given  in  Table  III.  §j 


D.  Conclusions 

We  first  derived  lower  bounds  on  the  conditional  total  squared  correlation  of 
overloaded  binary  signature  sets  grown  on  full-load  TSC-optimal  sets  for  any 
signature  length  L  and  set  size  K  =  L  +  N+n,  1  <  N  <  L+2,  n  €  {—1, 0, 1, 2}. 
Then,  we  presented  conditionally  optimal  designs  for  all  such  (K,  L)  set  sizes 
(which  establish  the  tightness  of  the  corresponding  bounds) . 

We  analyzed  the  TSC  performance  of  the  conditionally  optimal  ( L  +  N,  L) 
signature  sets  against  the  TSC  of  optimal  (L  +  N,  L)  sets  and  found  that  we 
are  always  within  -j  +  1  or  better.  Direct  comparison  of  the  conditionally  TSC- 
optimal  designs  with  the  TSC  bounds  and  optimal  sets  in  [2]  shows  that:  Case 
1  when  N  =  {1, 2,L  —  2,L  —  1  ,L},  Case  2  when  N  =  {1,  2,  L,  L  +  1,  L  +  2}, 
Case  3  when  N  =  {1,  2,  L  —  1,  L,  L  +  1},  and  Case  4  when  N  =  {L  —  3,  L  — 
2}  are  unconditionally  TSC-optimal  (that  is  overloads  near  0%  or  100%  are 
handled  TSC  optimally).  Numerically,  maximum  TSC  performance  degradation 
is  observed  when  the  sets  are  about  50%  overloaded.  The  observed  degradation 
value  is,  arguably,  a  reasonable  price  to  pay  for  the  gained  scalability. 
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E.  Appendix  A:  Derivation  of  TSC  for  scalable  L  =  1  (mod  4) 
designs 


Consider  S  = 


1 5/  m 

TSC (5')  -  TSC(Sj)  =  1  +  2 

r 

=  1  +  2 

=  1  +  2 


,  i  =  1,2, . . .  ,N,  where  nij  = 


—  1  'T 

Vl  m 


iT 


.  Then, 


m  '  ,  h  ';  -  T 


EZy —  1 
7  —  1 


m  ,  h. 


EjE1  r?  -  r  EfE1  m  i  hj 


L2 


_  1  m  u' 

L  2-jj—l  111  i  “i 


(126) 

where  the  last  equality  is  from  (114).  From  Lemma  1,  Ef=iX  m  i  hj  =  ^^TT- 


We  now  show  that  for  our  construction  Ej^]1  m  i  hy  =  Et  for  all  i  =  l, N. 
Expanding  and  simplifying, 


L—l  T  _  L—l 

\  ’  'T  /  L,  —  3  \  -  'T  / 

2^mihj  =  —j—+  2s  m  '  E 


(127) 


j=i 


i=l  j+i 


(m  i  hj  =  -fodi  since  d(m  j,h  j)  =  1).  Also  note  that  V)  ^  i  m  hj  =  ±|- 
because  d(m  ;,h  j)  =  EE  ±  1.  Assume  Y^j=i  m  i  hj  =  f°r  which  a 

necessary  and  sufficient  condition  is  EjEi^+i  m  ;  hj  =  ~2^~2^  or  m  hj  = 
TT  Vj  yf  i.  By  design,  since  we  multiplied  the  kth  column  of  HL_1  by  — 1, 
when  j  =  fe,  m  j  hj  =  |;  hence,  EjE+^i  m  i  hj  E  z2Lkz21_  Therefore,  from 
Lemma  1,  J2j=i  m  ;  hj  =  EE  Substituting  (127)  in  (126)  we  have  TSC(5  )  — 
TSC  (Si)  =  1  +  2  (L~2h2L-1l  for  our  choice  of  extension  matrix  Se  which  is  the 
minimum  possible  by  Theorem  1.  Then, 

F(Si,Se)  =  ^(L-2)(L-1)  (128) 

and  since,  from  Table  I, 

TSC(SE  =  N  +  N(N  ~  X)  (129) 


we  obtain  (122). 

F.  Appendix  B:  Derivation  of  (125) 

Case  1:  L  =  0  (mod  4),  1  <  N  <  L. 

Consider  an  (L  +  N,  L)  conditionally  optimal  TSC  set  S.  Calculate  and  de¬ 
fine  the  following  TSC  (S)  —  minTSC  difference  functions  where  minTSC  is  the 
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unconditional  TSC  bound  of  an  (L  +  A,  L)  set: 


Ai  (N)  =  L  +  3N-  (Z/+/)2 ,  N  =  0  (mod  4) 

A2(iV)  =  L  +  3N-  (L+/)2  -  N  =  2  (mod  4)  (130) 

A 3  (AT)  =  L  +  3N  —  (L+JV)2  -  ^1,  AT  =1  (mod  2). 

To  find  the  integer  value  of  N  £  {1,  2, ...,  i}  that  maximizes  the  functions  above, 
we  first  extend  their  domain  to  all  real  N  £  [1,L]  and  denote  the  extended 
functions  by  A1(AT),  A2(N),  and  A 3(N),  respectively.  Differentiating  each  of 
these  functions  with  respect  to  N  over  the  continuous  field  yields  maximum  at 
N  =  for  all  three  cases.  Since  A2(|r)  <  A3(:|)  <  At(j)  =  Ai(-j),  we  set 
Ai(-|)  =  j  as  the  upper  bound.  This  bound  is  tight  only  for  N  =  0  (mod  4). 
For  N  =  2  (mod  4)  and  N  =  1  (mod  2),  A2(tt)  and  A3(-j)  form  tighter  upper 
bounds,  correspondingly.  These  bounds  are  presented  in  Table  III. 

Case  2:  L  =  2  (mod  4),  1  <  N  <  L  +  2. 

Consider  an  ( L  +  N,  L)  conditionally  optimal  TSC  set  S.  Calculate  and  define 
the  difference  functions  TSC  (S)  —  minTSC 

Ai(AT)  =  frac(L  +  2 )2L  +  N  +  2jv(^~2)  +  2Jv^+2>  _  (£+iv+2)2  ? 

N  =  0  (mod  4) 

A  on  A  (i+2)2  ,  ¥  ,  2N(N—2)  2N(L+2)  ( L+N+2 )2  2(L-2) 

A2(1V)  -  — E - hiVH - £5 - 1 - X - z - — > 

N  =  2  (mod  4) 

A3(AT)  £  +  + 

jV  =  1  (mod  2) . 

(131) 

Then,  extending  the  domain  of  the  three  functions  to  real  N  £  [1,L  +  2]  and 
maximizing  over  N,  we  find  that  the  maximum  of  all  domain-extended  functions 
A^Ar),  *  =  1,2,3,  is  at  N  =  Since  A^(^)  <  Ag(^)  <  A)(^)  = 
Ai(^j^),  we  set  =  ^_i__4r_|_Ias  the  upper  bound.  This  bound 

is  tight  only  for  N  =  0  (mod  4).  For  N  =  2  (mod  4)  and  N  =  1  (mod  2), 
^2(^2^)  ancl  ^3(^2^))  correspondingly,  form  tighter  upper  bounds  presented 
in  Table  III. 

Case  3:  L  =  3  (mod  4),  1  <  N  <  L  +  1. 

We  calculate  the  differences  TSC(S)  —  minTSC  as  follows 

(L+iv+1)2 

N  =  0  (mod  4) 

(L+IV+1)2  2(L— l)2 

L  L 2  ’ 

N  =  2  (mod  4) 
(L+iV-H)2  _  L  — 1 

Ar  =  1  (mod  2) . 

(132) 


Ai(AT)  =  +  N  +  +  2jV^+b  _ 

A2(7V)  =  ^±1)1  +  A  +  AMfll  +  um±ll  _ 

A3(7V)  =  +  A  +  +  2W(i±i)  _ 
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Then,  extending  the  domain  of  the  three  functions  to  N  G  [1,  L+ 1]  and  maximiz¬ 
ing  over  N,  we  find  that  the  maximum  of  all  domain-extended  functions  A  j(-ZV), 
i  =  1,2,3,  is  at  N  =  A±I.  Since  A2(A±I)  <  A3(A±1)  <  A^Ail)  =  Ai(A±1), 
we  set  Ai(A±i)  =  T  _  _L  ^  +  1  as  the  upper  bound.  This  bound  is  tight 

only  for  N  =  0  (mod  4).  For  N  =  2  (mod  4)  and  N  =  1  (mod  2),  A2(A±1)  and 
A3(A±1),  correspondingly,  form  tighter  upper  bounds  presented  in  Table  III. 

Case  4:  L  =  1  (mod  4),  1  <  N  <  L  —  2. 

As  in  the  previous  cases,  we  calculate  TSC(S)  —  minTSC  to  be 

A1(N)  A  L  -  1  +  +  N  +  +  2W^2-3L+2)  -  ^-1+JV)2 , 

iV  =  0  (mod  4) 

A2(A)  A  L  -  1  +  +  TV  +  ^11  +  mL'-3L+2)  _  (L-l+iV)^ 

-2(Ll~1)2,  Aee  2  (mod  4) 
A3(A)  ^  £  -  1  +  +  A  +  ^11  +  ^(L^-3L+2)  _  (L-1+JV)’ 

—  ±=±,  N=  1  (mod  2). 

(133) 

Letting  A  G  [1,£  —  2]  and  differentiating,  we  find  that  the  maximum  of  all 
domain-extended  functions  A t(N),  i  =  1,2,3,  is  at  A  =  -^A.  Note  that 
this  value  of  A  lies  only  in  the  integer  domain  of  A3(A).  Since  A2(AA!)  < 

A3(^)  <  we  set  =  J  ~  Jl  ~  ih  +  3  as  the  uPPer  bound. 

Unlike  previous  derivations,  this  bound  is  not  tight  for  A  =  0  (mod  4).  For 
A  =  2  (mod  4)  and  A  =  1  (mod  2),  A2(-^A)  and  A3(-^A),  correspondingly, 
form  tighter  upper  bounds  presented  in  Table  III. 
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TABLE  I  from  [2] 

Underloaded  Signature  Set  (k  <  l) 


Processing  Gain 

Number  of  Users 

Lower  Bound  on  TSC 

L  =  0  (mod  4) 

Any  K 

I< 

L  =  2  (mod  4) 

K  =  0  (mod  2) 

K  +  2k^~V 

K  =  1  (mod  2) 

L  =  1  (mod  2) 

Any  K 

V  |  K(K-l) 

T  j  2 

Overloaded  Signature  Set  (k  >  l) 


Number  of  Users 

Processing  Gain 

Lower  Bound  on  TSC 

K  =  0  (mod  4) 

Any  L 

K“ 

_ L _ 

K  =  2  (mod  4) 

L  =  0  (mod  2) 

L  =  1  (mod  2) 

K  =  1  (mod  2) 

Any  L 

,  L-l 

L  ^  L 
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TABLE  II 


> 

o 

ERLOADI] 

MG  OF  TSC-OPTIMAL  SETS 

Signature 

Length 

BaseTSC- 

optimalSet 

Add  —  on 
Signatures 

Lower  Bound  on  TSC 

L  =  0  (mod  4) 

(L,L) 

1  <  N  <  L 

L  +  3N 

L  =  2  (mod  4) 

(L  +  2,  L) 

1  <  N < L  +  2 

(L+  W  ,  2(N—1)Z  ,  2N(L+2 ) 

L  +  N  +  - 12 -  +  L 

N  =  1  (mod  2); 

(i+2)2  ,  Ar  |  2N(N-1)  2N(L+2 ) 

L  1  iV  1  L2  1  L  5 

A  =  0  (mod  2) 

L  =  3  (mod  4) 

(L+l.L) 

1<N<L+ 1 

(i+1)2  ,  Ar  ;  N(N—1)  ,  2JV(L+1) 

r,  +  H  T2 - 1-  r. 

L  =  1  (mod  4) 

(L-l.L) 
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TABLE  III 

UPPER  BOUNDS  ON  TSC(S)  -  minTSC 

(tight  for  ( K  =  0  (mod  4),  L  =  0,  2,  3  (mod  4))  and  ( K  =  1  (mod  2),  L  =  1  (mod  4))) 


I\  =  0  (mod  4) 

K  =  2  (mod  4) 

K  =  1  (mod  2) 

L  =  0  (mod  4) 

L 

4 

L  9  L— 2 

_ 4 _ _ 

L  =  2  (mod  4) 

L  1  2,1 

4  L  L2  '  2 

L  =  3  (mod  4) 

L  1  1  ,  1 

4  4L  4L2  '  4 

Kl 

lisa 

L  =  1  (mod  4) 
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Fig.  14:  (19, 16)  conditionally  TSC-optimal  signature  set. 
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VI.  New  Bounds  on  the  Aperiodic  Total  Squared 
Correlation  of  Binary  Signature  Sets  and 
Optimal  Designs 

This  work  has  been  presented  at  IEEE  ICC  2007,  Scotland,  UK. 

In  code-division  multiple  access  (CDMA)  systems,  each  of  the  K  participat¬ 
ing  users  is  assigned  a  distinct  spreading  sequence  or  signature  s k  €  CL,  ||sfc||  = 
1,  k  =  1,2,...  AT.  Organizing  the  K  signatures  as  the  columns  of  a  matrix,  we 
define  the  signature  matrix 


S=  [Sl  s2...sA-]  gCLxK.  (134) 

In  synchronous  CDMA  communications  over  ideal  Nyquist  channels,  we  are 
interested  in  using  signature  sets  with  minimum  total  squared  correlation  (TSC) 

K  K 

TSC  (5)  =  EEKsif  (135) 

i=l  j= 1 

where  H  denotes  the  Hermitian  operator.  For  overloaded  (AT  >  L)  complex/real¬ 
valued  signature  sets  S  £  CLxK  or  MLxK,  TSC(<S)  >  [1]  (of  course, 

TSC  (S)  >  K  if  K  <  L).  Overloaded  ( K  >  L)  sets  with  TSC  equal  to 
have  been  known  as  Welch-bound-equality  (WBE)  sets.  A  few  algorithms  and 
studies  for  the  design  of  complex  or  real  WBE  signature  sets  can  be  found  in 
recent  literature. 

In  digital  transmission  systems,  however,  it  is  necessary  to  have  finite-alphabet 
signature  sets.  Recently,  new  bounds  were  derived  on  the  TSC  of  binary  antipo¬ 
dal  signature  sets  together  with  optimal  designs  for  arbitrary  signature  lengths 
and  set  sizes  [2]- [4].  The  sum  capacity,  total  asymptotic  efficiency,  and  maxi¬ 
mum  squared  correlation  of  minimum-TSC  optimal  binary  sets  were  evaluated 
in  [5] .  The  sum  capacity  of  several  other  signature  set  designs  under  potentially 
a  binary  or  quaternary  alphabet  was  examined  in  [6].  In  this  paper,  we  deal 
strictly  with  binary  antipodal  signatures  represented  as  (unnormalized)  column 
vectors  in  {±1}  . 

When  asynchronous  code-division  multiplexing  is  attempted  and/or  the  chan¬ 
nel  exhibits  multipath  behavior,  apart  from  the  total  squared  correlation  be¬ 
tween  signatures  we  are  also  concerned  about  the  individual  periodic  and  ape¬ 
riodic  correlation  values  between  the  signatures.  For  notational  simplicity,  we 
define  a^;  <E  {0,  ±l}lx(2L~1}  by 

aqi  =  [0_M)  sfc(l)  Sfc(2)...sfc(L)  0_Mf],  l  =  0, 1, . . . ,  2L  -  2,  (136) 

l  L-l-l 

to  represent  the  zero-padded-by-(L-l)  cycli-shifted-by-l  version  of  £  {±l}lxi, 
k  =  1,2,...,  A"  (T  denotes  the  transpose  operator  and  s  k(i),i  =  1,2,...,  A, 
k  —  1,2,...,  A",  refers  to  the  ith  element  of  vector  s*,).  Next,  we  define  the 
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zero-padded  cyclic  extension  matrix  Szpc  £  {0,  ±1}<'2L  1')xK(2L  b  0f  the  signa¬ 
ture  set  S  £  {±l}ixA  as 


A 


Szpc  ~  [al|0  a2 1 0  ■  •  -aA'|0  al|l  a2|l  •  •  -aX|l 
al|2L— 2  a2|2L— 2  •  •  •  aK|2A-2] 


(137) 


and  the  aperiodic  total  squared  correlation  (ATSC)8  of  the  signature  set  S  as 
the  total  squared  correlation  (TSC)  of  the  matrix  Szpc 


ATSC(S)  =  TSC(5zpc).  (138) 

Then, 

K  2L-2  K  K  2L-2 

ATSC(S)  =  E  E  \al\oak\l\  +(2L-1)E  E  E  \aT\oaJ\l\  ■  (139) 

k=  1  1=0  i=  1  j=l,i^j  1=0 


The  first  term  in  (139)  contains  all  aperiodic  auto-correlation  contribu¬ 
tions;  the  second  term  contains  all  aperiodic  cross-correlation  contributions.  For 
complex/real- valued  signature  sets  5  £  CLxK  or  RLxK,  ATSC(S)  >  K2L2(2L— 
1)  directly  by  the  TSC  bound  in  [1]  and  our  formulation  (4), (5).  Hence,  these 
constructions  establish  tightness  of  the  above  ATSC  bound  over  the  binary  do¬ 
main  for  those  specific  cases. 

In  our  work,  we  first  derived  new  lower  bounds  on  the  ATSC  of  binary 
signature  sets  for  all  possible  values  of  K  (number  of  signatures)  and  L  (sig¬ 
nature  length).  Then,  we  presented  optimal  designs  for  several  ( K,L )  pairs. 
The  designs  are  based  on  Hadamard  matrix  transformations  and  serve  as  proof- 
by-construction  for  the  tightness  of  the  corresponding  ATSC  bounds.  The  new 
bounds  and  designs  have  already  been  publiced  and  presented  in  the  2007  IEEE 
International  Conference  on  Communications  (ICC). 


A.  Discussion  and  Examples 

The  optimal  design  cases  that  we  produced  constitute  proof-by-construction  of 
the  tightness  of  our  corresponding  ATSC  bounds.  Under  the  KL  =  0  (mod  4) 
design  case,  we  relied  significantly  on  ACS  literature;  new  advancements  have 
been  made  under  the  KL  =  2  (mod  4)  and  KL  =  1  (mod  2)  design  cases.  To 
acquire  a  quantitative  feeling  of  the  coverage  of  the  presented  designs,  if  we 
restrict  the  domain  of  K,L  to  {1,2,...,  256}  (at  present,  it  does  not  appear  of 
practical  interest  to  consider  code-division  applications  outside  this  parameter 
range),  we  can  calculate  that  Underloaded  Cases  1  through  6  and  Overloaded 
Cases  1  through  4  together  represent  27.91%  of  all  possible  combination  pairs 
[K,  L)  £  {1,  2, . . . ,  256}“.  Maximum  coverage  in  the  number  of  users  is  available 

8  Another  signature  set  correlation  metric  of  interest  is  the  periodic  total  squared  correlation 
(PTSC)  which  is  lower  bounded  by  K2L 3  for  complex/real- valued  signature  sets  S  S  C L  x  K 
or  |P' x  K  [1] ,  [10] ,  [9] .  Recently,  new  lower  bounds  on  the  PTSC  of  binary  antipodal  signature 
sets  were  derived,  together  with  optimal  designs  for  many  signature  lengths  and  set  sizes,  in 
[9],  [10]. 
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in  overloaded  systems  when  L  G  he.  optimal  signature  sets  are  available 

for  all  values  of  K  >  L. 

Certainly,  tightness  of  the  bounds  and  optimal  ATSC  designs  under  the 
remaining  cases  is  an  important  open  research  problem.  Two  results  need  to  be 
highlighted  in  this  regard;  the  non-existence  of  Barker  sequences  for  60  <  L  < 
1022  [14]  and  the  non-existence  of  (2  ,L)  ACS  sets  if  all  prime  divisors  of  L  are 
not  congruent  to  1  (mod  4). 

Since  the  ATSC  of  each  signature  set  that  we  constructed  equals  the  corre¬ 
sponding  lower  bound,  by  our  own  Theorem  Number  2  all  designs  are  jointly 
ATSC-  and  PTSC-optimal.9  In  addition,  direct  comparison  of  our  ATSC- 
optimal  designs  with  the  TSC  bounds  and  optimal  sets  in  [2]  shows  that  Un¬ 
derloaded  Cases  1  and  2  when  K  is  a  power  of  2,  Underloaded  Case  3(i),  and 
Overloaded  Cases  1  through  4  are  triple,  ATSC-,  PTSC-  and  TSC-optimal. 
This  remarkable  conclusion  captures  and  demonstrates  the  inherent  robustness 
of  these  binary  designs. 

We  will  now  establish  that  the  familiar  Gold  and  “small”  or  “large”  Kasami 
sets  [15],  [16],  which  have  been  widely  used  for  their  correlation  properties  [8], 
are  not  ATSC  optimal  in  general.  It  suffices  to  find  at  least  one  ( K ,  L)  signature 
set  in  the  Gold/Kasami-compatible  range  that  has  lower  (preferably  minimum) 
ATSC  than  the  corresponding  Gold/Kasami  set  in  order  to  prove  sub-optimality 
of  the  latter.  Fig.  19(a)  provides  an  example  of  a  Gold  set  £31x24  with  ATSC  = 
43474212;  Fig.  19(b)  provides  an  optimal  signature  set  S]}}*  24  designed  under 
our  Underloaded  Case  3(ii)  with  minimum  ATSC  =  33765696.  In  Fig.  20(a) 
we  show  a  (2, 15)  signature  set  with  lower  (but  not  minimum)  ATSC  =  33292 
than  a  (2,15)  small-set  Kasami  design,  shown  in  Fig.  20(b),  with  ATSC  = 
34220.  We  follow  with  Fig.  21  which  provides  a  large-set  Kasami  design  /C[s5x8 
with  ATSC  =  441264  and  an  optimal  signature  set  8  designed  under  our 
Underloaded  Case  3(ii)  with  minimum  ATSC  =  417600. 

We  conclude  with  an  example  of  an  overloaded  ATSC-optimal  design  S°lx25 
given  in  Fig.  22.  The  set  is  designed  under  our  Overloaded  Case  2  and  has 
minimum  ATSC  value  ATSC  =  3308526. 

B.  Conclusions 

We  derived  lower  bounds  on  the  aperiodic  total  squared  correlation  of  binary 
antipodal  signature  sets  for  any  set  size  ( K ,  L) .  We  provided  optimal  designs  for 
a  range  of  (A',  L)  pairs  that  establish  the  tightness  of  the  corresponding  lower 
bounds.  The  constructions  include  underloaded  ( K  <  L)  and  overloaded  ( K  > 
L)  design  cases  and  cover,  as  an  example,  27.91%  of  all  possible  combinations 
of  K,  L  in  {1,2,...,  256}. 

Side  results  of  this  work  include  establishing  that  maximal-merit-factor  se¬ 
quences  (and,  hence,  Barker  sequences)  are  ATSC-optimal  and  that  neither 
Gold  nor  small-  or  large-set  Kasami  sequences  are  ATSC-optimal  in  general. 

''The  Barker  sequence  g4  applicable  to  Underloaded  Case  5  and  Overloaded  Cases  3  and 
4  is  double,  ATSC-  and  PTSC-optimal. 
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In  view  of  these  findings,  the  developed  ATSC-optimal  sets  take  precedence  in 
code-division  multiplexing  applications. 
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+  -  +  --  -  +  +  +  -  +  -  +  -  +  --  -  +  +  +  -  +  - 

+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 

-  +  +  -  +  +  +  --  -  +  --  +  +  -  +  +  +  --  -  +  - 

601  y94  —  -  +  -  +  +  -  +  +  +  —  -  +  -  +  +  -  +  +  +  —  —  — 

+  +  -  +  -  +  +  +  +  -  +  + 

+  +  +  -  +  -  +  +  +  +  -  +  + 
+  +  +  +  -  +  --  +  --  -  +  +  +  +  -  +  --  +  --  - 

-  +  --  -  +  -  +  +  -  +  +  -  +  --  -  +  -  +  +  -  +  + 

-  +  +  --  -  +  -  +  +  -  +  -  +  +  --  -  +  -  +  +  -  + 

-  +  +  +  --  -  +  -  +  +  --  +  +  +  --  -  +  -  +  +  - 

+  +  --  -  +  +  +  -  +  --  +  +  --  -  +  +  +  -  +  -- 

-  +  -  +  +  +  --  -  +  -  +  -  +  -  +  +  +  --  -  +  -  + 

+  +  +  +  +  +  +  +  +  +  +  +  --  --  --  --  --  -- 

+  --  +  --  -  +  +  +  -  +  -  +  +  -  +  +  +  --  -  +  - 

+  +  -  +  +  +  -  +  +  -  +  +  +  + 

+  -  +  --  +  --  -  +  +  +  -  +  -  +  +  -  +  +  +  --  - 

+  +  -  +  --  +  --  -  +  +  --  +  -  +  +  -  +  +  +  -- 

+  +  +  -  +  --  +  --  -  +  --  -  +  -  +  +  -  +  +  +  - 

+  +  +  +  -  +  -  +  -  +  +  -  +  +  + 


(b) 


Fig.  19.  (a)  03ix24  Gold  set  with  ATSC  =  43474212.  (b)  Optimal  signature  set 
‘-*31x24  designed  under  Underloaded  Case  3(ii)  with  ATSC  =  (24)2(31)2(2(31)  —  1)  = 
33765696. 
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(a)  (b) 

Fig.  20.  (a)  1C  i5X2  small-set  Kasami  with  AT  SC  =  34220.  (b)  Numerically  generated 
<Si5x2  signature  set  with  AT  SC  =  33292. 


+  +  -  +  --  +  - 

-  +  --  +  -  +  + 

+  +  -  +  -  +  -  + 

--  +  +  --  +  + 

---  +  +  +  +  - 

--  +  -  +  +  -  + 

+  -  +  -  +  -  +  - 

+  -  +  +  -  +  -- 

---  +  +  +  +  - 

+  -  +  -  +  - 

+  -  +  +  --  +  + 


’+  +  +  +  +  +  +  +“ 

+  -  +  -  +  -  +  - 

+  +  --  +  +  -- 
+  --  +  +  --  + 
+  +  +  +  +  +  +  + 

-  +  -  +  -  +  -  + 

+  +  --  +  +  -- 
-  +  +  --  +  +  - 
+  +  +  +  +  +  +  + 

+  -  +  -  +  -  +  - 

--  +  +  --  +  + 
-  +  +  --  +  +  - 
+  +  +  +  --  -- 
+  -  +  --  +  -  + 
+  +  --  --  +  + 


(a)  (b) 

Fig.  21.  (a)  JC  15x8  large-set  Kasami  with  AT  SC  =  441264.  (b)  Optimal  < S^g 

designed  under  Underloaded  Case  3(ii)  with  AT  SC  =  (8)2(15)2(2(15)  —  1)  =  417600. 


'+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +’ 

+  --  +  --  -  +  +  +  -  +  +  --  +  --  -  +  +  +  -  +  + 

+  -  +  --  +  --  -  +  +  +  +  -  +  --  +  --  -  +  +  +  + 

+  +  -  +  -  +  +  +  +  -  +  - 

+  +  +  -  +  --  +  --  -  +  +  +  +  -  +  --  +  --  -  +  + 

oopt  _  +  +  +  +  -  +  -  +  +  +  +  -  + 

°14x25  —  +  -  +  +  +  -  +  +  -  +  +  +  -  +  - 

+  --  +  +  +  -  +  --  +  -  +  --  +  +  +  -  +  --  +  -  + 

+  --  -  +  +  +  -  +  --  +  +  --  -  +  +  +  -  +  --  +  + 

+  +  -  +  +  +  -  +  -  +  +  +  -  +  - 
+  -  +  -  +  +  +  -  +  -  +  -  +  -  +  +  +  -  + 

+  +  +  +  +  +  +  +  +  +  +  +  - 
+  -  +  +  +  -  +  -  +  +  -  +  +  + 


Fig.  22.  Optimal  signature  set  S° 4^25  designed  under  Overloaded  Case  2  with  mini¬ 
mum  ATSC  =  (25)2(14)2(2(14)  -  1)  +  (2(14)  -  1)(3(14)  -  4)  =  3308526. 
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VII.  Upward  Scaling  of  Minimum-TSC  Binary 
Signature  Sets 

The  work  has  been  published  in  the  IEEE  Communications  Letters,  Nov.  2007. 

In  multiuser  communication  systems  that  follow  the  code-division  multi¬ 
plexing  paradigm,  multiple  signals  are  transmitted  simultaneously  in  time  and 
frequency.  Each  signal  -potentially  associated  with  a  distinct  user-  is  assigned 
an  individual  signature  (spreading  code) .  A  fundamental  measure  of  the  quality 
of  the  code-division  communication  link  is  the  total  squared  correlation  (TSC) 
[1]  over  the  set  of  assigned  signatures.  For  a  TT-signal  system  with  signature 
length  L,  if  the  signature  set  is  denoted  by  S  =  {si,s2,---  , s k},  ||sj||  =  L, 
i  =  1,  2,  •  •  •  ,  K,  then  the  TSC  of  the  signature  set  S  is  defined  as  the  sum  of 
the  squared  magnitudes  of  all  inner  products  between  signatures, 

TSC{S)  =  £X>M2-  (no) 

i=  1  J= 1 

In  the  theoretical  context  of  TSC  optimized  signature  sets  that  are  real  (or 
complex)  valued,  one  may  consider  the  early  work  of  Welch  [1]  followed  by 
representative  works  in  [2]- [13]. 

Findings  in  [1]- [13]  constitute  only  pertinent  performance  bounds  for  digital 
communication  systems  with  digital  signatures.  Recently,  new  bounds  on  the 
TSC  of  binary  signature  sets  were  presented  [14]  that  led  to  minimum-TSC 
optimal  binary  signature  set  designs  for  almost  all10  signature  lengths  and  set 
sizes  [14]- [14].  The  sum  capacity,  total  asymptotic  efficiency,  and  maximum 
squared  correlation  of  the  minimum-TSC  binary  sets  were  evaluated  in  [15].  The 
user  capacity  of  minimum  and  non-minimum-TSC  binary  sets  was  identified  and 
compared  in  [18].  A  procedure  to  find  minimum-TSC  binary  signature  sets  with 
low  cross-correlation  spectrum  was  presented  in  [19].  A  binary  code  allocation 
scheme  is  examined  in  [20]  and  tested  against  the  Karystinos-Pados  (KP)  bounds 
and  sequence  sets  [14]. 

The  technical  problem  that  we  consider  in  this  letter  is  upward  scaling  of 
an  overloaded  ( K  >  L)  min-TSC  binary  set.  Consider  the  min-TSC  optimal 
constructions  in  [14]- [14]  for  any  K  =  L  (fully  loaded  systems).  Subsets  of 
K  <  L  signatures  maintain  TSC  optimality  and  signatures  can  be  returned 
and  reassigned  without  loss  of  optimality.  This  is  not  the  case  unfortunately, 
in  general,  given  a  min-TSC  overloaded  set  ( K ,  L)  where  K  >  L.  Addition  of 
a  signature,  for  example,  may  require  complete  redesign/reassignment  of  the 
(K  +  1,  L)  set.  In  this  letter,  motivated  by  our  prior  work  in  [21],  we  develop  a 
novel  scheme  that  returns  a  binary  signature  for  the  new  signal  in  the  system 
that  lies  near  the  continuous-valued  arcs  of  least  TSC  increase.  The  quality  of 
the  new  binary  signature  design  is  tested  directly  against  the  KP  TSC  bound. 

1  °The  case  K  =  L  =  1  (mod  4)  remains  the  only  open  problem  at  present  in  min-TSC 
optimal  binary  set  design  [14]- [14], 
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The  rest  of  this  paper  is  organized  as  follows.  Section  II  presents  our  formu¬ 
lation  of  the  optimization  problem  under  consideration.  The  proposed  algorithm 
is  presented  in  Section  III.  Studies  and  comparisons  with  the  TSC  bound  are 
included  in  Section  IV. 

A.  Formulation 

We  consider  a  code  division  multiplexing  system  with  code  length  L  and  K  >  L 
signals  (overloaded).  The  K  signals  utilize  a  minimum  TSC  optimal  binary  sig¬ 
nature  set  S  designed  according  to  [14]-[14],  S  =  {si,s2,---  ,s k},  s*  €  {±1}L, 
i  =  l,---  ,  K.  When  a  new  signal  enters  this  system  with  signature  s/f+i  G 
{±1}L,  the  TSC  of  the  K  +  1  signatures  given  the  signatures  of  the  K  preex¬ 
isting  signals  is 

K+l  K+ 1 

TSCk+ i,k  =  E  E  ls^l2 

i=  1  3  =  1 

-  EE!  sfsjf  +  |sx+lS/f+l  2  +  2E  lSElsi|2 

i=  1  j= 1  i=  1 

K 

=  TSCk  L2  +  ^  ^  SjsfsK+i  (141) 

i=  1 

where  TSCk  denotes  the  TSC  of  the  K  preexisting  signals  in  the  system  that 
utilize  a  minimum  TSC  binary  signature  set.  If  we  denote  the  autocorrelation 
matrix  of  the  preexisting  K  signatures  by 

K 

RA-  =  Es*sf>  (142) 

i=  1 

(141)  shows  that  conditional  minimization  of  TSCk+i\k  with  respect  to  s^+i 
for  fixed  (min-TSC- valued)  TSCk  reduces  to 

sk+i  =  arg  min  stRa'S.  (143) 

Exhaustive  search  over  all  2L  vectors  in  {±1}L  to  find  the  one  that  minimizes 
(143)  is,  of  course,  unacceptable  computationally  even  for  moderate  values  of  L. 
Below,  we  propose  a  low  cost  search  algorithm  that  creates  a  signature  candidate 
list  of  size  linear  in  L. 

B.  Proposed  Algorithm 

Let  /( s)  denote  the  cost  function  in  (143)  that  we  try  to  minimize, 

/( s)  =  stHks. 
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(144) 


We  relax,  for  a  moment,  our  requirement  for  binary  antipodal  signature  (se¬ 
quence)  alphabet  and  assume,  instead,  that  s  is  real-valued  with  the  same  norm, 
s  £  sts  =  L.  Then,  the  real-field  optimization  problem  becomes 

siK+i,oPt  =  ar9  min  /( s)  (145) 

sGKl ,Sj  S =L 

where  the  superscript  (r)  indicates  that  t  is  real-valued.  The  optimiza¬ 

tion  in  (247)  is  carried  over  a  hypersphere  in  RL  of  radius  L  centered  at  the 
origin. 

Let  {qi,  q2,  •  •  •  ,  qi,}  be  the  L  eigenvectors  of  with  corresponding  eigen¬ 
values  Ai  <  A2  <  •  •  •  <  \l-  The  real- valued  sequence  that  minimizes  the  right- 
hand-side  of  (247)  is  well  known  and  equal  to  the  eigenvector  that  corresponds 
to  the  minimum  eigenvalue  of  the  matrix  R/<, 

sK)+i,oPt  =  ar9  min  /(s)  =  qi.  (146) 

’  y  sGlL,sHs =L 

Consider  now  the  following  L  —  1  lines  in  RL, 

qi-pqi7  i  =  2, • • •  ,L,  p£R.  (147) 

The  lines  in  (249)  lie  on  the  plane  that  is  tangent  to  the  searching  hypersphere, 
pass  through  the  real  minimizer  qi,  and  define  mutually  orthogonal  directions 
of  least  increase  in  the  cost  function  /( s),  s  £  RL,  from  the  optimum  point 

(r) 

sK+ 1  t  =  qi.  That  is,  algebraically,  for  fixed  p  and  i  =  2,3,  •  •  •  ,  L, 

q*  =  arg  min  {/( qi)  -  /(  Ql  ^ =)}  (148) 

venv  y'l  +  p2 

where  fiv  =  {v  :  vTqj  =  0,  j  =  1, 2,  •  •  •  ,  i  —  1,  vTv  =  1}. 

Projection  of  the  above  least-increase  lines  onto  the  searching  hypersphere 
results  in  the  slowest  increase  arcs  given  by 


i  =  2,-  •  •  ,L,  p  £] R. 


(149) 


The  slowest  increase  arcs  in  (250)  trace  the  searching  hypersphere,  extend  from 

— qi  to  qi,  i  =  2,  •  •  •  ,  L,  and  pass  through  the  optimum  point  opt.  On  each 

of  these  arcs,  the  function  /( s),  s  £  KL,  takes  values  in  [AiL,  AiL],  *  =  2,  •  •  •  ,L, 

respectively.  Our  objective  is  to  identify  the  binary  sequences  that  are  closest  in 

the  /2-sense  to  the  above  least-increase  arcs  [22].  It  is  important  to  note  that  for 

any  given  i  £  {2, 3,  •  •  •  ,  L},  the  binary  sequences  that  are  closest  in  the  /2-sense 

to  the  arc  qiTpq' ,  pgK,  can  be  expressed  as  sg?i(qi  —  pqd;  the  set  of  all  binary 
-sA+p2 

signatures  of  the  form  s(/n(qi  —  pqi),  p  £  R,  has  cardinality  at  most  L  +  1. 

Our  proposed  binary  signature  search  algorithm  is  described  below. 

Algorithm 

For  /  =  2,  •  •  •  ,P(P<L)  do: 
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Step  1 
Calculate 

Pn  —  Ql,n  /  Qi,m  ^  —  1,  *  ,  (150) 


where  <p>n,  n  =  l,-‘‘  , -L,  is  the  nth  element  of  the  vector  q^.  The  points  in 
(150)  define  non-overlapping  intervals  in  R.  such  that,  for  any  given  i,  signatures 
of  the  form  sgn{ qi  —  pqd  that  correspond  to  values  of  p  in  adjacent  intervals 
have  opposite  signs  in  exactly  one  coordinate.  Let  p[  <  p'2  <  •  •  •  <  p'L  be  a 
rearrangement  of  pi,  p2,  •  •  •  ,  Pl  in  ascending  order  and  let  u  be  the  index  of  the 
first  positive  element  in  the  ordered  sequence  p\ ,  p'2 ,  •  •  •  ,  p'L,  i.e.  p'u_1  <  0  <  p'u, 
u  G  {1,  2,  ■  ■  ■  ,Lj. 

Step  2 

Find  the  binary  sequence  that  is  closest  to  the  arc  for  each  interval  of 

p,  i.e.  for  p  G  (-oo,  p[),  (pi,  p2),  ■  ■  ■  ,  (p'L_i,  p'L),  (p'L,  oo).  There  are  L  +  1  bi¬ 
nary  sequences  in  total  denoted  as  s^*+1,  s^+1,  •  •  •  ,  s^z+1  that  can  be  computed 
recursively, 


S<M  _ 

sir+i 


=  sgn[  qi], 


i+l,t  _  l,i  _  r,  0,t  _  /  _  n  i  . 

»X+1  —  sjf+l  ^bK+l,u+ltlu+1’  1  —  ui  ri 


l,i  _  Z+l,i  r\  0 ,i 

sir+i  —  s/c+i  —  zsjr+i,u+ie“+0 


l  =  -1,-2, 


,L  —  u, 

■  ,  1  -  u, 


(151) 

(152) 

(153) 


where  s^+i  is  the  (u  +  element  of  s^*+1  and  eu+i  is  the  (u  +  /)-unit 
vector  in  RL. 

Step  3 

Evaluate  TSCk+ i  for  each  binary  signature  s^-+1  returned  by  Step  2,  l  = 
0,  •  •  •  ,  L,  and  choose  the  binary  signature  that  gives  minimum  TSC.  I 

The  above  procedure  identifies  and  composes  L  +  l  binary  sequences  for  each 
one  of  a  total  P  —  1  (P  <  L)  slowest  increase  arcs11;  we  select  the  signature 
that  gives  minimum  TSC  among  all  (P  —  1)L+  1  candidates.  Our  experimental 
studies  indicate  that  P  =  2  or  3  (i.e.  one  or  two  slowest  increase  arcs)  is  sufficient 
to  closely  approximate  the  performance  level  reached  when  all  possible  slowest 
increase  arcs  are  considered  (P  =  L). 


C.  Experimental  Studies 

We  consider  a  code-division  multiplexing  system  with  signature  length  L  = 
16.  We  assume  that  initially  there  are  K  =  16  users  (fully-loaded  system) 
that  utilize  a  minimum  TSC  optimal  binary  signature  set  (orthogonal  Walsh- 
Hadamard  in  this  trivial  case). 

In  Fig.  29,  we  plot  the  TSC  of  the  code-division  multiplexing  system  ( K  + 
1  ,L=  16)  for  K  =  16  up  to  31  where  each  ( K,L )  signature  set  is  optimally 
min-TSC  designed  by  [14]  and  the  K  +  1  signature  is  added  by  the  procedure 
in  the  previous  section.  The  quality  of  the  resulting  ( I\  +  1  ,L)  design  is  tested 
against  the  (K  +  1,  L)  TSC  bound  of  [14]. 

11The  binary  sequence  j  =  sgn[qi]  is  common  in  all  slowest  increase  arcs. 
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Fig.  30  repeats  the  study  for  the  (K  +  1,L  =  31)  system  set-up  with  K  = 
31  up  to  61.  As  with  the  studies  of  Fig.  1,  the  comparison  with  theoretical 
minimum  TSC  bounds  is  very  favorable.  Frequently,  the  resulting  sequence  set 
is  absolutely  TSC  optimal. 


Fig.  25:  TSC  of  ( K  +  1  ,L  =  16)  signature  set. 
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VIII.  Short-data-record  Adaptive  Detection 

The  work  has  been  presented  at  the  2007  IEEE  Radar  Conference,  Boston,  MA. 

In  the  context  of  an  array  radar  application  with  M  antenna  elements  (spa¬ 
tial  channels)  and  N  pulses  per  coherent  processing  interval  (CPI),  optimum 
signal  detection  in  the  presence  of  spatial  and  temporal  Gaussian  noise  requires 
joint  space-time  matched  filtering  in  the  MN  complex  vector  space.  Since 
prior  knowledge  of  the  noise  covariance  matrix  is  not  available,  the  maximum- 
likelihood  sample  average  estimate  is  often  used,  developed  from  K  noise  only 
vector  samples  that  correspond  to  distinct  range  cells.  This  is  the  so  called 
secondary  data  set.  Then,  the  inverse  of  the  sample-matrix  is  considered  as 
an  estimate  of  the  inverse  covariance  matrix.  This  approach  is  known  as  the 
Sample-Matrix-Inversion  method  (SMI)  and  in  [1]  it  was  shown  to  outper¬ 
form  the  recursive  least-mean-squares  (LMS)  adaptive  implementation  of  the 
matched  filter  (MF)  processor  in  terms  of  convergence  rate  and  small-sample 
output  SNR  characteristics.  Still,  it  was  found  that  K  >  2 MN  independent 
and  identically  distributed  (i.i.d.)  training  data  samples  are  needed  to  maintain 
with  probability  1/2  a  loss  lower  than  or  equal  to  3dB  compared  to  the  ideal 
MF.  Simple  variance  normalization  of  the  MF  decision  statistic  led  to  a  Con¬ 
stant  False  Alarm  Rate  (CFAR)  test  [2]  for  Neyman-Pearson  detection.  System 
optimization  in  the  generalized  likelihood  ratio  (GLR)  sense  was  pursued  in  [3] 
and  analyzed  in  [4].  The  resulting  test  statistic  offers  embedded  CFAR  behav¬ 
ior,  converges  to  the  ideal  MF  solution  at  least  in  probability  as  the  number  of 
secondary  data  K  grows,  and  it  is  shown  in  [2]  to  outperform  the  SMI  approach 
for  K  =  2 MN,  except  for  high  SNR  regions. 

It  is  important  to  note  that  the  SMI  and  GLR  tests  are  both  asymptotically 
optimal  to  the  extent  that  they  converge  in  a  probabilistic  sense  to  the  optimum 
ideal  MF  as  the  size  of  the  secondary  data  set  grows  to  infinity.  However,  for 
finite  sample  support  no  optimality  can  be  claimed  in  either  case  and  superior 
probability  of  detection  performance  for  a  fixed  false  alarm  rate  is  theoretically 
possible  by  other  filtering  means.  Moreover,  both  methods  share  the  need  to 
invert  the  sample  covariance  matrix  of  the  noise  process. 

Arguably,  in  airborne  surveillance  systems  the  training  data  size  require¬ 
ments  make  the  practicality  of  these  approaches  questionable  even  for  mod¬ 
erate  values  of  M  and  N.  This  is  particularly  true  if  we  consider  a  highly 
non-stationary,  non-homogeneous  [5]  operating  environment,  typically  encoun¬ 
tered  in  practice,  that  necessitates  brief  data  collection.  In  this  present  work, 
we  consider  the  auxiliary- vector  (AV)  iterative  algorithm  [7],  [8]  to  develop  an 
infinite  sequence  of  decision  statistics.  The  objective  is  linear  space-time  adap¬ 
tive  processing  with  sample  support  near  (or  below)  the  space-time  product 
MN.  The  AV  algorithm  is  a  non-invasive  procedure  where  no  explicit  matrix 
inversion/eigen-decomposition/diagonalization  is  attempted.  Mathematically, 
it  creates  an  infinite  sequence  of  filters  that  begins  from  the  target  vector  and 
converges  to  the  ideal  matched  filter.  The  development  of  the  iterative  algorithm 
is  founded  solely  on  statistical  signal  processing  principles. 

The  practical  motivation  behind  the  development  of  the  AV  algorithm  is 
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adaptive  signal  processing  where  the  input  autocorrelation  matrix  is  assumed 
unknown  and  it  is  sample-average  estimated  by  a  given  data  record.  When  the 
autocorrelation  matrix  is  substituted  by  its  estimate  in  the  recursively  gener¬ 
ated  sequence  of  filters,  the  corresponding  filter  estimators  offer  the  means  for 
effective  control  over  the  filter  estimator  bias  versus  (co-) variance  trade-off  [9]. 
Starting  from  the  zero-variance,  high-bias  (for  non-white  inputs)  target  vector, 
we  can  go  all  the  way  up  to  the  unbiased,  yet  high- variance  for  small  data  record 
sizes,  AMF  estimate  and  anywhere  in  between.  As  a  result,  adaptive  filters  from 
this  developed  class  are  seen  to  outperform  in  probability  of  detection  (for  any 
given  false  alarm)  all  known  and  tested  adaptive  detection  means  (including 
for  example  the  multistage  Wiener  filter  algorithm  [10],  [11]  which  is  known  [7] 
to  be  equivalent  to  an  orthogonalized  version  of  the  AV  algorithm  with  vector- 
optimum  AV  weights12).  It  is  worth  mentioning  that  the  familiar  trial-and-error 
tuning  to  problem  and  clata-record-size  specifics  of  the  real-valued  LMS  gain  or 
RLS  inverse  matrix  initialization  constant  or  SMI  diagonal  loading  parameter 
that  plagues  field  practitioners  is  now  replaced  by  an  automated  data-based 
integer  choice  of  one  of  the  recursively  generated  filters  [8].  Numerical  and 
simulation  results  included  in  this  work  support  the  theoretical  arguments  and 
promote  the  new  method  as  the  processor  of  choice  when  only  small  secondary 
data  training  sets  are  available. 

A.  Signal  Model  and  Background 

We  consider  a  narrowband  uniform  linear  array  radar  with  M  antenna  elements 
(subarrays  or  spatial  channels).  We  assume  that  each  element  collects  the  com¬ 
plex  (I/Q)  return  of  a  series  of  N  coherent  pulses  for  some  given  range  cell 
k  =  1, . . . ,  Kmax  =  TpRi/T,  where  Tppj  is  the  pulse  repetition  interval  and  T 
is  the  pulse  duration.  We  organize  the  received  data  in  the  form  of  a  matrix 
Xmxjv,  where  X(m ,  n),  m  =  1, . . . ,  M,  n  =  1, . . . ,  N,  denotes  the  m-element, 
n-pulse  signal  sample.  The  objective  is  to  cope  with  system  and  surrounding 
disturbances  and  detect  in  Xmxw  the  presence  of  a  desired  signal  of  unknown 
amplitude.  Without  loss  of  generality  and  for  notational  simplicity,  we  consider 
a  “vectorized”  form  of  Xmxn,  where  Vec (XmxJv)  =  xmjvxi  is  constructed  by 
sequencing  all  matrix  columns  in  the  form  of  a  vector. 

We  begin  by  casting  the  detection  problem  in  the  context  of  binary  hypoth¬ 
esis  testing.  We  denote  the  disturbance  only  hypothesis  by  H0  and  the  target 
plus  disturbance  hypothesis  by  Hi. 

Hn  :  x  =  i  +  c  +  n 

J  (154) 

H i  :  x  =  qv  +  j  +  c  +  n. 

In  (154),  j  represents  a  mixture  of  L  broadband  directional  interferers  (jammers) 

12The  orthogonal  version  of  the  AV  algorithm  with  conditionally  optimal  AV  weights  as 
presented  in  [12],  [13]  can  still  outperform  the  multistage  Wiener  filter  [10],  [11]  under  small 
sample  support. 
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in  the  far  field,  where  j  =  Vec(J mxn)  with 


E'o  /  i  \  Old 

Ji{n)  e^m~ i)— 


1=1 


n  =  1 , . . . ,  N,  m  =  1 , . . . ,  M. 


(155) 


We  assume  that  Ji(n),l  =  1 ,L,  is  complex  white  Gaussian  distributed  to 
account  for  channel  fading  at  the  pulse-repetition  frequency  as  in  the  model  of 
[14],  [15].  The  antenna  element  spacing  is  d  and  the  radar  carrier  wavelength  is 
A.  The  direction  of  arrival  (DOA)  Oi,l  =  1 ,L,  is  assumed  to  be  uniformly 
distributed  in  [— n/2, 7r/2].  We  find  it  convenient  to  define  the  spatial  frequency 

fi  =  Slnfid,  l  =  1  and  assume  that  fi  is  uniformly  distributed  in  [- 

0.5, 0.5]  with  proper  selection  of  d  and  A.  In  addition,  Cmjvxi  in  (154)  accounts 
for  colored  Gaussian  noise  with  covariance  matrix  Rc  and  corresponds  to  a 
radar  clutter  region.  Spatially  and  temporally  white  disturbances  are  denoted 
by  n.  The  signal  or  “steering  vector”  of  interest,  v,  is  present  in  x  under 
hypothesis  Hi  only.  Without  loss  of  generality  we  assume  that  vffv  =  1  (H 
denotes  the  Hermitian  operator)  such  that  all  energy  signal  characteristics  are 
absorbed  in  the  unknown  complex  amplitude  constant  a.  For  completeness,  if 
v  =  Vec (Vmxjv)  then 


V  (to,  n) 


1  ei27r(m-l)2iii^+j27r(n— 1)^1^- 

Vmn 

n  =  1 , . . . ,  N,  m  =  1, ,  M. 


(156) 


In  (156),  9S  €  [— 7t/2,  7t/2]  is  the  angle  of  arrival  of  the  signal  (target)  of  interest, 
v  is  the  target  radial  velocity,  and  /pr  is  the  radar  pulse  repetition  frequency. 
Once  again,  it  is  convenient  to  define  the  spatial  target  frequency  fs  =  Sln^sd 

and  the  “normalized  Doppler”  target  frequency  fo  =  aJFr  '  case  we 

assume  that  /s,/d  €  [—0.5, 0.5]. 

Given  a  data  vector  x  that  corresponds  to  some  range  cell  k  €  {1, . . . ,  A"max}, 
the  objective  is  to  decide  in  favor  of  H o  or  H\  in  a  way  that  maximizes  the 
probability  of  detection  Pd  =  Prjffi  decided|l?i  true}  subject  to  a  given  false 
alarm  constraint  Pfa  =  Pr{Pi  decided) //q  true}  <  p.  The  optimum  decision 
rule  is  well  known  [16],  [17]  and  of  the  form 

|wHx|  ^  r  (157) 

Ho 

where  r  >  0  is  the  threshold  parameter  to  be  determined  according  to  the 
condition  Pfa  =  P  and  w  is  the  linear  filter  defined  by 

w  =  6R_1v  (158) 

where  b  is  an  arbitrary  positive  scalar.  In  (158),  R  =  Ai/0{xxH}  where  E h,  { ■ } 
denotes  the  statistical  expectation  operation  under  Hi,i  =  0,1.  Then  the  vari¬ 
ance  under  Hq  of  the  test  statistic  wHx  is  Varpf0{wHx}  =  vffR_1v  which 
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implies  that  the  modified  test  statistic 


V  R  ~ 1  X 

Vv"R-!v 


is  zero-mean  unit-variance 


complex  Gaussian  distributed.  As  a  consequence 


v  R 


is  chi-square  dis- 


|  v/v-ffR_1v  | 

tributed  with  two  degrees  of  freedom  and  leads  to  the  familiar  CFAR  (constant 
false  alarm  rate)  optimum  decision  rule 


th 


R 


rH 


R 


‘x|2  § 


(A  >  0). 


(159) 


Substitution  in  (159)  of  the  sample- aver  age  covariance  matrix  estimate  R(AT )  = 
~k  zCfcl i  ^k^-k  from  K  data  samples  from  H0  defines  the  so  called  CFAR  adap¬ 
tive  matched  filter  (AMF)  detector  [2].  In  contrast,  we  recall  that  the  generalized 
likelihood  ratio  (GLR)  test  of  [3]  and  [4]  is 


|vff[R(Af)]~1x|2 

[R(iF)]-1v(  1  +  [R(if)]-ix)  Ho1 


(7  >  0). 


(160) 


In  the  following  section  we  revisit  the  auxiliary- vector  (AV)  algorithm  [7], 
[8]  to  develop  new  decision  statistics  that  maintain  the  principle  of  linear  fil¬ 
tering  for  signal  detection  in  Gaussian  disturbance.  Compared  to  AMF,  GLR, 
multistage  Wiener  filter  and  other  approaches  we  achieve  superior  adaptation 
performance,  especially  when  we  operate  with  small  secondary  data  sets. 


B.  Auxiliary- Vector  Detection 

The  AV  algorithm  generates  an  infinite  sequence  of  filters  {w„})(l0  that  is 
initialized  at  the  steering  vector  of  interest 

w0  =  v.  (161) 

At  each  step  k  +  1  of  the  algorithm,  k  =  0,1,2,...,  we  incorporate  in  w^,  an 
“auxiliary”  vector  component  gfc+i  that  is  orthogonal  to  v  (but  not  necessarily 
orthogonal  to  previously  generated  auxiliary  vectors)  and  weighted  by  a  scalar 
Hk- |—i  and  we  form  the  next  filter  in  the  sequence 

wfc+i  =  wfc  -  fik+igk+i-  (162) 

The  auxiliary  vector  gfc+i  is  chosen  to  maximize,  under  fixed  norm,  the  mag¬ 
nitude  of  the  cross-correlation  between  its  output  gjfVx  and  the  previous  filter 
output  w(^x  and  is  given  by 

gfc+i  =  Rwfc  -  (vffRwi)v.  (163) 


The  scalar  /J>k+i  is  selected  such  that  it  minimizes  the  output  variance  of  the  filter 
Wfc+i  or  equivalently  minimizes  the  MS  error  between  w^x  and  fi%+1  g£_1x.  The 
MS-optimum  j-ik+i  is 


Mfc+i  — 


gf+iRwfc 
gf+iRgfc+i ' 


(164) 
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The  AV  filter  recursion  is  completely  defined  by  (161)-(164).  Theoretical  anal¬ 
ysis  of  the  AV  algorithm  was  pursued  in  [7].  The  results  are  summarized  below 
in  the  form  of  a  theorem. 


Theorem  1  Let  R  be  a  Hermitian  positive  definite  matrix.  Consider  the  iter¬ 
ative  algorithm  of  eqs.  (161)-(16f). 

(i)  Successive  auxiliary  vectors  generated  through  (162)-(16f)  are  orthogonal: 
gf  gi+i  =  0,  i  =  1, 2, 3, . . .  ,  (but  in  general  gf  gj  ±  0  for  \i  -  j\  ±  1). 

(ii)  The  generated  sequence  of  auxiliary-vector  weights  {pn},  n  =  1,2,...,  is 
real-valued,  positive,  and  bounded:  0  <  -r-^—  <  un  <  -v-^— ,  n=  1,2, ... ,  where 
Amax  and  Am;n  are  the  maximum  and  minimum,  correspondingly,  eigenvalues  of 
R. 

(iii)  The  sequence  of  auxiliary  vectors  {g„},  n  =  1,2,...,  converges  to  the  0 
vector:  lim  gn  =  0. 

n— ►  oo 

(iv)  The  sequence  of  auxiliary-vector  filters  {w„} ,  n  =  1,2,...,  converges  to 

the  optimum  filter:  lim  w„  =  L  !  •  □ 

n— >oo  v  rt  v 


If  R  is  unknown  and  sample-average  estimated  from  a  data  record  of  K 
samples,  then  Theorem  1  shows  that 


w„(A") 


Woo  {K)  = 


R  (K) 


-1 


jH 


R  (K) 


(165) 


where  w00(A')  is  the  widely  used  filter  estimator  known  as  the  sample-matrix- 
inversion  (SMI)  filter  [1],  The  output  sequence  begins  from  Wo(A")  =  v,  which 
is  a  0- variance,  fixed- valued,  estimator  that  may  be  severely  biased  (wo(AT)  = 
v  ^  unless  R  =  a2 1  for  some  a  >  0.  In  the  latter  trivial  case,  w0(A')  is 

already  the  perfect  filter.  Otherwise,  the  next  filter  estimator  in  the  sequence, 
wi(A"),  has  significantly  reduced  bias  due  to  the  optimization  procedure  em¬ 
ployed  at  the  expense  of  non-zero  estimator  (co-)variance.  As  we  move  up 
in  the  sequence  of  filter  estimators  w n(K),  n  =  0, 1,2, . . .,  the  bias  decreases 
rapidly  to  zero  while  the  variance  rises  slowly  to  the  SMI  (w00(A’))  levels  (cf. 
'■  165)). 

For  the  proposed  AV  sequence  of  filters  wn  and  any  number  of  auxiliary  vec¬ 


tors  n  =  0,1,2,..,  ,  the  decision  statistics 
complex  Gaussian.  Then, 


\  'w'/  H  vv, 


W^Rw„  H0  ^ 


are  0-mean,  unit- variance 


(166) 


defines  a  CFAR  test  in  parallel  to  the  optimum  matched-filter  test  in  (159).  In 
adaptive  implementations  of  w„  where  the  covariance  matrix  R  is  estimated 
from  a  data  record  of  K  samples,  the  adaptive  auxiliary-vector  test  becomes 


|a  $  MK). 
wn(K)HHwn(K)  H0 


(167) 
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Fig.  23:  Pd  vs  SINR  (one-lag  temporal  clutter  correlation  0.4). 


The  test  is  asymptotically  CFAR  (as  AT  — >  oo,  A (K)  — >  A  in  (159)  for  the 
optimum  normalized  matched-filter  test  for  any  given  false  alarm  rate  Pfa)- 
For  finite  training  sets  it  qualifies  as  Cell- Averaging  CFAR  (CA-CFAR),  as  seen 
by  the  test  denominator. 

As  a  brief  summary,  the  sequence  of  calculations  for  the  design  of  the  joint 
S-T  auxiliary-vector  filter  estimator  is  as  follows.  The  data  record  of  AT  joint 
S-T  input  data  vectors  xi,X2, . . .  ,  x^-  is  utilized  to  obtain  the  estimate  R(A') 
of  R.  Then,  the  sequence  of  AV  filter  estimators  wo(A'),  Wi(A'),  W2(K), ...  is 
generated  by  the  recursive  algorithm  of  (161)-(164)  with  the  true  autocorrelation 
matrix  R  substituted  by  the  estimate  R(AT).  A  specific  AV  filter  estimator 
w„0(Aj  from  the  sequence  can  be  selected  by  the  cross-validation  criterion  of 
Qiao  and  Batalama  [8]  and  used  for  detection  according  to  (167). 

C.  Numerical  and  Simulation  Studies 

In  this  section  we  support  and  illustrate  the  preceding  theoretical  developments 
through  two  representative  case-studies  based  on  the  hypothesis  testing  problem 
in  (154).  In  all  cases  we  assume  presence  of  a  mixture  of  3  (three)  broadband 
interferers  with  Jammer-to-Noise- Ratio  JNR  =  -  «  35 dB,  l  =  1,2,3.  The 

jammers  follow  the  model  in  (155)  and  the  spatial  frequency  fi,  l  =  1,2,3,  is 
randomly  drawn  from  the  uniform  [—0.5,  0.5]  distribution.  The  “peak”  clutter 

(colored  Gaussian  noise)  to  noise  ratio  is  fixed  at  CNR  =  — f  =  40 dB.  The 
false  alarm  rate  is  set  at  Pfa  =  -01  and  for  simulation  purposes  threshold  and 
probability  of  detection  estimates  Pd  are  based  on  10,000  samples  from  H0 
and  Hi  respectively.  All  presented  results  are  averages  over  100  independent 
Monte-Carlo  runs  for  arbitrarily  chosen  target  vectors  and  jammers.  Pd  values 
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Fig.  24:  Pjj  vs  SINR  (one-lag  temporal  clutter  correlation  0.6). 


are  plotted  as  a  function  of  the  total  SINR  defined  by 

SINR  =  |  a  |2  vffR_1v.  (168) 

We  assume  a  radar  array  system  with  M  =  5  antenna  elements  (channels)  and 
N  =  12  pulses  per  coherent  processing  interval  (CPI),  the  GLR  test,  and  the 
adaptive  AV  filter  w„  in  (13),  all  for  secondary  data  sets  of  size  K  =  MN  = 
60.  All  studies  include  the  ideal  matched  filter  as  a  reference  point,  and  the 
adaptive  matched  filter,  GLR  test,  multistage  Wiener  filter,  and  AV  filter,  all 
for  secondary  data  sets  of  size  K  =  M N  =  60. 

Figures  27  and  28  present  probability  of  detection  versus  SINR  results  for 
two  different  clutter  scenarios  generated  as  in  [18]  with  normalized  one-lag  tem¬ 
poral  clutter  correlation  0.4  and  0.6,  respectively.  We  note  that  relatively  low 
temporal  clutter  correlation  (low  pulse-repetition  frequency  /pr  and/or  high 
clutter  velocity)  results  to  benign  eigenvalue  distributions.  We  observe  that  in 
both  cases  the  proposed  AV  filter  w„  in  (161)-(164)  outperforms  the  multistage 
Wiener  filter  and  offers  significant  performance  gains  over  the  adaptive  matched 
filter  and  the  GLR  test. 
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IX.  Subspace  Direction  Finding  with  an  Auxiliary- 
Vector  Basis 

This  work  has  been  presented  at  the  2005  SPIE  Defense  and  Security  Sympo¬ 
sium,  Orlando,  FL,  and  has  been  published  at  the  IEEE  Transactions  on  Signal 
Processing,  Feb.  2007. 

Solutions  for  the  classical  direction-of-arrival  (DOA)  estimation  problem 
can  be  broadly  categorized  into  maximum- likelihood- type  (ML)  algorithms  [1], 
which  are  based  on  techniques  for  the  maximization  of  the  probability  density 
function  of  the  received  signal,  and  subspace  algorithms,  which  are  based  on 
the  decomposition  of  the  autocovariance  matrix  of  the  received  signal.  Among 
the  most  successful  and  popular  subspace  algorithms  are  the  MUSIC  [2]  and 
ESPRIT  [3]  procedures.  In  general,  ML-type  algorithms  have  superior  perfor¬ 
mance  compared  to  subspace-based  techniques  when  the  signal-to-noise  (SNR) 
ratio  is  small  or  the  number  of  snapshots  is  small.  Also,  the  performance  of 
subspace-based  estimators  degrades  substantially  in  the  case  of  correlated  sig¬ 
nal  sources  as  compared  to  ML  schemes. 

In  this  work,  we  attempt  to  exploit  the  structure  of  the  received  data  auto¬ 
covariance  matrix  in  a  new  way.  When  K  distinct  signals  in  space  impinge  on 
M  antenna  elements  (K  <  M),  the  input  autocovariance  matrix  consists  of  a 
rank  K  signal  subspace  and  a  rank  M  —  K  noise  subspace.  Using  the  concept  of 
maximum  cross-correlation  auxiliary  vectors  (AV’s)  [4],  we  create  an  extended 
non-eigenvector  signal  subspace  basis  of  rank  K  + 1  (eigen  decomposition  is  not 
carried  out  at  all).  The  extended  signal  subspace  encompasses  the  true  signal 
subspace  of  rank  K  and  the  scanning  vector  dimension  itself.  Then,  the  pro¬ 
posed  DOA  estimation  algorithm  simply  looks  for  the  collapse  of  the  rank  of 
the  extended  signal  subspace  from  K  +  1  to  K  when  the  scanning  vector  falls  in 
the  signal  subspace.  Extensive  simulation  studies  demonstrate  that  significant 
resolution  performance  improvements  may  be  gained  over  MUSIC,  ESPRIT, 
and  ML  schemes  for  both  uncorrelated  and  correlated  sources. 

A.  Signal  model 

Consider  a  uniform  linear  array  (ULA)  with  M  elements  and  let  6  represent 
the  DOA  of  an  impinging  source  whose  array  response  vector  is  given  by 

se  =  [1  ,e-ndsine  ^e-i(M-l)fdsinO^H  (169) 

where  u>  is  the  carrier  frequency,  c  is  the  signal  propagation  speed,  d  is  the  inter¬ 
element  spacing,  and  H  denotes  the  Hermitian  operator.  When  K  narrowband 
signals  (K  <  M )  on  the  same  carrier  impinge  on  the  array  with  distinct  angles 
of  arrival  9 1, 02)  •  •  •  >  6k,  the  received  signal  snapshot  at  time  i,  rq  G  CM,  is  of 
the  form 

Tj  =  Ax,  +  n,,  i=l,...,N,  (170) 

where  A mxk  =  [s^, . . .  ,  s^*.],  x,  G  CK  is  the  zero  mean  composite  source 
signal  vector,  and  rij  G  CM  represents  Gaussian  noise  with  mean  zero  and  au¬ 
tocorrelation  matrix  E{ nnff}  =  cr2I  (£?{•}  denotes  statistical  expectation).  To 
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avoid  mathematical  -and  presentation-  ambiguities,  the  array  response  vectors 
s , . . . ,  sgK  are  considered  to  be  linearly  independent,  i.e.  the  inter-element 
spacing  is  no  greater  than  half  the  carrier  wavelength.  The  autocovariance  ma¬ 
trix  of  the  array  input  is 


R  =  E{vrH}  =  ARsAh  +a2I  (171) 

where 

Rs  =  E{xxh}.  (172) 

The  source  signal  autocovariance  matrix  Rs  is  diagonal  when  the  sources  are  un¬ 
correlated  and  is  non-diagonal  and  non-singular  for  partially  correlated  sources. 
Complete  knowledge  of  R  cannot  be  assumed  usually;  instead,  we  may  use  as 
necessary  the  sample-average  estimated  array  input  autocovariance  matrix  given 

by 

i  N 

*=^12^  (173) 

i=  1 

where  N  is  the  available  observation  data  record  size. 

B.  Auxiliary- vector  extended  signal  subspace  basis 

In  this  section,  we  consider  the  signal  model  presented  above  and  develop  a 
novel  non-eigenvector  basis  that  spans  the  signal  subspace  extended  in  dimen¬ 
sion  by  the  DOA  scanner  vector  sg,  9  £  (—90°,  90°).  Essential  theoretical 
background  for  this  derivation  is  provided  by  the  adaptive  filtering  work  in  [4] . 
Details  follow. 

In  contrast  to  the  work  in  [4],  here  we  define  the  initial  vector  in  our  basis 
calculations  vq(0)  as  follows: 


vo(0)  = 


Rsg 
IIR.se  || ' 


(174) 


When  the  sources  are  uncorrelated  and  it  happens  to  be  8  =  03  for  some  j  £ 

{i . /'!• 


K 

Rs^  =  (£{x2[j]}  M  +  o-2)  S6j  +  ^2  A{x2[fc]}  Sgksg.  s gk.  (175) 

k= 1,  k^j 

From  (175)  we  can  see  that  vo(0j)  is  a  linear  combination  of  the  K  array  response 
vectors  and  lies  in  the  true  signal  subspace  of  dimension  K.  However,  when 

0  7^  j  =  1,  ■  ■  ■ ,  K, 


K 

Rse  =  ^2  E  {x2[fc]}  Sgksg  Sgk  +  a2sg  (176) 

k=  1 

and  Vo(0)  is  a  linear  combination  of  the  K  +  1  array  response  vectors  and  lies 
in  the  extended  signal  subspace  of  dimension  K  +  1  which  encompasses  the  true 
signal  subspace. 

Having  defined  vq(0)  (the  first  vector  in  our  basis  to  be  formed),  we  now  seek 
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a  vector  gi(0)  (“auxiliary  vector”  as  we  call  it)  that  maximizes  the  magnitude 
of  the  statistical  cross-correlation  between  the  vo(0)  processed  received  data, 
Vq(6) r,  and  the  gi(0)  processed  data,  gff(0)r,  subject  to  the  orthonormality 
constraint  g^(0)vo(0)  =  0  and  g^(0)gi(0)  =  1: 


gi(0)  =  argmax  { \E  |v^(6»)rrffgi(6»)}  I } 

si(0) 

ef(9)v0(9)=0,llgi(9)||=l 

=  argmax  {|v"(0)Rgl(0)|}.  (177) 

Si(6) 

gf  (9)v„(9)=0, ||  gi  (8)11  =  1 


This  constraint  optimization  problem  was  first  posed  and  solved  in  [5]  yield¬ 


ing 


gi(0)  = 


(I-  vo(fl)v<f(fl))Rvo(0) 


(178) 


1(1  —  Vo(0)VQf  (0))RvO(0)||  ' 

The  vector  gi(0)  has  unit  norm,  is  orthogonal  to  vo(0),  and  is  a  linear  combi¬ 
nation  of  the  K  source  array  response  vectors  when  9  =  9j,  j  =  1 , . . . ,  K,  or 
K  +  1  array  response  vectors  including  s g,  when  9  ^  9 j. 

To  fill  in  the  proposed  basis,  we  define  the  intermediate  vector 


wi(0)  =  vo(0)  -  gi(9)gi(9)  (179) 

where  the  scalar  gi(9)  is  the  value  that  minimizes  the  output  variance  of  the 
wi(0)  processed  data,  E  j |w^(0)r|2  j  [6],  and  equals 


M  i(0)  = 


gf(g)Rv0(e) 

gf(0)RglW 


(180) 


Then,  recursively,  for  n  =  2, . . . ,  K  —  1  we  optimize  the  basis  as  follows:13 

(i  -  vo(g)v( f(g)  -  En1  giWgfW)  Rw„-,p) 

Sn  ||(I-  vo(0)v^(0)  -  Si(9)sf  (6»))Rw„_i(6»)||  ’ 


(181) 


Hn{9) 


Sn  (8)Rw„_i(8) 
g"(0)Rgn(0) 


(182) 


W n(9)  =  w„_i(9)  -  gn(0)gn(0).  (183) 

The  final  auxiliary  vector  g k  {&)  that  completes  the  construction  of  the  basis  is 
not  normalized: 


g  k{9) 


I  -  vq(9)vq  (9)  -  ^  gi(6>)gf  ( 6 )  Rw/f_i(0) 


(184) 


13 At  each  step,  g n(0)  maximizes  the  cross-correlation  magnitude  between  the  wn_i(0) 
processed  data,  w^_1(0)r,  and  the  gn(0)  processed  data,  g^(0)r.  Then,  Hn(0)  minimizes  the 

variance  E  j  |w^(0)r|2| . 
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When  9  =  dj,  j  =  1, . . . ,  K,  g k(9)  =  0  and 


span{v0(dj),gi(dj),---,gK-i(9j)}  =  span{sBl, . . . ,  seK  }.  (185) 

In  this  case,  the  vectors  (vo(0,),  gi(9j), . . . ,  gK-i(9j)}  form  an  orthonormal 
basis  for  the  true  signal  subspace.  When  9  ^  6j,  j  =  1, . . . ,  K,  however, 

span{v0(e),g1(d), . . .  ,gK(d)}  =  span{se,s01,...,seK}  (186) 


and  including  gk(9)  we  have  formed  an  orthogonal  basis  for  the  extended  signal 
subspace. 

Calculation  of  the  auxiliary  vectors  gi  (9) ,  g2  (9) , . . . ,  g k  {9)  as  presented  above 
in  (178)-(184)  requires  calculation  of  the  intermediate  vectors  w n(9)  and  scalars 
pn{9),  n  =  1, 2, . . . ,  K  —  1.  Proposition  1  below  shows  that  direct  second-order 
recursive  calculation  of  g2(0), . . . ,  gK-i(9)  is  possible  (see  also  [7]-[9]),  while 
calculation  of  the  final  unnormalized  auxiliary  vector  g k(9)  needs  in  addition  a 
(first-order)  recursion  on  pn(9).  14 

Proposition  1:  The  auxiliary  vectors  g2($),  g3  (#),•••  ,gif-i(#)  can  be  calculated 
as  follows: 


g  n(0)  = 


(I  ~  E"=~„1_2  g#)gf(g))  Rg,,-!  W 

|(I-E"rn1_2gIWgfW)Rgn_1W||: 


n  =  2, 3, . . . ,  K  —  1, 


(187) 


where  gi(0)  is  given  by  (178)  and  go  (9)  =  vo(0). 

The  final  (unnormalized)  auxiliary  vector  g k(9)  equals 

gK(0)=  -pk-i(0)[i-  ^  &(0)g? w)  Rgx-iW 


(188) 


where 


Pn(8) 


pn—  1  (9) 


gn(g)Rgn-t(g) 

g«(e)Rgn(0) 


and  p\{9)  is  given  by  (180).  || 


n 


2,3, . . .  ,K  —  1, 


(189) 


C.  DOA  estimation 

Having  defined  the  auxiliary-vector  basis  {vo(0),gi(0), . . .  ,g#(0)},  we  are 
now  ready  to  describe  the  proposed  DOA  estimation  procedure.  Let  A  be  the 
prearranged  angle  search  step  in  degrees  and,  without  loss  of  generality,  say  that 
180°/A°  is  an  integer.  Define  9^  =  nA°,  n  =  1,2,...,  180°/A°,  and 

S (0(n))=  [vo(0(n)),gi(0(n)),...,gx(0(n))]  •  (190) 

14 Alternatively,  the  proposed  basis  (vo($),  gi(#),  •  •  • ,  gAr($)}  can  be  calcu¬ 
lated  by:  (i)  Gram-Schmidt  orthonormalization  of  the  Krylov-type  basis  [10]- [12] 
{Rs0,  R2S0, . . . ,  that  is  modified  to  begin  with  the  first-order  vector  Rs#  or 

(ii)  through  the  first  K  + 1  stages  of  the  multistage  Wiener  filter  representation  [8],  [13]  when 
the  initialization  of  the  first  stage  is  changed  to  Rs^.  In  either  case,  the  last  vector  in  the 
basis  is  to  be  left  unnormalized. 
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To  analyze  the  behavior  of  Pav {9^)  and  its  suitability  for  DOA  estima¬ 
tion,  we  consider  first  the  case  where  the  input  autocovariance  matrix  R,  as 
it  participates  in  g k(9^)  and  S(d(™-1)),  is  perfectly  known.  When  9^  = 
9j,  j  =  g  k(O^)  =  0.  As  0(")  -  Oj,  j  =  1, . . . ,  K,  PAV(0 (n))  -  oo 

and  a  peak  in  the  spectrum  is  to  be  observed.  When,  on  the  other  hand, 
e{n)  ±  o j,  j  =  i,...,k,  g k(o^)  is  a  linear  combination  of  the  K  source 
signal  array  response  vectors  and  s 0(„).  The  columns  of  S(0("~R)  span  the 
source  signal  subspace  extended  by  the  se(„_i)  dimension.  Hence,  a  part  of 
g k(9^)  lies  in  the  subspace  of  S(0('™-1^)  and  g^(9^n'l)S(9^n~1^)  0.  Having 

Q{n-i)  _  Q..  j  =  1, . . .  ,K,  does  not  affect  these  findings. 

In  practical  scenarios  where  R  is  not  known  and  is  instead  estimated  from 
finitely  many  collected  data  (by  (173)  for  example),  the  above  arguments  become 
approximate  in  an  estimation  theoretic  sense.  When  j  =  1 ,K, 

g k(9^)  yf  0  with  probability  one.  We  expect,  however,  that  ||g/c(0^n')) ||  be¬ 
comes  increasingly  small  as  the  data  record  size  increases.  At  the  same  time, 
gic(0^)J-{vo(^rd)i . . . ,  gA'-i(0^)}  by  design  and  g k(9^)  lies  outside  the 
approximate  signal  subspace.  Hence,  g§(9^)S(9^n~^)  ~  0  and  we  expect 

to  see  peaks  in  Pav(9^)  as  9 ^  —>  Oj,  j  =  1, . . . ,  K. 

The  computational  complexity  of  the  algorithm  outlined  above  is  0(M2K ) 
per  test  angle  or  0(^-M2K)  in  total  (notice  that  no  direct  matrix  inversion 
operation  is  required).  As  a  reference  comparison,  MUSIC  and  ESPRIT  cost 
one  full  eigenvector  decomposition  of  0(AI3). 


D.  Simulation  studies 

For  illustration  purposes  we  consider  a  uniform  linear  antenna  array  of  ten 
elements  ( M  =  10)  and  two  narrowband  uncorrelated  binary  phase-shift-keying 
(BPSK)  signals  ( I\  =  2)  received  in  additive  white  Gaussian  noise.  The  inter¬ 
element  spacing  of  the  array  is  set  exactly  equal  to  half  the  center  wavelength 
of  the  signals  (d  =  (j). 

First,  we  fix  the  signal  angles  of  arrival  at  —1°  and  1°  and  the  SNR’s  at  7dB. 
In  Fig.  27,  we  examine  the  proposed  AV  basis  spectrum  when  the  observation 
data  record  size  is  N  =  60  as  compared  with  MUSIC  and  the  conventional 
matched-filter  (MF)  spectrum.  The  AV  spectrum  resolves  the  two  sources;  MF 
and  MUSIC  both  fail. 

In  Fig.  28,  we  ease  somewhat  the  problem  to  3°  of  separation  and  add  to  the 
comparisons  ESPRIT  and  the  widely  successful  and  popular  grid-ML  estimator 
of  [14,  15]  that  intends  to  approximate  the  true  ML  solution  [16,  17]  by  searching 

15The  proposed  use,  herein,  of  the  modified  orthogonal  AV  algorithm  for  DOA  estimation 
via  (191)  is  in  sharp  contrast  to  the  use  in  [13]  of  the  standard  multistage  Wiener  filter 
representation  algorithm  of  [8]  for  conventional  reduced-rank  ML-type  DOA  estimation. 
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over  a  set  of 


selected  grid  points.  Fig.  28(a)  presents  a  probability  of 

resolution16  versus  SNR  study  when  N  =  50.  In  Fig.  28(b),  the  SNR  values 
are  fixed  at  OdB  and  the  probability  of  resolution  is  plotted  against  sample 
support.  The  AV  basis  scheme  significantly  outperforms  grid-ML  at  somewhat 
lower  computational  cost.17  MF,  as  one  might  expect,  fails  to  resolve  at  all  SNR 
values  in  Fig.  28(a).  MUSIC  fails  as  well  when  SNR=0dB  for  all  data  record 
sizes  in  Fig.  28(b). 

Next,  we  repeat  the  study  of  Fig.  28  for  a  scenario  that  involves  random 
Gaussian  correlated  sources  si,  S2  generated  as  follows: 

si  ~  Af( 0,of)  and  S2  =  psi  +  \/l  -  p2  s3  (192) 

where  S3  ~  7V(0,  a 2)  and  p  is  the  correlation  coefficient.  Alongside  standard 
MUSIC  and  ESPRIT,  we  consider  MUSIC  and  ESPRIT  upon  spatial  smooth¬ 
ing  [18].  Interestingly,  Fig.  29  demonstrates  that  for  correlation  values  as  high  as 
p  =  0.7,  the  AV-basis  estimator  continues  to  outperform  grid-ML  in  probability 
of  resolution  (as  well  as  MUSIC  and  ESPRIT  with  or  without  spatial  smooth¬ 
ing).  Fig.  29(a)  plots  the  probability  of  resolution  as  a  function  of  the  signals’ 

SNR—  for  fixed  data  record  size  N  =  50.  Fig.  29(b)  plots  the  probability  of 
resolution  as  a  function  of  N  when  SNR=5dB.  At  correlation  value  p  =  0.9,  the 
AV  and  grid-ML  probability  of  resolution  curves  are  seen  to  intersect  (Fig.  30). 
At  small  sample  support,  AV  estimation  continues  to  significantly  outperform 
grid-ML  (Fig.  30(b)). 

E.  Conclusion 

We  derived  a  new  basis  for  the  extended  signal  subspace  that  does  not  involve 
eigenvectors  and  includes  the  search  vector  dimension.  Borrowing  terminology 
from  related  recent  literature,  we  called  the  elements  of  the  basis  auxiliary  vec¬ 
tors  (AV’s).  In  this  report,  we  used  the  new  AV  basis  for  direction-of-arrival 
estimation  with  remarkable  success  relative  to  (grid-)ML-type  estimation  for 
uncorrelated  or  even  highly  correlated  sources  (the  latter  being  especially  true 
under  small  sample  support  operation).  Intuitively,  one  may  attribute  the  suc¬ 
cess  to  the  specific  basis  in  use  together  with  the  joint  effect  of  rank  collapsing 
on  the  norm  of  the  last  AV  and  its  projection  onto  the  previously  calculated 
scanner-extended  signal  subspace. 


16  Two  sources  with  DOA  0 i  and  62  are  said  to  be  resolved  if  the  respective  estimates  6\ 
and  62  are  such  that  both  \6\  —  6\\  and  |02  —  02 1  are  less  than  \6\  —  ^2 1 / 2  [14]. 

17The  AV  estimator  is  run  with  A  =  0.5°  for  a  computational  cost  of  the  order  of 

O(360M2 K).  Grid-ML  has  complexity  O  M2K^j.  Here  (Fig.  28),  N  =  50,  K  =  2, 

and  (^)  =  1225. 
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(a) 


(b) 


Data  record  size  N 


Fig.  29:  (a)  Probability  of  resolution  versus  SNR  (separation  3°,  N  =  50,  cor¬ 
relation  70%).  (b)  Probability  of  resolution  versus  sample  support  (separation 
3°,  SNRi=SNR2=5dB,  correlation  70%). 
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(a) 


(b) 


Data  record  size  N 


Fig.  30:  (a)  Probability  of  resolution  versus  SNR  (separation  3°,  N  =  50,  cor¬ 
relation  90%).  (b)  Probability  of  resolution  versus  sample  support  (separation 
3°,  SNRi=SNR2=5dB,  correlation  90%). 
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X.  An  8x8  Quasi-Orthogonal  STBC  Form  for 
Transmissions  over  Eight  or  Four  Antennas 

This  work  has  been  presented  at  IEEE  ICASSP  2007,  Honolulu,  Hawaii. 

Orthogonal  space-time  block  codes  (O-STBC)  [1]  achieve  full  transmit  di¬ 
versity  and  allow  single-symbol  (two  real  symbols)  maximum  likelihood  (ML) 
decoding.  The  drawback  of  O-STBCs  is  that  full-rate  codewords  do  not  exist 
for  more  than  two  transmit  antennas.  For  the  case  of  four  transmit  antennas, 
the  rate  limitation  of  O-STBCs  was  overcome  by  quasi-orthogonal  (QO)  STBCs 
at  the  expense  of  diversity  loss  [2]- [4].  Full-rate  full-diversity  quasi-orthogonal 
codewords  for  4-transmit-antennas  were  then  presented  in  [5],  [6]  by  retaining 
the  code  structure  of  [2] ,  [3]  and  modifying  the  constellation  of  some  of  the  sym¬ 
bols.  ML  decoding  of  the  QO-STBCs  in  [5],  [6]  requires  joint  detection  of  two 
symbols  (four  real  symbols).  Interleaving  real  and  imaginary  parts  of  different 
symbols  enables  single-symbol  decoding  of  full-rank,  full-diversity  QO-STBCs 
for  the  4-transmit-antenna  case  [7],  [8],  at  the  expense  of  some  performance  loss 
in  comparison  with  joint  two-symbol  clectection  [5]. 

The  codewords  in  [1]-  [8]  partition  the  symbols  into  orthogonal  sets  and  ML 
detection  requires  only  joint  decoding  of  the  symbols  in  each  orthogonal  set  in¬ 
dividually.  Since  the  complexity  of  the  ML  decoder  increases  exponentially  with 
the  number  of  symbols  in  each  orthogonal  set,  a  trade-off  between  rate/diversity 
and  decoding  requirements  is  taking  shape,  especially  for  large  number  of  trans¬ 
mit  antennas.  In  [9],  QO-STBCs  for  8  transmit  antennas  that  attain  full  di¬ 
versity  and  full  rate  were  presented  that  require,  however,  joint  detection  of 
four  symbols  (eight  real  symbols).  Reduction  in  complexity  was  achieved  for 
the  8-antenna  case  through  the  process  of  interleaving  the  real  and  imaginary 
parts  of  different  symbols  [10],  [11];  the  codewords  can  be  partitioned  into  four 
orthogonal  sets  and  hence  require  joint  two-symbol  decoding  only. 

In  this  paper  we  present  a  new  two-symbol-decodable,  full  rate,  full  diver¬ 
sity  order,  8x8  QO-STBC  form  that  can  be  applied  across  either  4-transmit 
or  8-transmit-antenna  systems.  The  proposed  QO-STBC  employs  constella¬ 
tion  rotation  (CR)  and  symbol-interleaving  similar  to  [10],  [11],  however,  the 
8x8  codeword  can  be  divided  into  two  4x4  codewords  without  the  omis¬ 
sion  of  any  symbols.  We  initially  concentrate  our  efforts  on  the  codeword  for 
the  8-transmit-antenna  system  and  evaluate  conditions  on  the  rotation  angles 
necessary  for  the  codeword  to  achieve  full  diversity  order.  Since  several  rota¬ 
tion  angle  pairs  may  exist  that  maximize  the  diversity  order,  we  investigate 
different  criteria  for  rotation  angle  selection  to  further  improve  error-rate  per¬ 
formance.  Common  choice  for  rotation  angle  optimization  of  QO-STBCs  is  the 
maximization  of  the  diversity  product  which  in  turn  leads  to  minimization  of  the 
pairwise-error-probability  (PEP)-upper-bound  at  (asymptotically)  high  signal- 
to-noise  ratios  (SNR)  [12].  However,  for  space-time  codes  with  large  diversity 
order  and/or  large  number  of  transmit  antennas,  diversity  product  maximiza¬ 
tion  may  not  provide  satisfactory  PEP-bound  minimization  and  error-rate  per¬ 
formance  over  operable  SNRs  [13].  Contrary  to  the  popular  choice  of  diversity 


114 


product  maximization,  we  pursue  other  means  of  optimizing  the  rotation  angles 
to  improve  system  performance.  Within  the  class  of  two-symbol  decodable,  full 
rate  8x8  QO-STBCs  we  examine  four  different  rotation  angle  optimization 
criteria:  (i)  We  find  a  new  expression  for  the  rotation  angles  that  maximize  the 
diversity  product  of  the  suggested  codeword,  (ii)  We  show  that  sum-eigenvalue 
maximization  as  proposed  in  [13]  is  irrelevant/non-applicable  to  the  8-transmit- 
antenna  QO-STBCs.  (iii)  We  suggest,  instead,  and  solve  minimum-eigenvalue 
maximization,  (iv)  Finally,  we  use  directly  the  PEP-upper-bound  to  obtain  new 
true  PEP-upper-bound  optimal  rotation  angles. 

The  rest  of  the  paper  deals  with  codeword  design  and  rotation  angle  opti¬ 
mization  for  the  4-transmit-antenna  system.  For  the  4-transmit-antenna  sys¬ 
tem  we  allow  the  channel  coefficients  of  the  first  4x4  codeword  block  to  be 
correlated  with  the  channel  coefficients  of  the  next  4x4  codeword  block.  We 
show  that  as  long  as  the  correlation  is  less  than  100%,  a  diversity  order  of  8 
is  achieved  using  only  4  antennas.  Using  the  PEP-upper-bound  results  for  fast 
fading  time-correlated  channels  in  [14],  we  re-evaluate  the  PEP-upper-bound  to 
incorporate  the  channel  correlation.  We  then  optimize  the  rotation  angles  using 
the  four  criteria  discussed  above  for  the  8-transmit-antenna  case.  We  find  that 
the  rotation  angles  that  maximize  the  diversity  product  for  the  uncorrelated 
8-transmit-antenna  system  also  maximize  the  diversity  product  for  the  corre¬ 
lated  4-transmit-antenna  system.  Arguably,  this  may  be  somewhat  surprising 
because  it  did  not  seem  possible  in  the  past  to  make  diversity  product  rotation 
angle  optimization  independent  of  the  correlation  coefficient  [14].  For  the  other 
three  criteria,  we  show  that  sum-eigenvalue  maximization  as  proposed  in  [13] 
is  again  irrelevant/non-applicable  to  the  4-transmit-antenna  case,  while  the  re¬ 
maining  two  criteria  (minimum-eigenvalue  minimization  and  PEP-upper-bound 
minimization)  yield  complicated  rather  intractable  optimization  equations  as 
they  become  dependent  on  the  correlation  coefficient. 


A.  Code  Structure  and  Transceiver  Model 

Let  Nt  be  the  number  of  transmit  antennas,  Nr  the  number  of  receive  anten¬ 
nas,  and  T  the  number  of  time  slots  over  which  the  code  is  transmitted.  We 
denote  the  number  of  transmitted  symbols  by  K.  The  eight  symbols  au,  k  = 
1, . . . ,  K(=  8),  to  be  transmitted  are  formed  by  mapping  the  incoming  bits  onto 
known  constellations,  e.g.  quadrature-amplitude-modulated  (QAM),  while  their 
corresponding  constellation  rotated  version  a*,,  k  =  1, . . . ,  K(=  8),  is  created  by 

Cirri  —  [&mR  T  fctm/)e  TTl  —  1,  2,  5,  6,  /,  QT 

an  =  («nfi  +ianI)eie,  n  =  3,4, 7,8, 

where  akR  and  ciki  denote  the  real  and  imaginary  part  of  the  symbol  a*,,  respec¬ 
tively,  and  (f) ,  9,  are  the  rotation  angles  to  be  optimized.  The  symbols  au  are 
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interleaved  to  form  Xk,  k  =  1, . . . ,  K{=  8), 


Xl  =  (ilB  +  *057) 

x3  =  a3R  +  ia7i, 

X5  =  «577  +  7017, 
X7  =  a7R  +  ia3I, 


X2  =  0277  +  iaei, 
Xa  =  «477  +  7087, 
Xe  =  0677  +  *027, 
X3  =  087?  +  7047- 


We  now  form/define  Xi  and  X2  as  shown  below 


Xi 

X2 

X3 

X\ 

x5 

x6 

X7 

x8 

-x2 

Xi 

—x4 

^3 

,  X2  = 

—Xq 

x5 

-x8 

x7 

X3 

X\ 

Xi 

X2 

X7 

X8 

x5 

x6 

—X4 

X3 

-X2 

xi  ] 

L  ^8 

-*6 

X5  } 

(194) 


(195) 


Consider  now  availability  of  either  Nt  =  8  or  Nt  =  4  transmit  antennas.  Below, 
we  describe  our  transceiver  model  for  each  case. 


Case  1:  Eight  transmit  antennas 


The  transmitted  codeword  X  is  of  the  form 


X  = 


Xl  04x4 

04x4  X2 


(196) 


and  the  Nt(=  8)  x  Nr  received  signal  matrix  Y  is  given  by 

Y  =  p\/^rXH  +  N  (197) 

where  A  is  the  received  signal  energy  at  each  receive  antenna,  H  is  the  8  x  Nr 
channel  matrix  and  N  is  the  8  x  Nr  noise  matrix.  The  elements  of  H  and  N  are 
modeled  as  independent  and  identically  distributed  complex  Gaussian  random 
variables  of  zero  mean  and  unit  variance  without  loss  of  generality;  p  =  \f2  sat¬ 
isfies  the  energy  constraint  E{||X|j|,}  =  TNt  (-&{•}  is  the  expectation  operator 
and  ||  •  |||,  denotes  Frobenius  norm  of  a  matrix).  From  an  implementation  point 
of  view,  reduction  of  the  large  peak-to-average-power  ratio  (PAPR)  created  by 
the  zeros  during  transmission  of  (196)  is  achieved  by  multiplying  X  by  an  8  x  8 
normalized  Hadamard  matrix  before  transmission.  This  effect  is  reversed  at 
the  receiver  by  multiplying  the  received  signal  matrix  by  the  transpose  of  the 
Hadamard  matrix.  The  net  transceiver  model  still  takes  up  the  form  in  (196) 
and  (197)  by  virtue  of  the  orthogonality  of  the  Hadamard  matrix. 


Case  2:  Four  transmit  antennas 

Assume  now  that  we  have  available  only  Nt  =  A  antennas.  We  transmit  Xi  in 
the  first  four  time  slots  and  X2  in  the  next  four  time  slots.  The  4  x  Nr  received 
signal  matrices  Y 1  and  Y 2  are  given  by 

Yi  =  pj— X1H1  +  Ni,  Y2  =  pj— X2H2  +N2. 

'  V  Nt  ’  1  \  Nt  ~ 
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(198) 


Similar  to  the  8-transmit-antenna  case,  the  elements  of  the  4  x  Nr  matrices, 
Hi,  H2,  Ni  and  N2  are  modeled  as  identically  distributed  complex  Gaussian 
random  variables  (with  zero  mean  and  unit  variance) .  p  =  1  satisfies  the  energy 
constraint  i£{||Xi|||,}  =  ||X2 1|^}  =  TNt.  Individually  within  a  matrix, 

the  entries  are  taken  to  be  independent  from  each  other;  Ni  as  a  whole  is 
independent  of  N2,  Hi  and  H2  are  independent  from  Ni  and  N2,  but  Hi  and 
H2  are  characterized  as  potentially  correlated  with  each  other  via  a  correlation 
coefficient  p.  To  express  the  correlation  between  the  elements  of  the  two  channel 
matrices  we  define  the  extended  channel  matrix  He  ==  [H^  H^] T  and  obtain 


E{hejh%}  = 


I4  pU 
pU  I4 

085 


j  =  l 

otherwise, 


j,  l  =  1,..., Nr 


(199) 


where  heJ,  j  =  1, 
ficient.  In  similar  fashion  to  H 

T 


Nr,  are  the  columns  of  H, 


Ne=  [N 


and  p  is  the  correlation  coef- 
e,  the  extended  noise  matrix  may  be  defined  as 


N2  ] 


E{nejn^}  = 


with 


j  =  h 

otherwise, 


j,l  =  l,...,Nr 


(200) 


where  nej,  j  =  1, . . .  ,Nr,  are  the  columns  of  Ne.  Combining  Yi  and  Y2  in 
(198)  into  a  single  equation  we  have 


Y e  =  P\j  ^XHe  +  Ne  (201) 

where  Ye  =  [Yf  Yj]T  and  X  is  defined  in  (196).  It  is  important  to  note  that 
(198)  represents  a  system  with  time-correlated  channel  coefficients  while  (201) 
represents  a  system  in  which  the  channel  coefficients  are  correlated  across  space 
(channel  coefficients  of  the  four  actual  antennas  correlated  with  the  channel 
coefficients  of  the  four  virtual  antennas).  This  distinction  is  important  when 
obtaining  the  PEP-upper-bound  for  the  4-transmit-antenna  system  later  in  the 
presentation. 

We  now  proceed  with  the  description  of  the  ML  detector  for  both  cases. 
Let  A  represent  the  symbol  constellation  and  Z  denote  the  set  of  \A\K  symbol 
vector  points  (|  •  |  denoting  cardinality  of  a  set)  in  the  complex  X-dimensional 
space.  If  /(•)  represents  the  one-to-one  mapping  of  the  symbol  vector  a  €  AK 
into  X,  for  the  8-antenna  case  the  ML  estimate  of  the  symbol  vector  assuming 
perfect  channel  state  information  at  the  receiver  is 


aMI 


argmin  llY  -  p w^XH||f, 
k£AK,±=f(a)  V  _ 

argmin  Ejii  {  4^hf  ^ ±H±hj  -  2pJ^Re  { yfXh, }  } 

\  J  J 

(202) 
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( H  is  the  conjugate  transpose  operator  and  Re{-}  denotes  real  part)  where 
yj,  h j,  j  =  1, . . .  ,Nr,  are  the  columns  of  Y  and  H  respectively.  Expanding 
±H±  we  have 


±H±  = 


al2  bl2 
bl2  al2 


0 


0 


cl2  dl2 
dl2  cl2 


(203) 


where 

a 

A 

EfcLi  \xk\2, 

b 

A 

Efcfl  2-Re{x££/s+T£y2}, 

c 

A 

d 

A 

Efc=Tf+i  2 J?e(ar^a:fc_|_7S7/2  } ,  K 

(204) 


Substituting  (193)  and  (194)  in  (204)  we  obtain 


a  —  alR  +  &2R  +  R  +  a5I  +  a6I  +  ®7/^+  ®8/i 

b  =  2aina3R  +  2a2Ra4R  +  2a,5iaTi  +  2aeiasi, 

a2  ^ 2  ^ 2  a2  a2  ^ 2  ^ 2  ^2  ( ZUD ) 

c  =  +  ®1J  +  ®2J  +53b~l"  ®4/’ 

d  =  2a5fla7/{  +  2a6_R08fl  +  20.1/037  +  2a2/a4/. 

From  (205)  we  observe  that  inter-symbol-interference  occurs  among  the  follow¬ 
ing  symbol  pairs  (di,  S3),  (a2,  S4),  (S5,  a 7),  and  (de,  as).  In  (202),  since 
multiplying  XHX  by  hj4  and  h j  does  not  create  any  additional  cross-terms  be¬ 
tween  the  symbols  dk,  ML-decoding  of  the  symbol  vector  a  can  be  reduced  to 
jointly  decoding  the  symbol  pairs  {(ai,  S3),  (a2,  S4),  (S5,  S7),  (a6,  as)}  indepen¬ 
dently.  Since  dk  is  obtained  by  rotating  the  constellation  of  Sfe,  k  =  1, . . . ,  8, 
decoding  the  symbols  pairs  {(Si,  S3),  (S2,  S4),  (S5,  S7),  (S6,  Ss)}  independently  is 
ML-optimum.  For  the  4-transmit-antenna  case  we  replace  H  by  He  and  Y  by 
Ye  in  (202)  to  obtain  a  similar  result. 

As  a  concluding  remark,  we  can  show  that  our  suggested  code  structure  in 
(196)  for  the  8-transmit-antenna  case  is  equivalent  to  the  one  in  [10].  However, 
we  do  favor  the  code  in  (196)  for  two  reasons:  (i)  Its  simplicity  in  code  con¬ 
struction  enables  our  theoretical  analysis  on  rotation  angle  optimization  and, 
(ii)  unlike  the  codewords  in  [10],  [11],  the  codeword  in  (196)  may  be  applied  to 
a  4-transmit-antenna  system  across  eight  time-slots  without  omitting  any  sym¬ 
bols.  Having  shown  already  that  the  code  supports  ML-optimal  two-symbol 
decoding,  we  now  evaluate  its  diversity  order  and  seek  the  rotation  angles  that 
improve  error-rate  performance. 


B.  Diversity  Order  Calculation:  Eight  Transmit  Antennas 

The  probability  of  receiving  the  codeword  X  when  X  ^  X  is  transmitted  is 
upper  bounded  by  [12] 


where  R  is  the  rank  of  (X  —  X)  and  A/,  i  =  1 ,...  ,R,  are  the  eigenvalues  of 
(X  —  X)ff  (X  —  X).  The  minimum  value  of  R  obtained  over  all  codeword  pairs 
is  called  the  diversity  order  of  the  system  [12]  and  dictates  the  slope  of  the  error- 
rate  curves  at  asymptotically  high  SNR  values.  To  obtain  the  diversity  order  of 
the  codeword  in  (196)  we  require  the  eigenvalues  A  *  of  (X  —  X)ff(X  —  X).  Set 

AX  =  X  —  X  and  Ax^  =  —  Xfc,  k  =  1, ...  ,8;  then, 


where 


AXff AX  = 


Aal2  Abl2 
Abl2  Aal2 

AcI2  Adl2 
Adl2  AcI2 


Aa  =  _  Ef=i|A^|2, 

A  b  =  J2k=l2Re{AxkAxk+K/2}, 

Ac  ±  Y^K+ iM2, 

A d  =  J2k=Ki2i2MAxtAxk+K/2},  K  =  A. 


(207) 


(208) 


The  block  diagonal  nature  of  AXflAX  in  (207)  allows  us  to  calculate  its  eigen¬ 
values  as  {(Aa  —  Ab),  (Aa  +  Ab),  (Ac  —  Ad),  (Ac  +  Ad)}  which  exist  with 
multiplicity  of  two.  Expanding  and  simplifying,  the  eigenvalues  are  as  shown 
below: 


{  (A«i_r  —  Ad3//)2  +  (A d2R  —  Au^r)2  +  (Ads/  —  Abu)2  +  (A d6/  —  Adg/)2, 
(Adi_R  +  A  d3/j)2  +  (Aa2R  +  Ad///)2  +  (Ads/  +  Adz/)2  +  (. Aaei  +  Adg/)2, 
(Ads//  —  Adz//)2  +  (AbeR  —  AdS//)2  +  (Adi/  —  Ad3/)2  +  (Ad2/  —  Abu)2, 
(Ads//  +  Aa-m)2  +  (Aden  +  Aci8r)2  +  (Abu  +  Ad3/)2  +  (Ad2/  +  Abu)2  }  ■ 

(209) 

As  each  eigenvalue  is  a  summation  of  squares,  the  set  of  minimum  eigenvalues 
over  all  possible  codewords  pairs  is 


{(AaiR  —  A a3//)2,  (A am  +  Aa3n)2,  (Aau  —  Aa3 /)2,  (Aau  +  Aa3/)2}  (210) 

and  represents  the  worst  case  scenario  for  the  upper  bound  in  (206).  As  long 
as  all  eigenvalues  in  (210)  are  non-zero,  the  codeword  X  achieves  the  maximum 
diversity  order  of  8.  We  now  attempt  to  identify  conditions  on  the  rotation 
angles  that  have  to  be  satisfied  to  allow  the  codeword  X  to  achieve  full  diversity 
order. 

It  is  easily  observed  that  all  the  eigenvalues  in  (210)  are  non-zero  if  and  only 
if  the  product  of  the  eigenvalues  in  (210)  is  non-zero,  which  in  turn  is  the  square 
root  of  the  minimum  value  of  the  determinant  of  the  AX^AX.  We  have, 

min  det( AXHAX) 

=  (Adi//  —  Ad3//)2(Adi_R  +  Aa3n)2  (Aau  —  Aa3/)2  (Aau  +  A  d3/)2] 

=  (Aaf/j  -  Ab2R)(Aa2u  -  Ad27)]4 

[(Aautcos((/))  —  Aau  sin((j)))2  —  (A  a3Rcos(0)  —  A  a3i  sin(9))2] 
x  [( Aamsin(<l >)  +  Aaucos(<j)))2  —  (Aa3nsin(Q)  +  A a3icos(0))2] 
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For  fixed  <j>,  9 ,  we  categorize  (211)  into  five  cases  depending  on  the  values  of 
Aaifl,  Aa3R,  Aai/,  Aa3i.  Note  that  we  cannot  have  Aa\R  =  Aa3R  =  Aai/  = 
A  (i3/  =  0. 

Case  1:  Only  one  non-zero  value 

min  det( AXffAX)  =  [(Aai/;cos(0))2(Aai/jsm(</>))2]4.  (212) 

If  any  of  cos{(j)),  sin{(j)),  cos(9),  sin{9)  is  zero,  full  diversity  order  is  not  achieved. 

Case  2:  Two  non-zero  values  both  from  same  symbol 

min  det( AXH  AX)  =  [((Aaiflcos(0)  -  Aai/sm(</>))2)  ,  , 

({Aa1Rsin((j))  +  Aaucos(4>))2)}4. 

Let  A  denote  the  constellation  of  the  symbols  a*,,  k  =  1, . . . ,  K„  A^  and  Ae 
denote  the  constellations  of  symbols  formed  by  rotating  A  by  (f>  and  9 ,  respec¬ 
tively.  If  any  two  symbols  chosen  from  A ^  or  As  have  the  same  real  part  or  the 
same  imaginary  part,  full  diversity  order  cannot  be  achieved. 

Case  3:  Two  non-zero  values  both  either  real  or  imaginary 

min  det(AXH AX.)  =  ([(Aa1Rcos((j)))2  -  (Aa3Rcos(0))2] 

[(Aa1Rsin((j)))2  -  (Aa3Rsin(9))2])4 . 

If  9  =  (f>,  full  diversity  order  cannot  be  achieved. 

Case  4:  Two  non-zero  values,  one  from  real  part  of  one  symbol  and 
other  from  imaginary  part  of  other  symbol 

min  det(AXH AX)  =  ([(AaiRcos((f)))2  -  (Aa3Isin(6))2] 

[(Aa1Rsin((j)))2  -  (Aa3/cos(6»))2])4  . 

If  9  =  </>  ±  7t/2,  full  diversity  order  cannot  be  achieved. 

Case  5:  All  values  are  non-zero 

If  the  square  distance  between  the  real  parts  of  any  two  symbols  from  A 0  is 
equal  to  the  square  distance  between  the  real  parts  of  any  two  symbols  from  As 
or  the  square  distance  between  the  imaginary  parts  of  any  two  symbols  from  A^ 
is  equal  to  the  square  distance  between  the  imaginary  parts  of  any  two  symbols 
from  As ,  full  diversity  order  is  not  achieved. 

For  a  given  constellation,  several  pairs  of  </>,  9  may  exist  that  satisfy  the  ex¬ 
clusion  conditions  listed  above.  In  the  next  section,  we  investigate  four  different 
criteria  for  rotation  angle  optimization. 
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C.  Rotation  Angle  Optimization:  Eight  Transmit  Anten¬ 
nas 

We  now  consider  four  different  criteria  to  optimize  the  rotation  angles  and  for 
each  criterion,  we  show  how  <f> ,  9  can  be  obtained  for  the  codeword  in  (196)  using 
(210).  For  ease  in  reference,  our  findings  (values  of  (f>,9  for  each  criterion)  are 
summarized  in  Table  IV. 


Case  1:  Diversity  product  maximization 

At  high  SNR  values  assuming  full  transmit  diversity  order  and  Nr  =  1,  (206) 
can  be  approximated  by 

__  i  /Ai  \  1  /  A  \ ~Nt 

/Mx-X)<-(nA,j  (-— )  .  (216) 


Worst-case  minimization  of  the  bound  in  (216)  is  equivalent  to  maximization 
of  the  minimum  product  of  the  eigenvalues  (determinant  of  AX^  AX)  over  all 
possible  codeword  pairs,  which  in  turn  is  commonly  represented  by  the  diversity 
product 


c 


1 

2  VA 


min 


x/x 


det  [AXAXh]|V(2T)  . 


(217) 


Diversity  product  maximization  was  used  as  the  rotation  angle  design  criterion 
in  [10],  [11].  For  the  codeword  in  (196),  the  minimum  determinant  of  AX^AX 
over  all  codeword  pairs  is 


min  det( AXHAX)  = 

[(Aaincos((f>)  —  Aausin{(j)))2  —  ( Aa3Rcos(9 )  —  Aa3j  sin(9))2]4 
x  [( Aa\Rsin{(j) )  +  Aaucos(<f))2  —  ( Aa3Rsin{9 )  +  Aa3icos(9))2]4 

(218) 

and  (j>,9  should  be  chosen  to  maximize  (218).  We  evaluate  and  conclude  (see 
Table  IV)  that  the  codeword  in  (196)  (and  the  proposed  codewords  in  [10],  [11]) 
achieve  diversity  product  of  0.1747  with  4-QAM  constellation  and  0.1071  with 
the  non-rectangular  8-QAM  constellation  for  the  8-transmit-antenna  case. 


Case  2:  Minimum  sum-of-eigenvalues  maximization 

In  [13],  it  was  proposed  that  for  high  diversity  order  systems  (number  of  transmit 
antennas  greater  than  four)  the  minimum  trace  of  AXff  AX  is  to  be  maximized 
over  all  codeword  pairs.  From  our  expression  (209),  we  observe  that 

min  (Tr{ AXff  AX))  =  2  [(Aa^  -  Ad3R)2  +  {Ad1R  +  Ad3R)2 

+(Adu  —  Aa3i)2  +  (Adi/  +  Ad3/)2]  .  , 

=  +  Ad2fl  +  Ad27  +  Ad|7] 

=  4lAal  R  +  Aa3R  +  Aai  i  +  A«!/] 

which  is  independent  of  <p,  9. 
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Hence,  for  the  proposed  codeword  and  that  of  [10],  [11],  the  sum-of-eigenvalues 
criterion  in  [13]  is  not  relevant. 


Case  3:  Minimum  eigenvalue  maximization 

If  r  <  Nt  eigenvalues  of  AXflAX  are  significantly  less  than  1,  then  even  for 
large  SNR  values,  1  +  ~  1  and  (206)  is  approximated  by 


Pr(X  -*  X)  < 


1 


V  41V J 


(220) 


The  system  seems  to  lose  diversity;  for  the  codeword  in  (196)  and  the  proposed 
codewords  in  [10],  [11]  the  occurrence  of  the  eigenvalues  in  pairs  causes  loss 
of  diversity  in  steps  of  two.  In  such  circumstances,  it  appears  reasonable  to 
consider  rotation  angle  choices  that  maximize  the  minimum  possible  eigenvalue 
over  all  pairs  of  codewords.  For  our  code  structure  in  (196),  the  rotation  angles 
that  maximize  the  minimum  eigenvalue  are 

(</>,  9)  =  argmax  min  { (Adm  -  Ad3R)2,  (AdiR  +  A a3R)2, 

<t>,e  (221) 

{Adu  —  Ad3 1)2,  ( Adu  +  Aa3i )2}  . 

The  solution  is  listed  in  Table  IV. 


Case  4:  PEP-bound  Minimization 


We  now  show  that  for  all  STBCs  that  employ  CR  and  have  no  more  than  two 
unique  eigenvalues  of  AXff  AX  over  all  possible  codeword  pairs,  maximization  of 
the  diversity  product  is  equivalent  to  minimization  of  the  upper  bound  on  PEP 
for  all  SNRs.  Minimization  of  the  bound  in  (206)  is  equivalent  to  maximizing 

nf=i  {l  +  If  Ai  and  A2  represent  the  2  unique  eigenvalues  of  AXH  AX, 


then 


argmax  II*  1  (l  +  3^7) 


argmax 

argmax 


+  (1  +  w)]9 

1  +  4Vf(Al  +  A2)  + 


(AlA2)J 

(222) 


where  q  denotes  the  multiplicity  of  the  eigenvalues.  Since 


q(  Ai  +  A2)  =  tr(  AXflAX)  =  ||AX|||. 


is  independent  of  rotation  angles18  we  need  to  maximize  only  the  product  of  the 
eigenvalues  to  minimize  the  bound  in  (206).  In  the  case  of  a  single  unique  eigen¬ 
value  (as  in  O-STBCs  for  example),  the  STBC  is  independent  of  the  rotation 
angle. 

18  Constellation  rotation  or  interleaving  does  not  change  the  transmitted  energy  of  the  STBC 
codeword. 
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While  the  case  of  two  unique  eigenvalues  applies  to  the  4x4  QO-STBCs 
proposed  in  [5] ,  [7] ,  [8]  and  their  choice  of  rotation  angle  is  PEP-bound  optimal, 
for  the  8x8  codewords  four  unique  eigenvalues  exist  and  maximizing  the  eigen¬ 
value  (diversity)  product  over  all  possible  codeword  pairs  does  not  necessarily 
minimize  the  maximum  bound  in  (206). 

We  now  directly  find  the  rotation  angles  that  minimize  the  maximum  (worst 
case)  PEP-upper-bound.  Substitution  of  (210)  in  (206)  gives  us  the  worst  case 
scenario  for  all  codeword  pairs.  We  need  to  optimize  (f>1  9  such  that 


(</>,  9)  =  argmax  n£=i  [  f1  + 

a  a  \  r  / 


cj>,6 

A*  =  {(Aaifl  —  A a3F>)2,  (Aam  +  A a3jt)2,  (A au  —  Aa3/)2,  (Aau  +  Aa3/)2}  , 

*  =  1,2, 3, 4. 


(223) 

Suitable  values  for  A  can  be  chosen  such  that  AXi  >  1  for  all  i  =  1,2, 3, 4.  The 
solution  is  listed  in  Table  IV. 

In  the  next  section,  we  aim  to  extend  the  results  that  we  have  obtained  for  the 
8-transmit-antenna  case  to  the  correlated  4-transmit-antenna  system. 


D.  Diversity  Order  Calculation  and  Rotation  Angle  Op¬ 
timization:  Four  Transmit  Antennas 


Due  to  the  correlation  that  exists  between  elements  of  Hi  and  H2,  the  PEP- 
upper-bound  in  (206)  requires  reevaluation.  The  new  PEP-upper-bound  can  be 
obtained  either  by  considering  the  time-correlated  channel  model  in  (198)  or 
the  space-correlated  channel  model  in  (201).  We  begin  with  the  time-correlated 
model  in  (198). 

Defining  Xt  =  [Xf  X%]T  and  AXf  =  Xt  -  Xt  where  Xt  ^  Xt,  the  PEP  for 
a  time-correlated  fast  fading  channel  (channel  coefficients  changing  every  time 
slot)  is  upper-bounded  by  [14] 


Pr(Xt  -  Xt)  < 


2 rtNr  -  1 
rtNr 


-Nr 


—  rtNr 


(224) 


where  rt  is  the  rank  of  the  matrix  (AXt  AX^)  o  Rt  (o  is  the  Hadamard  product 
operator),  Ati,  i  =  1 ,rtl  are  the  non-zero  eigenvalues  of  (AXtAXf)  o  Rf, 
and  Rt  is  the  time-correlated  channel  matrix  to  be  evaluated.  Conforming  with 
the  signal  model  in  [14],  with  T  =  4  and  Nr  =  1  the  received  signal  vector 

yei  =  [yf  y^]T  (yi  and  y2  are  the  first  columns  of  Yi,  Y2  in  (198))  is 


yei  = 


•  • ,  Djvt]h(  +  nf 


(225) 


where  Dj  =  diag{xti}  are  the  2 T  x  2 T  matrices  created  by  the  columns  of  Xt, 
x(j,  i  =  1,2,  •  •  •  ,Nt.  The  2 T  x  1  noise  vector  nt  satisfies  E{ntn^}  =  I2 t- 
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The  2 TNt  x  1  channel  vector  is  ht  =  [h^,  . . . ,  hfNt]T  where  h ti  is  the  vec¬ 
tor  corresponding  to  the  ith ,  i  =  1  transmit  antenna.  If  heij,  i  = 

1, . . . ,  2 Nt,  j  =  1, . . . ,  Nr,  represent  the  elements  of  He  =  [Hf  Hf  ]T,  then  hti  = 

[hei,  1  ^eM  hei  l  heiA  hei+Nt  l  hei+Nt  l  hei+NtA  /lei+jvt.i]  •  The  time-correlated 
channel  matrix  Rt  can  now  be  defined  as 


Kt^Eihubg} 


14  pi  4 
pi  4  I4 


(226) 


where  I4  is  a  4  x  4  all-one  matrix.  Having  evaluated  Rf,  we  calculate 


(AXtAXf )  o  Rf 


AXiAXf 

pAX2AXf 

AX4  O4 

O4  ax2 
axrsaxh 


pAXiAXf 

AX2AXf 


I4  pl4 

'  AXf  04 

pl4  I4 

04  AXf 

(227) 


with  Rs  =  E{hejh*}  = 


14 

pl4 


pl4 

14 


for  j  =  l  from  (199).  Since 


ig( RsAXff  AX)  =  eig( AXRsAXff)  [15],  expanding  RsAXffAX  we  have 


RsAXffAX  = 


Aal2  Afol2 
A6I2  Aal2 
pAaI2  pA&I2 
pA6I2  pAaI2 


pAcI2  pAdI2 
pAdI2  pAcI2 
Acl2  Adl2 
Adl2  AcI2 


(228) 


This  concludes  the  analysis  work  for  the  PEP-upper-bound  for  the  time-correlated 
model  in  (224). 

When  we  consider  the  space-correlated  model  in  (201)  we  obtain  the  PEP- 
upper-bound  as 

(229) 

where  ry  and  Atj  have  the  same  definition  as  before.  For  clarity  and  brevity  in 
presentation,  the  derivation  of  (229)  is  shifted  to  the  Appendix.  Even  though  the 
obtained  expressions  in  (224)  and  (229)  differ,  they  both  are  primarily  functions 
of  A  ti  which  need  to  be  evaluated  for  rotation  angle  optimization.  Due  to  the  loss 
of  the  block  diagonal  structure  (cf.  (228)),  individual  eigenvalues  of  RsAXflAX 
do  not  yield  simplified  expressions  as  was  the  case  for  the  8-transmit-antenna 
system.  However,  the  product  and  sum  of  the  eigenvalues  can  still  be  evaluated. 


Case  1:  Diversity  product  maximization 

We  have 

Nt 

fj  A  u  =  det{  Rs  AX^AX)  =  det(Rs)det(AXH  AX)  (230) 

i=  1 
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and  the  diversity  product  for  the  4-transmit-antenna  case  is  given  by 


c 


TO*nx^x  |  det  [AXAXH]  det[Rs] 


2^Wt 

(l_p2)2/T 

2  VNl 


-mm 


x/x  ^  [AXAX 


H  l|l/(2T) 


(231) 


We  conclude  that  the  rotation  angles  (f> ,  6  that  maximize  the  diversity  prod¬ 
uct  for  the  uncorrrelated  8-transmit-antenna  case  also  maximize  the  diversity 
product  for  the  correlated  4-transmit-antenna  case  (for  any  correlation  coeffi¬ 
cient  value)  and  provide  maximum  diversity  order  of  819.  This  might  seem, 
arguably,  surprising  in  the  context  of  the  findings  in  [14]  where  for  fast  fading 
time-correlated  channels  optimal  diversity-product  rotation  angles  are  a  func¬ 
tion  of  the  correlation  coefficient  p.  We  can  also  observe  that  when  p  =  1,  the 
diversity  product  in  (231)  becomes  zero  and  the  maximum  achievable  rank  of 
RsAXhAX  is  4  (Aa  =  Ac  and  A b  =  Ad)  which  seems  intuitively  satisfying 
as  X2  experiences  the  same  fading  as  Xi.  Hence,  fluctuation  in  the  channel 
conditions  are  actually  beneficial  to  the  proposed  transmission  scheme. 


Case  2:  Minimum  sum-of-eigenvalues  maximization 

Substituting  from  (205), 


Nt  8  8 

=  Tr(RsAXHAX)  =  2(Aa+Ac)  =  2^  Ad2iR+Ad2u  =  2^Ao?H+A  a2u 

i—  1  i=i  i=i 

(232) 

which  is  independent  of  the  rotation  angles  <fi,  6.  Hence,  once  again,  maximiza¬ 
tion  of  the  minimum  sum  of  eigenvalues  (criterion  suggested  in  [13])  is  irrelevant 
to  this  system.  Obtaining  cf>,  0  for  the  other  two  criteria  (minimum  eigenvalue 
maximization  and  direct  PEP-bound  minimization)  proposed  in  Section  IV  re¬ 
quires  knowledge  of  each  individual  eigenvalue  of  RsAXflAX.  Even  though 
the  eigenvalues  can  be  evaluated,  they  yield  complicated  expressions  that  are 
functions  of  p. 

As  a  concluding  remark,  most  important  aspect  of  the  proposed  transmission 
scheme  for  the  4-transmit-antenna  case  is  that  it  allows  for  transmit  diversity 
order  of  8  with  no  additional  physical  hardware  at  the  transmitter  (4  antennas 
only)  and  requires  only  joint  two-symbol  decoding  at  the  receiver.  In  the  bel¬ 
lowing  section,  we  evaluate  the  error-rate  performance  of  the  codeword  X  along 
with  the  various  criteria  for  rotation  angle  optimization  proposed  in  Section  IV. 

E.  Simulation  Studies 

We  divide  the  presentation  of  our  simulations  studies  into  the  4  and  8-transmit- 
antenna  case.  In  all  studies  we  set  Nr  =  1. 

1 ''However  ,  due  to  the  correlation  between  the  corresponding  elements  of  Hi  and  H2.  the 
codeword  X  experiences  a  loss  in  diversity  product  by  a  factor  of  (1  —  p2)2/T. 
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Case  1:  Four  Transmit  Antennas 

To  evaluate  the  performance  of  the  8x8  QO-STBC  transmitted  by  4  anten¬ 
nas  only,  we  compare  its  error-rate  performance  against  the  4x4  QO-STBC 
codeword  in  [5] .  Note  that  to  our  best  knowledge,  for  a  system  with  4  transmit 
antennas  the  code  in  [5]  has  the  best  error-rate  performance  of  all  QO-STBCs 
that  jointly  decode  two  symbols  or  less.  Since  our  proposed  codeword  has  di¬ 
versity  order  of  8  and  rate  1  over  8  time  slots,  we  may  compare  against  two 
different  4x4  QO-STBC  transmission  schemes.  In  one  scenario,  we  transmit 
the  4x4  QO-STBC  codeword  over  the  first  four  time  slots  and  then  repeat 
the  same  codeword  over  the  next  four  time  slots  to  maintain  diversity  order  of 
8  with,  however,  code  rate  1/2  (4  symbols  over  8  time  slots).  In  the  second 
scenario,  we  transmit  a  4  x  4  QO-STBC  codeword  over  the  first  four  time  slots 
and  then  transmit  a  new  4x4  QO-STBC  codeword  over  the  next  four  time  slots 
to  maintain  code  rate  of  1  (8  symbols  over  8  time  slots)  with,  however,  diversity 
order  4. 

The  elements  of  the  channel  matrices,  h\i  j  in  Hi  and  h2i  ,  in  H2,  are  cor¬ 
related  and  generated  as  follows: 

j  ~  CA//0, 1)  and  h2iij  =  phUj  +  \/l  -p2zitj,  i  =  1, . . . ,  Nt,  j  =  1,  (233) 

where  Zij  ~  CM(0, 1)  and  p  is  the  correlation  coefficient  (CM  stands  for  com¬ 
plex  Gaussian  random  variable).  In  Fig.  1,  we  plot  the  block-error-rate  as  a 
function  of  the  received  signal-to-noise-ratio  (SNR)  for  (i)  our  4-antenna  code¬ 
word  transmission  scheme  in  (195),  (ii)  transmission  rate  1/2,  diversity  order  8 
transmission  of  the  4x4  QO-STBC  of  [5],  and  (iii)  transmission  rate  1,  diversity 
order  4,  4  x  4  QO-STBC.  We  consider  three  channel  fading  correlation  scenarios, 
p  =  0,  p  =  0.5,  and  p  =  0.8.  For  rate  1  codewords,  we  select  the  symbols  from  a 
4-QAM  constellation,  while  for  rate  1  /2  codewords  the  symbols  are  chosen  from 
a  16-QAM  constellation  to  maintain  equal  spectral  efficiency  for  all  transmis¬ 
sion  schemes.  Under  diversity  product  minimization,  the  rotation  angles  for  our 
proposed  codeword  are  {4>,6}  =  {37.9,  21.4},  while  for  the  4x4  QO-STBC  code¬ 
word,  the  rotation  angle  is  7r/4  as  calculated  in  [5].  Since  the  4x4  QO-STBC  has 
only  two  unique  eigenvalues,  the  angle  of  7r/4  is  also  PEP-bound  optimal  (not 
the  case,  however,  for  the  proposed  4-transmit-antenna  transmission  scheme). 

From  Fig.  31,  we  observe  a  substantial  gain  in  performance  achieved  by 
transmitting  the  proposed  8x8  codeword  over  four  antennas  when  compared 
with  the  two  other  4x4  QO-STBC  schemes.  The  best  achievable  performance 
from  the  4x4  QO-STBC  codeword  of  [5]  is  the  direct  rate  1,  diversity  order 
4  transmission,  which  the  proposed  transmission  scheme  outperforms  handily 
even  in  the  extreme  scenario  of  80%  channel  correlation.  In  fact,  to  that  respect 
it  may  be  worth  noting  that  our  matrices  Xi  and  X2  may  not  need  to  be 
transmitted  directly  one  after  the  other  in  time.  Separating  the  transmissions 
of  the  two  codeword  blocks  over  several  time-slot  blocks  (block  interleaving) 
may  help  reduce  the  effective  coefficient  value  of  p  when  allowed. 
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Case  2:  Eight  Transmit  Antennas 

We  now  evaluate  the  performance  of  8  x  8  QO-STBC  in  (196)  under  8-element 
antenna  transmission  and  minimum  eigenvalue  CR  optimization  by  (221),  di¬ 
versity  product  CR  optimization  by  (218),  and  the  proposed  direct  maximum 
PEP-bound  CR  optimization  by  (223).  In  Fig.  32,  we  plot  the  block-error-rate 
versus  SNR  when  the  symbols  are  chosen  from  a  4-QAM  constellation.  For  di¬ 
rect  PEP-bound  optimization  of  <f>,  9 ,  we  solve  for  A  to  obtain  a  received  SNR  of 
20dB.  We  observe  that  the  rotation  angles  that  maximize  the  diversity  product 
and  the  rotation  angles  that  minimize  the  maximum  PEP  bound  provide  the 
best  results,  with  the  latter  having  indeed  better  performance.  The  exact  angle 
values  are  shown  in  Table  IV  (along  with  the  resulting  diversity  product  and 
minimum  eigenvalue) . 

In  Fig.  33,  we  repeat  the  studies  of  Fig.  2  for  symbols  chosen  from  an  8- 
QAM  non-rectangular  constellation  [5].  To  obtain  the  PEP-bound  optimal  </, 6 , 
values  we  solve  for  A  that  corresponds  to  received  SNR  of  30dB.  Again,  all 
calculated  values  are  given  in  Table  IV.  For  reference  purposes,  we  include  in 
our  comparisons  the  8x8  single-symbol  decodable  STBC  in  [7].  Since  the  code 
in  [7]  contains  only  two  unique  eigenvalues  the  rotation  angle  of  tan-1  (1/2)  is 
PEP-bound  optimal  for  that  code;  since  its  rate  is  3/4,  we  select  symbols  from  a 
16- QAM  constellation  to  ensure  equal  spectral  efficiency  for  all  codewords  under 
comparison.  Our  PEP-bound  optimized  codeword  in  (196)  offers  a  gain  of  about 
1  dB  over  the  single-symbol  decodable  STBC  in  [7].  The  minimum  eigenvalue 
optimized  version  performs  almost  similarly  well.  As  argued  in  Section  IV. C, 
due  to  decreased  values  of  the  minimum  eigenvalues  as  compared  to  the  4-QAM 
scenario,  the  maximum  diversity  (eigenvalue)  product  optimized  system  seems 
to  lose  diversity  over  the  operable  SNR  range.  Similar  performance  loss  was  also 
observed  in  [10]  when  the  rotation  angles  were  chosen  to  maximize  the  diversity 
product. 

F.  Conclusions 

We  proposed  an  alternative  representation  of  the  8x8  two-symbol  decodable 
quasi-orthogonal  space-time  block  code  (QO-STBC)  that  can  be  used  on  both 
8  and  4-transmit-antenna  systems.  For  the  8-transmit-antenna  case,  we  derived 
three  different  sets  of  rotation  angle  values:  maximum  diversity-product  opti¬ 
mal,  maximum  minimum-eigenvalue  optimal,  and  minimum  maximum-pairwise- 
error-probability-bound  (PEP)  optimal. 

Most  importantly,  the  proposed  codeword  doubles  the  transmit  diversity 
order  of  a  4-transmit-antenna  system.  By  obtaining  new  expressions  for  the 
PEP-upper-bound  for  correlated  channels  we  were  able  to  find  rotation  angles 
that  maximize  the  diversity  product  to  4-antenna  transmission  systems  and  also 
proved  that  the  rotation  angles  are  independent  of  the  correlation  coefficient. 
As  a  by-product,  we  showed  that  for  the  proposed  codeword,  maximization  of 
the  sum  of  eigenvalues  is  an  irrelevant/non-applicable  criterion  (for  both  8  and 
4-antenna  transmissions). 
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Appendix:  Derivation  of  the  PEP-upper-bound  in  (229) 

The  PEP  for  a  fixed  channel  realization  He  is  given  by  [12] 


Pr(X  -  X|He) 


(234) 


where  Q(x)  —  f^°  e  1 2  dy  and  AX  =  X  —  X.  Using  the  approximation 

1  _  A. 

Q(x)  <  ^e~  2  we  have 

Pr(X  ->  X|He)  <  Ie"3iv7l|AXHc|1^.  (235) 


To  evaluate  Pr(X  — >  X)  we  observe  that  the  channel  coefficients  are  indepen¬ 
dent  and  identically  distributed  (i.ixl)  across  receive  antenna  space  (columns  of 
He  are  i.i.d)  and 


PH{Pr(X-^X)}  < 


^H{n^e-^l|AXs"2} 

\  (PH{e_3l^l|AXh¥|2})Wr 


(236) 


where  Ph{  •}  is  the  expectation  operator  over  the  channel  coefficients  and  hej ,  j  = 
1, . . . ,  Nr,  are  the  columns  of  He.  Correlation  in  space  among  the  four  real  an¬ 
tennas  and  the  four  virtual  antennas  is  given  by  (199)  and  denoted  by  Rs. 
Performing  Cholesky  decomposition  on  R„  =  QQff,  we  obtain  h6j  =  Qv  • 
where  Vj  ~  CN( 0,1s).  From  (236)  we  have 

||AXhe.||2  =  h^AXHAXhej.  =  vfQHAXHAXQVj. 


Performing  eigen-decomposition  on  QH  AX.H AXQ,  it  can  be  shown  [12] 

/  D  /•  \  \  Nr 


Pr(X  ->  X)  <  \  (  [] 


1 


1+  Ak 

1  '  4Nt  . 


(237) 


where  A  j  are  the  eigenvalues  of  QH  AXH  AXQ.  Since  eig(QH  AX.H  AXQ)  = 
eig(QQHAXH AX)  =  eig(RsAXH AX)  =  ei5(AXRsAXk)  [15],  thus  we  have 
the  PEP-upper-bound  in  (229). 
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TABLE  IV 
OPTIMAL  ANGLES 


QAM 

Criterion 

(<M) 

Diversity 

Product 

Min. 

Eigenvalue 

4 

Diversity 

Product 

(37.9,21.4) 

0.1747 

0.0093 

4 

Min. 

Eigenvalue 

(30.9,13.3) 

0.1623 

0.0524 

4 

Max. 

PEP  Bound 

(28.5,40) 

0.1352 

0.0112 

8 

Diversity 

Product 

{tan~l{  2)/2, 
tan-1(l/2)/2) 

0.1071 

2.35  x  10"4 

8 

Min. 

Eigenvalue 

(3,  tan~1{  2)) 

0.0732 

0.0022 

8 

Max. 

PEP  Bound 

(7.2,25.1) 

0.0792 

2.7  x  10"4 
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Fig.  31:  Block-error-rate  versus  SNR  (4-transmit-antenna  system) 


Block  Error  Rate 


Fig.  32:  Block-error-rate  versus  SNR  for  4-QAM  constellation  (8-transmit 
antenna  system). 
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Fig.  33:  Block-error-rate  versus  SNR  for  non-rectangular  8-QAM  constellation 
and  16-QAM  for  the  rate  3/4  code  in  [7]  (8-transmit-antenna  system). 


XI.  Code-Division  MAC  in  Wireless  Sensor  Networks  by 
Adaptive  Binary  Signature  Design 

This  work  has  been  presented  at  the  2007  SPIE  Defense  and  Security  Sympo¬ 
sium,  Orlando,  FL. 

Wireless  sensor  networks  (WSN)  are  characterized  by  their  dense  and  large 
node  population,  severe  energy  constraints  and  asymmetric  many-to-one  data 
flows  [1]  — [4]  .  In  terms  of  node  multiplexing  for  efficient  medium  sharing,  the 
distributed,  decentralized  nature  and  objectives  of  sensor  networks  do  not  fa¬ 
vor  time-division  multiple-access  (TDMA)  developments.  Similarly,  frequency- 
division  multiple-access  (FDMA),  would  require  broad  frequency  tuning  capa¬ 
bilities  of  sensor  transceivers  and  may  be  used  for  control  and  data  channel 
separation  only.  On  the  other  hand,  code-division  multiple-access  (CDMA) 
seems  to  serve  the  particular  needs  of  wireless  sensor  networks  well. 

In  this  work,  we  focus  on  CDMA  channel  design  and  allocation  schemes 
which  allow  many  sensor  nodes  to  transmit  simultaneously  over  the  common 
shared  medium  and  possibly  collaborate  with  each  other.  The  problem  of 
CDMA  signature  set  optimization  over  the  real/complex  field  was  studied  in 
[5]  —  [8]  .  In  the  context  of  optimum  design  of  binary  signature  sets,  past  work  in¬ 
cludes  the  recent  developments  in  [9]  — [12]  that  provide  the  first  designs  of  binary 
signature  sets  that  exhibit  minimum  total-squared-correlation  (TSC)  for  almost 
all  signature  lengths  and  set  sizes.  The  performance  of  other  non-minimum  TSC 
binary  signature  sets  was  studied  in  [13]’  [14]  .  In  this  work,  we  deal  with  the 
problem  of  pursuing  dynamic  code  division  MAC  operations  via  fully  adaptive 
design  of  binary  spreading  codes  that  maximize  the  signal-to-interference-plus- 
noise  ratio  (SINR)  at  the  output  of  the  maximum-SINR  receiver.  The  optimal 
code  is  a  function  of  the  disturbance  autocovariance  matrix  and  its  evaluation 
over  the  binary  field  is  NP-hard. 

A.  Code-divison  MAC:  CDMA  Channel  Design  and  Al¬ 
location 

We  consider  code-division  MAC  along  the  following  basic  steps  (Fig.  34).  We 
note  that  our  focus  is  on  CDMA  channel  design  and  allocation;  other  network 
services  such  as  organization,  synchronization  and  routing  are  built  at  a  higher 
layer. 

In  step  1,  node  1  uses  an  arbitrarily  chosen  signature  and  broadcasts  a 
signal  that  indicates  intention  to  transmit.  The  broadcast  signal  is  received 
by  node  2  in  the  presence  of  other  node  interference  and  noise.  Node  2  sends 
an  acknowledgement  to  node  1  indicating  availability  to  receive  (and  possibly 
process  and  forward)  future  messages  from  node  1.  In  this  acknowledgement, 
node  2  submits  an  “optimally”  designed  binary  signature  (sopt).  Subsequent 
communication  from  node  1  to  node  2  (step  3)  is  accomplished  by  means  of 

Sopt' 

Identifying  an  appropriate  signature  at  step  2  allows  subsequent  data  trans¬ 
mission/reception  that  is  impaired  the  least  by  the  present  disturbance  (other 
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node  interference  and  noise).  The  above  basic  one-hop  transmission  is  gener¬ 
alized  later  in  this  paper  to  various  multi-hop  network  configurations  and  the 
performance  of  such  schemes  is  evaluated  in  the  simulations  Section. 

In  this  work  we  are  interested  in  designing  binary  signature/spreading  codes 
that  maximize  the  SINR  of  the  output  of  the  max  SINR  receiver  which  is  an 
NP-hard  optimization  problem.  We  begin  our  developments  by  introducing 
formally  the  signal  model  and  pertinent  notation. 

We  consider  a  multinode  CDMA-type  environment  where  K  nodes  transmit 
asynchronously,  simultaneously  in  frequency  and  time,  over  different,  in  general, 
multipath  fading  channels. 

Assuming  synchronization  with  the  signal  of  the  node  of  interest  k,  k  = 
1, 2, . . . ,  K  upon  carrier  demodulation,  chip  matched-filtering  and  sampling  at 
the  chip  rate  over  a  presumed  multipath  extended  data  bit  period  of  L  +  N  — 
1  chips,  where  L  is  the  signature  length  and  N  is  the  number  of  resolvable 
multipaths,  the  received  data  vector  r fc(m)  £  ^l+n- i  takes  the  following 

general  form, 

rfe(m)  =  \fE~k  bk(m)Hksk  +  zk+ik  +  nk,  m  =  0,1,...  .  (238) 


In  (238),  with  respect  to  fctli  node,  k  =  1,2,..., A",  bk(m )  €  {±1}  is  the 
mth  data  bit,  m  =  0,1,2,.. .,  E k  represents  transmitted  energy  per  bit  period, 
Hfc  £  (£(l+n-i)xl  js  {.jjg  node  k  channel  matrix  of  the  form, 


hk,  l 

0 

0 

hk,  2 

hk,  i 

0 

TT  ^ 

(_L+iv-i)xr 

hk,N 

0 

hk,N- i 

hk,N 

0 

hk,  i 

0 

0 

•  •  hk,N 

(239) 


with  entries  /ifcjra,n  =  1, . . . ,  IV,  considered  as  complex  Gaussian  random  vari¬ 
ables  to  model  fading  phenomena,  sk  £  {±1}L  is  the  signature  vector  to  be  de¬ 
signed.  zk  £  CL+N~1  represents  comprehensively  multiple-access-interference 
(MAI)  to  node  k  by  the  other  K  —  1  nodes,  ifc  £  <CL+N~1  denotes  multi- 
path  induced  inter-symbol- interference  (ISI)  to  node  k  by  its  own  signal,  and 
rifc  is  a  zero-mean  additive  Gaussian  noise  vector  with  autocorrelation  matrix 

Information  bit  detection  of  node  k  is  achieved  via  linear  minimum-mean- 
square-error  (MMSE)  filtering  as  follows  (for  one-shot  detection  the  data  bit 
index  can  be  suppressed), 


bk  =  sgn  {Re{wMMSE,krk}) ,  k=l,2,---,K  (240) 


where  w MMSE,k  =  cRfe  1HfcSfc  £  CL+N  \  Rfc  =  E{ rfcrf},  c  >  0,  Re{  } 
denotes  the  real  part  of  a  complex  number,  and  E{-}  represents  statistical  ex- 
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pectation.  We  recall  that  the  output  SINR  of  the  filter  w MMSE,k  is  given  by 

E 

SINRMMSE,k(Sk )  =  — 

E 

=  Ek  sfHfR^HfcSfc  (241) 

where 

Rfc  —  E  |(zfc  +  ifc  +  rife)  (zfc  +  ifc  +  n^)  |  (242) 

is  the  autocorrelation  matrix  of  the  compound  channel  disturbance.  Since  the 
effect  of  ISI  is  practically  negligible  for  CDMA  communication  systems,  in  the 
rest  of  this  paper,  we  will  approximate  Rfc  by  Rfc  «  E  |(z^  +  n^)  (zk  +  nk)H 

Our  objective  is  to  find  the  signature  sk  £  {±1}L  that  optimizes  (maximizes) 
SIN RMMSE,ki  he. 

S k,opt  =  arg  max  sffHf RjT1Hfcs.  (243) 

sG{±l}i 

In  the  rest  of  this  section  we  present  an  algorithm  that  performs  exactly  this 
optimization  and,  upon  eigenvector  decomposition,  exhibits  linear  complexity 
in  the  signature  length. 

For  notational  simplicity  we  define  the  L  x  L  matrix 

Qfc  =  HfR^Hfc  (244) 

Then,  the  SINR-optimum  sequence  (cf.  (243))  is  given  by 

Sfc  opt  =  arg  max  sHQks.  (245) 

sG{±1}L 

The  optimization  problem  in  (245)  is  equivalent  to 

sk  opt  =  arg  max  sHQkrs  (246) 

se{±i}L 

where  Qkr  denotes  the  real  part  of  the  complex,  in  general,  hermitian  matrix 

Q k  (Q kr  —  Rc{Qfc}). 

If  we  relax,  for  a  moment,  our  constraint  that  the  signature  (sequence)  al¬ 
phabet  is  binary  and  assume,  instead,  that  s  is  real-valued  (s  £  RL)  with  same 
norm  (sHs  =  L ),  then  the  corresponding  optimization  problem  becomes 

sfc  Lt  =  ar9  max  sHQkrs  (247) 

’  y  sGKl,shs —L 

where  the  superscript  (r)  indicates  that  s^opt  is  real- valued.  The  optimization 
in  (247)  is  carried  over  a  hypersphere  in  the  RL  space  of  radius  L  and  centered 
at  the  origin. 


{ 

WMMSE,k  (\/^fc^fcHfcSfc) 

1 

{ 

wMMS£,ic  (Zfc  +  ifc  +  nfc) 
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Let  {qjfc,i,qfc,2,  ■  ■  ■  ,  qfc.z.}  be  the  L  eigenvectors  of  Qkr  and  {Afc>i,  A fej2,  •  ■  •  ,  Afc,z,}, 
the  corresponding  eigenvalues  such  that  A^q  >\kt2  >  •  •  •  >  A k,L-  The  real¬ 
valued  sequence  that  maximizes  the  right-hand-side  of  (247)  is  well  known  and 
equal  to  the  eigenvector  that  corresponds  to  the  maximum  eigenvalue  of  the 
matrix  Qkr,  i.e. 

skLt  =  ar9  sHQfers  =  qM.  (248) 

’  y  sgRI,sHs —L 

Since  {q^q,  q/c.3,  •  •  •  ,  q/t.z,}  are  orthogonal  to  q^i  by  definition,  the  L  —  1  lines 
defined  by 


qk,i+PQk,i,  i  =  2,  -  ■  ■  ,  L,  peR 


(249) 


lie  on  the  plane  that  is  tangent  to  the  searching  hypersphere,  pass  through  the 
real  maximizer  q^q,  and  define  mutually  orthogonal  directions  of  least  decrease 
in  the  cost  function  sHQkrs  ,  s  £  RL  from  the  optimum  point  sk  lipt  =  q^q.  Pro¬ 
jection  of  the  above  least  decrease  lines  onto  the  searching  hypersphere  results 
in  the  slowest  descent  cords  given  by 


qfc,i  +  PQfcp 
Vl  +  P2  ’ 


*  =  2,-  •  •  ,L,  p  £  R 


(250) 


The  slowest  descent  cords  in  (250)  trace  the  searching  hypersphere,  extend  from 
— q k,i  to  q h,u  i  =  2,  •  •  •  ,  L,  and  pass  through  the  optimum  point  s(  opt.  On  these 
cords,  the  function  sHQkrs  ,  s  £  ML,  takes  values  in  [A kpL,  XkpL]  respectively. 
Our  objective  is  to  identify  the  binary  sequences  that  are  closest  in  the  I2  sense 
to  the  above  least  decrease  cords.  Use  of  least-decrease-driven  search  for  finite- 
alphabet  solutions  has  been  considered  in  [15]  to  convert  general  maximum  like¬ 
lihood  estimates  to  maximum  likelihood  decision.  Some  important  observations 
now  follow. 


(i)  It  is  straightforward  to  show  that  for  any  given  i  £  {2,3,- ••  ,  L},  the 
binary  sequences  that  are  closest  in  the  1 2  sense  to  the  cord  2^==^,  p  £ 

R,  can  be  expressed  as  sgn{ q^i  +  pqk,i)-  The  set  of  all  binary  signatures 
of  the  form  sgn( q^p  +  pq^j),  p  £  R,  has  cardinality  L  +  1. 

(ii)  The  I2  distance  measure  of  a  binary  sequence  x  £  {±1}L  and  a  real 
sequence  y  £  RL  is  conventionally  defined  as  Yhi=i  (xi  ~  2/i)2 ,  where  Xi 
and  yi,  i  =  1,  •  ■  •  ,  L,  are  the  binary  and  real  elements  of  the  vectors  x  and 
y  respectively. 

(iii)  The  binary  sequence  x  £  {±1}L  that  is  closest  in  the  1 2  sense  to  a  given 
real  sequence  y  £  RL  has  elements  x i  =  sgn(-yi),  *  =  !,■••  ,  L. 


(iv)  The  binary  sequence  x  £  {±1}L  that  is  closest  to  a  cord,  e.g.  2^==^, 
p  £  R,  i  £  {2,---  ,  L},  is  the  sequence  that  satisfies  minxg^±1ji,  ||x  — 


qt.i+pqfr,/- 

y/^+P2 


1 2,  Vp  £  R,  where  ||  •  H2  denotes  /2-norm. 
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Fig.  34:  Code  division  MAC. 


Our  proposed  binary  signature  search  algorithm  is  outlined  below. 

Algorithm 

For  *  =  2,  •  •  •  ,  P  (P  <  L)  do  : 

1.  Divide  the  range  of  p  into  non-overlapping  intervals  such  that  adjacent 

intervals  have  opposite  signs  in  exactly  one  coordinate.  The  intervals  are  defined 
by  the  points 

Pn  =  Qk,l,n/ Tl  =  1,  *  *  *  ,  A,  (251) 

where  qk,i,m  n  =  1,  •  •  •  ,  L,  is  the  n-th  element  of  the  vector  qfcy.  Then,  arrange 
{pi,  P2,  •  •  ’  .  Pl}  in  ascending  order  {pi,  p'2,  ■■■  ,  p'L},  pj  <  p'2  <  ■  <  p'L. 

Find  the  index,  say  J,  of  the  first  positive  element  in  the  ordered  sequence 
Pd  P2 )  ‘  ‘  ‘  )  Pli  i-e-  Pj-i  <  0  <  Pj>  J  G  {1, 2,  •  •  •  ,  L}. 

2.  Find  the  binary  sequence  that  is  closest  to  the  line  qj^i  +  pqfe.i  for  each 

interval  of  p,  i.e.  for  p  G  (-00,  pi), (pi,  p'2),-  ■  ■  ,(p’L-i,  Pl)>  (Pl.°°)-  There  are 
L  + 1  binary  sequences  in  total  (denoted  as  s°,  s£,  •  •  •  ,  s^)  that  can  be  computed 
recursively, 

s°  =  sgn[  qfc)1]  (252) 

s[+1=s[.-2S°>J+ieJ+i,  Z  =  0,1,---  ,L  —  J  (253) 

si  =  s[+1  -  2 slJ+leJ+i,  l  =  -1,  -2,  •  •  •  ,  1  -  J  (254) 

where  J+l  is  the  ( J  +  Z)th  element  of  s°  and  ej+i  is  the  ( J  +  Z)-unit  vector  in 
Rl. 

3.  Evaluate  SINR( slk)  for  each  binary  signature  s(.  returned  by  Step  2, 
l  =  0,  •  •  •  ,  L,  and  choose  the  binary  signature  that  gives  maximum  SINR.  ■ 

B.  Simulation  Studies 

We  consider  five  basic  network  configurations  (Fig.  35(a)-(e))  to  demonstrate 
the  performance  of  our  proposed  code  division  MAC  scheme. 

In  particular,  Figs.  35(a)  and  35(b)  represent  a  one-hop  and  a  two-hop  one- 
to-one  transmission,  respectively.  Figs.  35(c)  and  35(d)  represent  a  one-hop, 
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Fig.  35:  Sensor  network  configurations. 


two-to-one  and  a  one-hop  one-to-two  transmission  correspondingly,  while  Fig. 
35(e)  describes  a  4-node  network  configuration  where  communication  of  nodes 
1  and  4  is  achieved  via  node  2  and  3  that  simply  decode  and  forward.  In  every 
network  configuration,  intended  (desired)  signals  are  received  in  the  presence  of 
other  node  interference  and  noise. 

We  compare  the  performance  of  the  proposed  signature  optimization  algo¬ 
rithm  against  the  following  benchmarks:  (i)  The  maximum  complex  eigenvalue 
eigenvector  of  Qfc,  denoted  in  the  figures  as  “complex  max-EV”,  which  is  the 
theoretical  maximum  SINR  signature  solution  over  the  complex,  CL ,  field  and 
(ii)  the  maximum  real  eigenvalue  eigenvector  of  Qfc,  denoted  as  “real  max-EV”, 
which  is  the  theoretical  maximum  SINR  signature  solution  over  the  real,  KL, 
field. 

We  consider  the  multipath  DS-CDMA  signal  model  of  the  previous  section 
with  spreading  gain  L  =  16  and  we  assume  the  presence  of  K  =  8  active  nodes. 
Each  node  signal  experiences  N  =  3  independent  paths  and  the  corresponding 
channel  coefficients  are  assumed  to  be  zero-mean  complex  Gaussian  random 
variables  of  equal  energy.  Then,  following  the  notation  of  the  previous  section, 
the  total  average  received  SNR  for  node  k  is, 

A  ^EliWlWI2}  W||hfe||2} 

SNRk  = - 1 - L  = - ^ - L  (255) 

cH  aA 

We  set  SNRi-3  =  8 dB,  SNR.4_q  =  9 dB,  and  SNRj_g  =  10 dB.  We  initialize 
the  signature  set  arbitrarily  and  execute  each  signature  set  design  scheme  se¬ 
quentially  node-after-node  in  what  we  call  a  multinode  adaptation  cycle.  Several 
multinode  adaptation  cycles  are  carried  out  (all  algorithms  are  seen  to  converge 
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in  about  three  cycles).  A  fixed  Walsh-Hadamard  signature  assignment  is  also 
included  in  the  study  to  challenge,  potentially,  the  notion  of  signature  adaptiv¬ 
ity.  As  we  can  see  from  Figs.  36-40,  the  static  (non-adaptive  or  non-controlled) 
Walsh-Hadamard  signature  assignment  exhibits  rather  poor  performance,  as 
expected,  while  the  proposed  adaptive  binary  signature  design  and  assignment 
scheme  performs  very  close  to  the  theoretical  complex-EV  and  real-EV  bench¬ 
marks. 

C.  Conclusions 

In  this  work  we  focused  on  the  problem  of  CDMA  channel  design  and  allocation 
for  code  division  MAC  of  wireless  sensor  networks  operating  in  multipath  fad¬ 
ing  environments.  We  considered  binary  signature  alphabets  and  thus  pursued 
digital  signature  optimization.  We  proposed  and  studied  a  suboptimal  CDMA 
channel  design  and  allocation  algorithm  that,  upon  eigenvector  decomposition, 
exhibits  linear  computational  complexity  with  respect  to  the  signature  length. 
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:  Net  4:  One-hop,  one-to-two  transmission  (L=16,  K=8) 
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Fig.  40:  Net  5:  Four  node  configuration  (two  decode-and-forward  nodes,  L=16, 


